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COMPUTER PROGRAM FOR SUPERSONIC KERNEL- FUNCTION 
FLUTTER ANALYSIS OF THIN LIFTING SURFACES 

By Herbert J. Cunningham 
Langley Research Center 

SUMMARY 

This report describes a computer program (program D2180) that has been prepared 
to implement the analysis described in NASA TN D-6012 lor calculating the aerodynamic 
forces on a class of harmonically oscillating planar lifting surfaces in supersonic poten- 
tial flow. The planforms treated are the delta and modifed-delta (arrowhead) planforms 
with subsonic leading and supersonic trailing edges, and (essentially) pointed tips. The 
resulting aerodynamic forces are applied in a Galerkin modal flutter analysis. The 
required input data are the flow and planform parameters including deflection- mode data, 
modal frequencies, and generalized masses. 

INTRODUCTION 

References 1 and 2 present the analytical background for, and results from, the 
flutter analysis of thin lifting surfaces based on a supersonic kernel function procedure . 
The present report describes the current version of the computer program used to obtain 
such results, and the user should have access to references 1 and 2 for an adequate imder- 
standing of the program . An appendix describes the evaluation of the supersonic kernel 
function that is used. 



•^0 


Si 


SYMBOLS 

weighting factor for the term (denoted by subscript n,m) in series for Ap^ 
semichord length at root or plane of symmetry 
modal- independent damping coefficient 

coefficient of structural, solid-friction damping for mode i 
amplitude of natural mode-shape deflection for mode i 
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Ap.Ap- 

Qij,Qy 

qi,qj 

Qn 

R 

Re( ),Im( 
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w 

X 
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transfer matrix 

elements of integrating matrices for chordwise Integrations 
kernel function of integral equation (1) 
reduced frequency with reference length Pq, k = 
distribution function in lifting pressure series, n = 1, 2, . . . 

Mach number of stream flow 
generalized mass 

lifting pressure, general, and for mode j, respectively 

dimensional and nondime nsional generalized aerodynamic forces 

generalized coordinates of motion for modes i and j 

integrals (n = 1 to 5) (see eqs. (A 7) to (A12)) 

region of integration on wing (eq. (1)) 

real and imaginary parts of ( ) 

value of spanwise coordinate at right-hand wing tip 

time 

velocity of undisturbed stream flow 

Instantaneous downwash at wing surface, positive with z axis 
chordwise coordinate 

nondimensional local section chordwise coordinate, referred to local chord 

X coordinate of points at which integrands are evaluated for numerical 
chordwise integration, where c = 1, 2, . , . 
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local values of chordwise coordinate at leading and trailing edges, respectively 


Xq = X - 4 


y spanwise coordinate 

yo = y - ^ 

y^ y coordinate of span stations at which chordwise integrations are done for 

subsequent spanwise integration, where a = 1, 2, , . , 

Q! mass of air contained in volume 47rbQ3^ a = 4irpbQ^ 

|3 = ^m2 - 1 
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width of subregion HI (see fig. 3 of ref. 1) 
dummy variable for coordinate y 
limits of integration (ref. 1) 
chordwise coordinate (ref. 1) 
dummy variable for coordinate x 
local value of 4 leading edge 
air density 

dummy variable of integration 
circular frequency of oscillation 
chosen base or reference frequency 
natural frequency of mode i 
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PROBLEM DESCRIPTION 


For oscillating and steady thin lifting surfaces in supersonic flow, the linear inte- 
gral equation that relates the distributions of downwash and lifting pressure on such sur- 
faces is (see eq. (1) of ref. 1 and eq. (6) of ref. 2) 


w(x,y,t) 

V 


(2bn)^ fT 

4^^ j]e K(M,k.X(|,yo) d| d. 


( 1 ) 


where the region of integration R is the area of the wing surface bounded by the forward- 
facing Mach cone with its apex at x,y. Generally, the downwash w(x,y,t) is known but 
the lifting pressure distribution Ap(4,tj,t) is not known. A finite series is assumed for 
Ap and each term in the series includes an unknown weighting factor. (See equations (11) 
to (17) of ref. 1.) Equation (1) is solved by collocation or by a least-squares solution 
technique. 

The surface integrations and the singularity extractions indicated in equation (1) 
are performed mainly by numerical quadrature, with some aid from closed-form inte- 
gration. The surface integrals are evaluated and the simultaneous equations are 
solved for the weighting factors and hence the pressure distribution under the control 
of subprogram D2181. 

Subprogram D2182 uses the results from subprogram D2181 plus the input modal 
deflection data to compute the generalized aerodynamic force elements (see eqs. (3), (5), 
and (11) of ref. 2) 


r® Ap^(x,y,t) 

= ^^0 Jo J hi(x,y) dx = 4,7pV2boQ* 

U x^g e 

> ( 2 ) 

Ze J 

and then, for a sequence of input values of an air density parameter, to solve the flutter 
stability equations (see eqs. (A9) to (A13) of ref. 2) 
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(i = 1, 2, . . .) (3) 
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in searching for the value(s) of air density that correspond to neutrally stable motion(s). 
The value that corresponds to the lowest dynamic pressure is on a critical boundary, that 
is, a flutter boundary. Program execution continues for any particular case until a spec- 
ified number of air density parameters is used. Succeeding cases (other planforms, 
modes, k-values, M-values, and so forth) are executed until an end of file is reached. 

PROGRAM ORGANIZATION 
Overlay Structure 

This program is organized in overlays to reduce the required field length. The 
overlay structure is described by the following diagram: 

OVERLAY (JMF,0,0> 

PROGRAM MAIN 


0VERLAYIJMF,1,0) OVERLAYt JMF,2,0> 

PROGRAM D2181 PROGRAM 02182 

MAIN Labelled COMMON 

The following listing contains the FORTRAN variables that appear in labelled 
COMMON: 


COMMON label 

FORTRAN 

variable 

Description 

/TRA2182/ 

XM 

Mach number of free stream 


XK 

reduced frequency k 


TANLE 

tangent of leading-edge sweep angle 


TANLTE 

tangent of trailing- edge sweep angle 


S 

ratio of wing half span to root- 



chord length 


NSYM 

integer quantity: = 0 for spanwise 



antisymmetry, ^ 0 for spanwise 
symmetry of lift distribution 


IDENT 

identification information, up to 



80 Hollerith characters 

/DEFAULT/ 

HMLTPLY 

two-word array, if nonzero, multi- 


plies real and imaginary parts of 
downwash input, respectively 
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D2180 Subprogram Descriptions, Flow Charts, and Listings 

Program D2180 is overlaid in order to reduce the maximum required field length. 
The main overlay 0,0 consists of subprogram MAIN. Primary overlay 1,0 is controlled 
by subprogram D2181, and primary overlay 2,0 is made up of subprogram D2182. There 
are no secondary overlays. 

This section of the report presents a brief description of each subprogram, its 
function, flow chart, and listing. A list of the subprograms and their descriptions 
follows: 


Subprogram 


Description 


MAIN 

D2181 

PARTI 

KERNEL 

MTXMPY 


SSSSS 


GAUSS 

LCOMP 

PRT2 


GLSP 

SIMEQ 


Reads a NAMELIST input parameter, writes a heading, and directs 
the access to overlay levels 1,0 and 2,0. 

Writes headings, reads and writes IDENT input, reads and writes 
NAMELIST input, and calls PARTI and PRT2. 

Computes elements of Rnm matrix (eq. (32) of ref. 1) 

Computes K (KBAR) (eqs. (3) and (4) of ref. 1) 

Computes matrix in (eqs. (28) to (31) and table I of ref, 1), multi- 
plies as in equation (28) of reference 1 to get n^m by row, 
and writes III matrix on TAPE 9. 

Stores span stations at which chordwise integrations are to be 
obtained and also stores the integrating factors for the span- 
wise integrations to obtain elements of . 

For each subregion of surface integration, computes the span sta- 
tions and spanwise Gaussian quadrature factors. 

Computes the (eq, (13c), ref. 1) 

Optionally multiplies the real parts -(ahj^/ax] and the imaginary 

parts -2hi of the downwash input array by HMLTPLY(l) and 

HMLTPLY(2), respectively; multiplies the imaginary part also 

by k; solves equation (32) of reference 1 for the al^2, matrix; 
til 

writes the matrix and the matrix of residuals from the 

illrl /.V 

least-squares solution; optionally punches the a|^ matrix 
elements. 

Solves N X M system of equations (N g M) by least squares 
(for N < M), or as a linear system (for N = M) 

Solves set of linear simultaneous equations 
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Subprogram 


Description 


D2182 Prints heading, reads and writes IDENT input, reads and writes 

NAMELIST input, reads the III and a^^ matrices from 
TAPE9, and computes the of equation (lib) of reference 2; 

optionally writes matrices (2), (3), (4), (5), (6), (7), (8), and (8x) of 
equations (14), (16), and (18) of reference 2; optionally punches 
Qy matrix and seven associated parameters; writes matrix of 
hj^ and span stations, a^, m^j, g.; writes Q^; solves 

equation (A13) of reference 2 for eigenvalues and eigenvectors; 
for each a, it writes a, flutter determinant, eigenvalues, eigen- 
vectors, eigenfrequencies, g values, and stiffness parameters. 

OVERLAY (JMF,0,0); MAIN 

This is the only subprogram in the 0,0 overlay and its function is to read a single- 
parameter NAMELIST input IP2180, to write that parameter in a sentence, then to call 
overlays 1,0 and 2,0 in sequence if KPROG < 1, or to call only overlay 2,0 if KPROG = 1. 
The flow chart for subprogram MAIN follows. 

MAJOR STEP FLOW CHART 



/ 
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40 


30 


OVERLAYS JMF, 0,0) 

A 

1 

* SEgUENCfc ♦ 

♦ FUR MODS ♦ 

iOOOOO 

PROGRAM MAINilNPuT=l,0UTPUT=l,PUNCH=2,TAPty*l,TAPE3=lNPUT, 

A 

U 

200000 

1TAPE6*0UTPUTJ 

A 


300000 

PROGRAM CAKU SPECIFIES MINIMUM BUFFERS TO KtUUCt FIELD LENGTH 


A 

A 

2 

3 

400000 

300300 

DIMENSION IDtNTtSi, HMLTPLYI2J 


A 

A 

4 

5 

600000 

700000 

COMMON /TRA2lti2/AM, XK, TANLE • TANLTE , S ,NMOOE,NSYM,KPKaG 


A 

6 

BOOOOO 

NAMELIST /lP21tiO/KPROG 


A 

7 

900000 

CaHMON/UEFAULT/HMLTPLY 


A 

8 

1000000 

HMLTPLYI LJ*0.,0$hMLTPLY 12) *0^.0 


A 

9 

1 IOOOOO 

READ iD. IP2180) 


A 

10 

1200000 

IF (ENDFILE 3) 20*30 


A 

11 

1300000 

STOP 71 


A 

12 

1400000 

PRINT 40, KPKOG 


A 

13 

1300000 

FORMAT UH1,///31H ♦**6EGlN PROGRAM D2180, NAMELIST 1P2160 IS 

iKPk 

A 

14 

1600000 

10G«f Ilf IH)J 


A 

15 

1700000 

IF (KPROG.EU*!) GO TO 50 


A 

16 

IbOOOOO 

CALL OVERLAY 


A 

IT 

1900000 

CONTINUE 


A 

18 

2000000 

CALL OVERLAY I3HJMFf2,0) 


A 

19 

2100000 

GO TO 10 


A 

20 

2200000 

END 


A 

21- 

2300000 


OVERLAY (JMF,1,0) 

D2181 ,- Subprogram D2181 is the controlling subprogram in overlay 1,0. Its func- 
tion is to direct the calculation of the weighting factors a^^ of the lifting pressure 
series for each mode j of the downwash as described in reference 1. The subprogram 
first calls a NAMELIST input, and then calls PARTI and PRT2. Matrices III and anm 
are written on TAPE 9 for subsequent communication to D2182. The flow chart for sub- 
program D2181 follows. 


HAJOk STEPS IN PROGRAM D21B1 


OVERLAY l«0 \ 

PROGRAM DZiai J 


/PR INI heading/ 

/read, PR I NT INPUT UDENTl/ 

T 

/READ, PRINT NAMELIST INPUT lP218l/ 


fREWiND* TAPE9 | 



\PRT2 / 


/ return N 

^ !TO MAIN) J 


\ 
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UVtKLAyi * 1,0) ^ 

PkGuKAM u2ltil 

,Y( ^8) tlOENT(S) 

OlMtNSIUN 0HXUb,6), HU8,6J, HMLTPLV12) 

COMPLEX XMAT214ci,l6J 
OlMbNSIUiM rtTi'Jb) 

C0M<'4CN/TKA2l62/XMt XK, T A NLE , T A.>iLT Et S ,NHUOfc, NS t h, KPkuo 

CJMMUN/OtEAULT/HMLTPLY 

CUMuHUn ♦UtL ,a, Y,DHX,HtkT 

CUMMUN aA, A b,XC,XN^,XNSt XN6,XN7, XNE , XN9 , BE T ASw , 62 Y 01 f xC , i 0 , XK Ob 2 , 
lXKMb2iXxAY02,X0TtaMtSKAY0,XKYOSg ,XKREAL,XKiMAC,NKKGKi,XMAT2,XMPT 
NAMcU sr / IP2i6l/XK,XM,XAME,TANLT£, S ,NXy, 
lNSYM,JtL,X, r,UHX,H,NMDDE ,WT,HMLTPLV 
rtKlTE IbiiO) 
mRITE 

I-UkMAT I IHO. ///3>X , A5HNA SA - LANGLEY PESEAKcn CENT EK - HAMPTON, VA 
U///) 

FORMAT I ///lOX, oBHJEAN FOSTER FOR HER8 CUNNiNOrtArt- Ucl bO ( I NCLUUESO 
12181 AND D2162)///) 
aRITE lo*30) 

Format uox, /ehonsteady lifting-surfacl theuky oy the suPcksuNic x 
ILKNEL FunCTICn MbTHOD AND A/ , 2 OX , 76hoAL ERKi N moUmL FlJTTLK ANALVSl 
2S Of AR^UrtH£AU PlANFORMS WTTM SUBSCMC LEAul N<// 20a , 3 pHL Oob S ANU SU 
BpERSQnIC TKaIlIn^ EDGES//, 20X,43H«*^REF. 1 = aIaA JUJr^-, NuY^i^bb, 

^ P,iS}ol-l3i>b/,20A, 26H*^»REf • 2 = NASA TN O-o0i2 /2 Ox , t»0h**»REF . 3 B 
NASA TM THE DESCRIPTION CF THIS PKUOkAH 

DEFAULT VALdLS UF OEL.NSYM, WT 
DEL = 0.0AiNSVH=lirt'r*-1.0 
KEAU i Pt^O) iOENT 
FORMAT loAlO) 
rtklTt (bfPU) I dent 
FORMAT l/aAlOJ 
read ip,iP^:lbli 
*vRl TE lo , oO) 

60 FORMAT IZ/iX, j 2H ECHO OF NAMELIST IP2181 FuLLUrtS > 
f^RITE loilPZIbl) 

RtRiND V 

TAPE9 Ib A LCMMuN ICATl ON CHANNEL TC OVERLAY ^,0 

CALL PAKTl 

call FkT2 

RETURN 

END 
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B 

0 
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B 
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b 

b 

fa 

B 

fa 

B 

B 

b 

B 


1 

2500000 

2 

2600000 

3 

2700C00 


ZfaOOOOO 

5 

2900000 

6 

3000000 

7 

3100000 

a 

3200000 

9 

3300000 

10 

3400000 

11 

3500000 

12 

3600000 

13 

3700000 

14 

3800000 

15 

jSQOOOO 

16 

4000000 

17 

4100000 

la 

4200000 

IS 

4300000 

20 

4400000 

21 

4300000 

22 

•♦oOOOUO 

23 

4700000 

24 

4600000 

25 

4S0C00C 

2b 

5000000 

2/ 

5100000 

2B 

5200000 

2S 

5300000 

3C 

5403000 

31 

5500000 

32 

5o00000 

33 

5700000 

34 

5600000 

35 

pSOOGOO 

36 

dOOOOOO 

37 

dIOOOOO 

3b 

□200000 

3S 

O3000CG 

40 

6400000 

41 

6500000 

42- 

6 oOOO 00 


PARTI .- The function of PARTI is to calculate the Ilnm matrix of equation (28) 
of reference 1. By using a number of called subprograms, it first calculates the subdi- 
vision of the region of surface integration as shown in figure 3 of reference 1, determines 
for subregions I, n, and IV whether 6-, 10- , or 16-point Gaussian quadrature is employed, 
carries out the chordwise and spanwise quadrature to calculate the matrix, and 

multiplies the latter by the III matrix of equation (28) of reference 1. The flow chart 
for PARTI follows. 
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COMPUTE THE H SPAN STATIONS AND INTEGRATING FACTORS 
OF EO* 27 FOR SU6REGION III 





Cifo 3 


COMPUTE BOUNDARY 
BETHEEN SUBREGIONS I AND IjJ 


I THERE IS NO SUBREGION II | 



^=T I — ^ 

USE 10-P0INT| I USE 16-POINT 

GAUSS RULE TO COMPUTE SPAN STATIONS 
ANO INTEGRATING FACTORS OF SUBREGION II 



NC 


GAUSS RULE TO COMPUTE SPAN STATIONS AND 
INTEGRATING FACTORS IN SUBREGION I 


CALL SSSSS 
I STORE ABOVE RESULTS! 
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OUMPuTt UNM matrix OF EOS- Id, 26, ANU ^2 AKU uUIkuT iT 
COMPLEX aM2STk(4o> ,XMAT 2(48tl61 
UlMtNSIuN IQLNTld) 

UlMtNSiur^ X< 4b*. Y(48) ,P.NTA8{ U) , TdLlOd 0 ) , ItiLo < o J , 

I TbLlb< lo) ,KTBL6(6) . RTB LI 0 ( 10 J , RT BLl6 il 6 ) , 

2ETKA(N< loJ iKlKANi L6> ,ETASTN( D?> * ^OHSTN i b 0 i » Ko T t id ) , Tn t STN i io J , 
3AL( o) ,EiM;:>PK lo,3i , FNSPI{6,3 ) , FS n 6) ,FSR16 J ,FlPKlo J , 

4FlPUbJ , Aix( 1C J fFF IR( l"'t6 ) ,FF II I lOtCJ rANSRl^oi ,AN5n4Qi , 
tjALMC idJ , TAbRhi lo> ,THEWGH(16J 
UIMENSIUN UHAl40,6},Ht4d,6J 
U 1 MENS luN wT ( J 

C JKMUN/TKA2Ib2/XM,XK, TANLF.TANLT E,S,NM0DE,ivi:>VM,KFKUG 
COMMuN NXY , ObL , Af Y ,DHX tH ,V^T 

o J^MUN AA, Xfci,XC,AN4,XM5,XN6,XN7,XKdTXN9* 6tTA^g,b4:YUl,Ag,YO,XKUo2, 
1 AKMb2,XXAYC2,AOTLRM,SKAY^»XKYOSOi XKREALf XKlMMi^.NKRukVf aMaT2,KMFT 
OATA kNTAu/ 

1 t>.y3330^3CL-02, 4,66287?40E-01, 2 -3783 77 7 Ot-ui , 1 .4 jd229d Ot ♦Ou, 

2 7.9b4 43odCt+C0,-2- 240 82 8srE<*Oi, 7 • 982436oCt ♦ 00 1 1 -4idz ^9o0fc + 00 , 

3 2.37d377/Ot'0i t 4 • 662 87 340E-U1 , 3- 9333 09:>gE-w2/ 
uATA TBLo/ 

I 9 • J 240 9 31 4t“ 01 , 6 - 61 2 09 "*8 6E -0 1 , 2 • 3 b6 1 91oot'”01 # “2 #3ool 91 6Ot-01 1 
2'd-oU093d6L-Cit-9, 32469 514B-0I/ 
data IbLiO/ 

i 9 « ?39 Oo 32 y E^O 1 f 8 • 63063 366F-01 , 6 • 794093i>oL-"u 1 » 4* 3 j-39 3394L" 01 1 

4 i-t8b7433bt-01,-l4 488 74?3aE-01»^4^33393394L“0l t-o. /9^09i>o6E-0I , 
3-d, D30o33o6b-01 ,-9,73906528E-0l/ 

UaTA T8L1o/ 

1 9*d9^0u934t-01, 9.44575Q23E-01 , 8.636312J2L-0i • 7 .334044 Odfc-01 t 

2 6.17b7o244t-01, 4.58316777E-01, 2 .d IdO 3 33»^t-0 1 1 9, 3012 309dE-C2 , 
3-9.3C12 3 09db-J2,-2.816 0^3 5C'E-ai«-4. 5 801677 7i:-01 ,-d - 1 7dV 6244fc-0i i 
4'“ / •3 34 0't4 0dt~01 ,-6.65631 2 02 E~01t -9.4457 50^ 3L- U 1 , ■“ 9. o 94C C93 4t • 0 1 / 

jaTA KTbLd/ 

1 1. n324^92b-0L, 3.60761573E-01, 4 .679 139^4t-Ul , 4 .o 79 13934t-0 1 , 

2 3.6076l373L-Ci, 1 . 7 1324 492 E-0 1/ 

OaTA KTbLlO/ 

1 ij * ot>7 1 3 *<■-►3 b—u 2 , 1 . 49 4 5 1 349E— 0 1 1 2. 1 908 o3u4c-'j 1 , 2.o9boo71 9E— 01 , 

2 2.93324224L-01, 2. 93324224E-J 1 , 2 .69260 / I9c-ul , 2 .i90du3o 2t-01 , 

3 i.494>l-^‘#9b-01, 6.66713443E-02/ 

UATA kToLid/ 

1 2.7l32<»y94L-02, 6.225 35 2396-02, 9 . 5 156 5 Uoc-02 , i . 44D2d97 It- Cl , 

2 1.493939odb-01, 1 . 69 1 56 5 19E*0 1 , 1 . 82603413b -01 , i .6945061 OE-01 , 

3 i. 6943061 OL-Ol , 1 .d260?4l5E-31, 1 .69 1565l9b-3l , 1 .4939 39 tidb-01 , 

4 1 # 2 4o2o9 7 ib — 0 1 1 9. 51 5 05 1 LdE— 02 , 6 - 2253 5239t— 02 , 2 . 7 i3i.45 94b— 02/ 
wRITb (u.lOJ 

FOKMAl ( IHO, 43X, '♦IFi PART I OF 02181 GLNLkaTL UNM MAThlXJ 


print 20 

FORMAT 1 IHO/ 116HD 

K 

TANLit 

tanla 

S 

MACH 

UbLNUM 

NSYMMETRY 


2i 

FKINT 30, XA, * aNLTE, TANLE f StXM,a£L,NSYM 
FUKMAT iotl7-d,U2) 

PRINT 40 

FUKMAT UH0,24H NO. PF C . P TS . ( =NXYi J 

print 50* Nxv 

FOMAT *120) 

aA=-.329 

XB=-1.40c>7 

AC = -2. 9 

XNU-.lOl 

XN2=-. 099 

XN3=-. 09460933 

0LL*P=LlbL 

XM2= AM««2 

BbTA=^SyRT * XM2-1 .0 ) 

Xn2=XK*2.0 

bT ANLE = 6tT Af TANLt 
dbTASu = DtT A**2 
XK0d2=XK2/ bfc TASW 
AAMB2=XM*AK0b2 
00 oO C=i.NAV 

00 oO J = l, Id 


C 2 
C 3 
C 4 
C 5 
C d 
C 7 
C b 
C 9 
6 1C 
C il 
C 12 
C 13 
C 14 
C 15 
C 16 
C 1 7 
C 18 
C 19 
C 2C 
C 21 
C 22 
C 23 
C 24 
L 25 
C 26 
C 27 
C 26 
C 29 
C 30 
C 31 
C 32 
L 33 
C 34 
C 35 
C 36 
C 3/ 
C 38 
C 39 
C 40 
C 41 
C 42 
C 43 
L 44 
L 43 
L 40 
C 47 
C 4b 
C 49 
C 50 
C 51 
C 52 
C 53 
C 54 
C 35 
C 56 
C 57 
C 56 
C 59 
C 60 
C 6l 
C 62 
C 63 
C 64 
C 65 
C dd 
C 67 
C 68 
C 69 
C 70 
C 71 
C 72 
C 73 


obooooo 

6SC0000 
7000000 
7100000 
7^00000 
7300000 
7^^00000 
730uaco 
7O000C0 
7700000 
7800003 
7900000 
dUOQoau 
oiOOOOC 
6200000 
6300000 
64000CO 
630000U 
odOOOCO 
6700000 
dtiOOOCO 
8900000 
9000U0U 
9100000 
9200000 
9300000 
9400000 
9300000 
9o00000 
9700000 
960000G 
9900UU0 
lOOOOOOC 
iOlGUOOO 
10200000 
iOBuOOOO 
lO^tOOOCO 
10300000 
LOoOOOOC 
10700000 
10800000 
10900000 
1 lOOOoOO 
lilOOOOC 
11200 000 
11300000 
11400000 
115OCOO0 
11600000 
1I70000C 
11800000 
11900000 
12 0000 CO 
12100000 
12200000 
12300003 
12^00000 
12300000 
12600000 
12700000 
ItbOOOOO 
12900000 
IjjCOOOCO 
13100000 
13200000 
13300000 
13400000 
13:>000C0 
136OOU0O 
1:^700000 
13800000 
13900000 
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60 

AMAT21L, Ji =(0- »0. J 

c 

7A 

l^OOOOOO 


NxyCNT=l 

c 

75 

14100000 

c 

BEGIN CONTROL-POINT Uu LOuP ♦♦♦»» 

c 

76 

1^200000 

c 

SUBOlVlOb THb REoION OF INTEGRATICN INTO 6UoKt6loNS 

c 

77 

14300000 

c 

AS IN FIG. 3 OF REF., 1. COMPOTE SPAN STATIONS ANt> SPAN 

c 

78 

14400000 

c 

integrating FACTORS. 

c 

79 

14300000 

c 


c 

80 

14600000 


uQ 7b0 hJ=l»NXT 

c 

81 

14700000 


L)U 70 JN = l»^o 

c 

82 

14800000 


ANSKIJN J =0.0 

c 

82 

14900000 

70 

ANSI UN) = 0.C 

c 

84 

13000000 


IC2=1 

c 

83 

13100000 


YTERM=yI MOI 

c 

do 

13200000 


iiVMO=bETA^y I 

c 

87 

13300000 


ETAL=I byMu-XCMO) J/BTANLE 

c 

88 

13400000 


ETAK=( BYMvJ*-aI MUi i /BTANLE 

c 

89 

13300000 


OtLN2=DEL/2.0 

c 

90 

15600000 


aE= 1 tTAk-Y( MGi J/OELN2 

c 

91 

13700000 


it IAt-10.01 60,d0t90 

c 

92 

13800000 

dO 

OtL2=t ETAK-YIMO) )/10,'^ 

c 

93 

13900000 

C 

PKEPAKt SObRtGIQN !V 

c 

94 

1600 0000 


ETA3=LTAK 

c 

93 

16100000 


tTA^sy IMU ) *-01:12 

c 

96 

16200000 


ICGJNT=iO 

c 

97 

16300000 


NPT0Sk=3 

c 

98 

16400000 


CALL GAUSS U 0 # tT AA ♦ ET AS , T 8L 10 ( 1 ) , RTbLl 0 ( U , V i EKH, if ET kANI IJ »RTRAN 

c 

99 

lo30000U 


inn 

c 

100 

16600000 


GU TO 130 

c 

101 

16700000 

yO 

IF (At-20.0i iOOflOOtUO 

c 

102 

16800000 

LOO 

uEL2=UELN2 

c 

1C3 

16900000 


ETAA=YtMw) *-OfcL2 

c 

104 

170000UO 


t T AS =LT A R 

c 

103 

17100000 


GU TO UO 

c 

106 

17200000 

110 

OfcL2=lETAR-Y IMO) )/20-? 

c 

107 

17300000 


tTA^=ylMG) *OtL2 

c 

ICb 

17400000 


t TA>=i:TAK 

c 

109 

17300000 

120 

ICUU NT = lt> 

c 

110 

1 76000CO 


NPTOSR=ti 

c 

111 

L7700000 


call GAUS6 UbftT AAfCTA5*TeLl6a>tRT0Ll6m,yTtRM, it t IRANI DtkTKAN 

c 

112 

17600000 


I ( m 

c 

lU 

17900000 

130 

0SR=bTA3-bTA^ 

c 

114 

IdOOOUUU 


CALL SSSSS ( lC2t ICQUNTtETASTNi ll* WGHSTNI IJ ,hTKAN(lJ fkTkANi lU 

c 

113 

18100000 

C 

PKEPAKc SObREGlCN III 

c 

116 

18^00000 


IC2=ICUUNT*l 

c 

ill 

18300000 


tTA3 = Y IMi^i-DtL2 

c 

118 

18400000 


OTri=ET A4-tTA3 

c 

119 

1 8300000 


W = RlOTH/iG.O 

c 

120 

iboOOOCO 


KG3STI bl = YlMU) 

c 

121 

187OCOO0 


MK=7 

c 

122 

18800000 


Jx = P 

c 

123 

18900000 


TAtiRNio) = ( l.C7UtL2 )*RNTA0 16 I 

c 

124 

19GU00U0 


uO 1 AO LX= It 3 

c 

123 

19100000 


NK = JK*-1 

c 

12b 

19200000 


ML=rtK-l 

c 

127 

19300000 


RGJSTIJKJ =ko3STlNKI*-W 

c 

128 

19400000 


KG3ST4 MK^ = kG3STI ML I^W 

c 

129 

19300000 


TABkNUKJ = ( 1.0/OEL2)«RNTAe<JK> 

c 

130 

IVbOOOOO 


TAbKM Mki = l 1.0/OtL2)«PNTAe IMK) 

c 

131 

1970000U 


JK=JK-1 

c 

132 

19800000 

UO 

f-iX-MXt 1 

c 

133 

19900000 


ICQUNT = 1CuUNT*-11 

c 

134 

2000000U 


CALL iiSSSS 1 IL2, ICQUNTfETASTNI DfWGhSTNlli fKGbSlm rTAbKNiiJi 

c 

133 

201000CO 

C 

PKEPARb SUbKbGICN II 

c 

136 

20200000 


IC2=ICOUNT*-l 

c 

137 

2O300CO0 


OL=tTA3-bTAL 

c 

138 

20400000 


IF U0L/0EL2l“2a.:>) 15'>tl5C,l70 

c 

139 

20300000 

C 

iMbRb IS NO SUBREGION I 

c 

140 

20600000 

130 

ETA2=tTAL 

c 

141 

20700JOO 


J1L = 0 

c 

142 

20800000 


02L=bTA3-ETA2 

c 

143 

20900000 


IF 4 lOL/otL^i-io.r' j isotieotiao 

c 

144 

21000000 
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ibO 


c 

170 

IbO 


190 


C 

200 


213 


220 

230 


240 


2 ij 0 

2t»0 

270 


2 b 0 

290 


300 


310 

C 

c 

c 


320 


CALL UOr fcTA2,ETA3 ,T8L 10( 1 ) t RTBLl 0U> * VT LKM* 1 1 LTk AN ( 1> iKTKAN 

im ) 

lC0UNr=JL0UNT+10 
Ni^TD2L = b 
NPTo1l=0 
OU TQ 190 

bUUNOAKV Ui- SUBREGION I AND JI 
ETA2=t (MOJ-l9*0*O6L2 

CALL ^AUSS < 16 ,fcT A2 f ET A3 tTBLlbU ) f RTdL 16U) t iTEKfl* It t TK AN ( 1 J f KTR AN 
1 m j 

NPTL ) 2 L = ti 
NPTU1l=0 

ICUUNT- UUUNT4'16 

CALL SSSS;> ( 1 C2 1 I COUNT • £ T AST Nil ) . V^iGHSTN < U , tT KAN 1 1 i • KTK AN I in 
lC2=iCUUNT*l 

IE UL)L/OtL2)-2U.0I 26C»26Ct200 
PKLPAKt SUttKtGlON I 
tTAl=tTAL 
01L-ETA3-LTA2 
02L- tT A^-tTAl 

IE MU2L/U1LI--7DI 2l0,210t220 

CALL GAUSS (Qtfcl Al,ETA2tTBL6{lltRT6Lt>I U » V 7tkM , 1 tL TKAN4 iJ tKTKANUl 
I ) 

1 CG0NT=1 CUUNT+o 
NPT02L=d 
NPTU1L=*3 
GO TO 2p0 

lE i ID2 L/l)LL m.3) 230f230t240 

CALL uAUSS 4 10. tTAl ,ETA2.TBL10U ItRTbLlCmt^fEKEj, ItETKAM IJ tKTRAN 

i ( in 

lC0UN7=iCUuNT+10 
NPT J2L=b 
NPTOlL = l) 

GO rO 2:?0 

CALL GAUSS do , fc T A I tE TA2 ,TBL 16 ( 1 ) * RTBL16 (1 n iT tKM, 1, LTRAN ( II t KTKAN 

Kin 

1 COUNT ■= i COUN T-f I u 

NPTU2L-S 

NPT01l=6 

CALL SSSSS lic2. i COUNT .ET AST N( 1 i , V^GHSTNi 1) t tTKAN t il tKTRAN 4 U J 

USL-01L+U2L 

rtRITE 4ot2/uJ 

EQRMAT UH0/9C.H X Y USL 

1 UlL D2L i>SKJ 

WRITE <6t2bOJ XlrtQI tY( MQJ tDSLtUlLtOZL.USR 
EORMaT tscl7-di 
WRITE 4ot290J 

FORMAT l‘t7H0 PTSDIL PTSD2L PTSOSk J 

WRITE 46.3001 NPTDlL.NPTD2LtNPTDSR 
FORMAT 419.2117J 
WRITE 46.310) UtL2 

FORMAT i^lHu HALF-WIDTH OF CONTRCL POINT STRIP E17.6) 

COMPOTE gUANTlTibS INVARIANT WITH CHUROwISE LuUROlNATE 
THElE^O 

THEsTM D^^IhELE 
00 7&0 IU=1*1CU0NT 
IDCAN=10-1 

AETA=AbS4ETASTN( ID)) 

Y0=T IMU J -ETaSTNI I D) 

AYO^ABSI TO) 

psile=aeta*tanle 

PSlMC=Xi M u)-dETA*AVO 
PSIT£^KOfAETA*TANLTE 
PPP=PSlTE+PSiLE 
BB0=PSITE-PS1L£ 

Bb03=- 5*db0 

DUM1 = I PPP-2^ C*PSIMCI /BBO 
IF iOUMl) 320.3d0t333 
nEg=-1 
UO TO 340 


C 143 
C 146 
C 147 
C 146 
C 149 
C 150 
C 151 
C 132 
C 153 
C 154 
C 155 
C L5o 
C 157 
C 156 
C 159 
C 160 
C 161 
C 162 
C 163 
C 164 
C 165 
C l66 
C 167 
C 166 
C 169 
C 170 
C 171 
C 112 
C 173 
C 174 
C 175 
C 17o 
C 177 
C 176 
C 179 
C 160 
C 181 
C 162 
C 163 
C 164 
C 165 
C 166 
C 167 
C 166 
C 169 
C 190 
C 191 
C 192 
C 195 
C 194 
C 195 
C 196 
C 197 
t 196 
C 199 
C 200 
C 201 
C 202 
C 203 
C 204 
C 205 
C 206 
C 207 
C 206 
C 209 
C 210 
C 211 
C 212 
C 213 
C 214 
C 215 
C 216 


21100000 

21200000 

213000D0 

2 W 0 U 000 

21^00000 

2 I 0 OOOOO 

2 U 00000 

21600000 

21900000 

220000CO 

22100000 

22200000 

22300000 

22400000 

225000C0 

22600000 

22700000 

22600000 

22900000 

230000CO 

23100000 

23200000 

23300000 

23400000 

23500000 

23600000 

23700000 

23600000 

23900000 

24000000 

24100000 

24200000 

24300000 

24400000 

24500000 

24600000 

24700000 

24600000 

24900000 

25000000 

25100000 

25200000 

25300000 

25400000 

25500000 

25600000 

25700000 

25600000 

25900000 

26000000 

26100000 

26200000 

26300000 

26-tOOOOO 

26500000 

266000C0 

26700000 

26600000 

26900000 

27000000 

27100000 

27200000 

27300000 

27400000 

27500000 

27600000 

27700000 

27800000 

27900000 

2600C000 

26100000 

26200000 
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330 

NEO=l 

C 

217 

20300000 

340 

AOUrt=AttbiOUMU 

c 

210 

28400000 


II- { AdS( AUUM--1 cO J-o00?001) ?50 • 3 » 300 

c 

219 

26500000 

3!>0 

le INbOi 3o0*360,370 

c 

220 

28600000 

3&0 


c 

221 

287COOOO 


00 TO 360 

c 

222 

208OOOCO 

370 

0UML=1 -0 

c 

223 

26900000 

3d0 

SJU.'il=AbS(SQkU i * O-OUM 1**2 ) ) 

c 

224 

29000000 


THE(HC=ATAN2i SUUMLfDUMl ) 

c 

225 

29100000 


ANAl'0=XK*Ay0 

c 

226 

29200000 


XKAY02=2,0*aKaY0 

C 

22 7 

29300000 


AKAt04=4 . 0*AKA¥0 

c 

22 8 

29400000 


XXY0S0=AKAY02**2 

c 

229 

29500000 


XN4= UM/ tXA**2* aKYOSO) ) ♦XKAYO 

c 

230 

29600000 


XN:>=lXN2/4AO**2tXKY0SQ ) )*XKAY? 

c 

231 

29700000 


XNt>=3. 14L3926+XKAYC2 

c 

232 

29800000 


XN7=3, 141 6926-XKAYQ2 

c 

233 

29900000 


XL2-XC**2 

c 

234 

30000000 


XN0=(aN3/12*O»(xC2^XN6 **2> ) )«XKAY0 

c 

235 

30100000 


XN9=lXiV3/ 42.0*1 AC2*XN7**2) >)*XKAYC 

c 

236 

30200000 


62V0i=l. 0/ <8ETAS0*AY0) 

c 

23 7 

303000U0 


CALL uAUSi ( i0,THELE*THEh1C*TBL10t 1I,KTBL101 IJ »OtO,tTKANU) ,HTRaN(1 

C 

23b 

3C4000CO 


U > 

c 

239 

30500000 


NPSI =3 

c 

240 

30600000 


jC2=3 

c 

241 

30700000 


JCUUNT=i2 

c 

242 

30000000 


CALL SSSSS ( JC2»OCOUNT,THESTN( 1> ,T»-fcWGH4 1>,LTKAM1J tKTKAN( in 

c 

243 

30900000 


THESTNi 2)=THfcMC 

c 

244 

31000000 


IF lYOJ j90*400,390 

c 

245 

3UOOOGO 

C 

FUR t06. 22*23 REF. 1 *♦* 

c 

24b 

31200000 

390 

X3*X (rluJ '-PS IMC 

c 

247 

31300U0O 


aY*6trASu*( Y0**2J 

c 

24b 

31400000 


THAT = 1 XK2*XS*AM2I /BETASO 

c 

249 

3I50000O 


COST=CDSlThAT) 

c 

250 

31600000 


6lNr=6IN4ThATi 

c 

251 

31700000 


GU TO 4iO 

c 

252 

31800000 

400 

BY-0 

c 

253 

31900000 

410 

00 330 I£=2,13 

c 

254 

32000000 

c 

m Uu LOOP 16T TIME FOR XSUE i 2ND TiMt FOR XSlMCi 3KU TO 

c 

255 

32100000 

c 

12TH TIME FOP 10-PT. GAUSS CUADKATURE *♦♦ 

c 

256 

32200000 


NLC.«^ = lfc'l 

c 

257 

32300000 


SiNTh=SlN4 THtST(^4NL0W) ) 

c 

250 

32400000 


PSI = .6*4PPP- 4 660) *COS< THE STN<NLOwJ J ) 

c 

259 

325UUOOO 


XO-X IM0)-PS1 

c 

260 

3^600000 


6 rHfc = THbSTN4NLUw) 

c 

261 

32700000 


IF lNL0«-2) 420.430*420 

c 

262 

32000000 

4^0 

X0TERM-A0/SURT4 X0**2-6Y) 

c 

263 

32900000 


XHEKt=XOTtRM*CGST 

c 

264 

33000000 


THEKE=xOTtKM*SlNT 

c 

265 

33100000 

c 

XbILE TERMS ARE COMPLETE *♦* 

c 

266 

33200000 


IF 4NL0W-1) 430i!>30*430 

c 

267 

33300000 

430 

CALL LCUMP 4PS1 . TANLE,PS!LE*AL41 ) ,STHE) 

c 

268 

33400000 


IF 4(NLUrt-2) b30*450.440 

c 

26 9 

33500000 

440 

CALL NEKNEL 

c 

270 

33b0C000 


IF lNLU«-3) 530,47Ct490 

c 

271 

337UOOOO 

4i>D 

00 **oO J I-I • b 

c 

272 

330OOOCO 

460 

ALMCiJD^ALi Jl) 

c 

273 

33900000 

c 

xsiMC Terms are complete *** 

c 

274 

24000000 


GU to !>30 

c 

275 

341000CG 

470 

UG 4bO LOw=I.o 

c 

276 

34200000 


FNSPKI LOW. j 1=^0 

c 

277 

34300000 


FNSPI4LUW. 3) =0 

c 

27b 

34400000 

4d0 

CCNT INOE 

c 

279 

34500000 

490 

CUMT ItsiOE 

c 

2B0 

34O00000 


00 d20 LUw=1.o 

c 

20 1 

34700000 


IF 4Y0) :>C0* DlO, D'^O 

c 

282 

34800000 

C 

ACCUMULATE iivlTEGRAL OF EQ 23. REF. 1 *** 

c 

203 

34900000 

t)00 

FNSPKlLOn*3)-FNSPR4LOW.3J + 4TH£WGH4NLQv«i)>*4 aakEAL*AL4 LuW ) ♦■ ALMC ( Low ) 

c 

204 

35000000 


1 *XhtkE ) *S1 nTh 

c 

205 

35100000 


FNSPI4 LUW.3 J =FNSPI ( LOW .3 ) + ITHEWGH INLOW ) ) ♦ 4 AK 1MA0*AH LOh ) - ALMC 4 lUw ) 

c 

286 

35200000 


l*rHERtJ *SlNTh 

c 

267 

35300000 


GO TO f>20 

c 

280 

354JU0UU 
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c 

ACCUMULATE INTEGRAL OF- EO 24» REE- 

1 


c 

269 

35500000 

310 

ENSPkiLUrt,3i =ENSPR ( LOW • 3 1 > ( THt WGH ( NLQW 

') )*XKhEAL*AL Iluw) ♦S INTh 

c 

290 

35600000 


ENSPIl LUiV, 3 1 = FN5P I <LOW,2 > ♦( T HE WGH ( NLOW ) ) ♦XKIMau^'AL (LUWJ *6lNTH 

c 

291 

35700000 

520 

CCNTINUE 



c 

292 

35600000 

d3 3 

CCNT INUt 



c 

29 3 

3590Q0C0 


IE (YOJ 3o0*54'^ 



c 

294 

SoOOOOOO 

C 

FOK eg. 22t kEF. 1 



c 

295 

36100000 

540 

XT^lXl Mg)-PilLbJ*^2 



c 

296 

3o 200000 


TRY=i>wKnxr-by ) 



c 

297 

36300000 


whk- iKy^c jst 



c 

296 

36400000 


WHI=TKY«SINT 



c 

299 

36500000 


UU d 50 KAN-lt6 



c 

300 

3o600000 


F5R{ lCANI-=-ULMC( tCAN) ) >^WHR 



c 

301 

36700000 

5b0 

FSU ICAN)=(ALMCl lCAN)i *WHI 



c 

302 

bobOOOOO 


UU rO 560 



c 

303 

36900000 

:>o3 

UQ 570 lUAN^i.o 



c 

304 

37000000 


FSKl lUANJ =0 



c 

305 

3710000C 

^70 

FSi 1 lUANJ =0 



c 

306 

37200000 

5d0 

UU 590 JC-l.o 



c 

307 

37300000 


FlP^iJCi =6B03«FNSPRUC *3UFSft( JC J 



c 

30b 

37400000 

390 

FlPl (JC)^bti0 5^FN$P!( JCt3)*FSU JO 



c 

309 

37500000 

C 

CrlLlKJrtlSt INTEGRAL F(ETA> OR F(Y) 

CQMPLtTt ^*** 

c 

310 

37600000 

c 

bEUlN MULTIPLICATION BY POWERS OF 

ETA AS IN tu. 16, ktF. 1 

c 

311 

37700000 


IF (NSYMi 6^Ct600.620 



c 

312 

37800000 

6Q0 

Nl N=1 



c 

313 

37900000 


00 olO Iu^l,l9t2 



c 

314 

38000000 


ANtNINi = ItTASTN* 10) 



c 

315 

38100000 

610 

N1N=NI N*i 



c 

316 

38200000 


GO TO 670 



c 

317 

38300000 

620 

M1N=L 



c 

31b 

38400000 


JU 660 1M=1»19*2 



c 

319 

38500000 


IF 1 tTASTNF li>) ) o5C»630f650 



c 

320 

38600000 

63 0 

IF (IM-i) 65Ctb‘*0t65O 



c 

321 

38700000 

640 

AN<M1N)=1.0 



c 

322 

38800000 


G Li re 66 0 



c 

323 

38900000 

U30 

ANIMIN)={lTaSTN( ID) )*«(IM-1) 



c 

324 

39000000 

6o 0 

MI N*MI 1 



c 

325 

3910000U 

670 

JKN = 0 



c 

326 

39200000 

C 

COMPUTE 40 TtKMS OF I INMSTAft, EO. 28 1 

KtF.l 

c 

327 

39300000 


DU 750 JLAN=1,6 



c 

328 

39400000 


IF FJCAN-2) 660*o90*7C0 



c 

329 

39500000 

c 

10 POWERS OF ETA FOP N^l AND 3t S 

POwtkS FQk *N-2 ANU 4f 

c 

330 

39600000 

c 

o FuK N=5f 4 FOR N-6. 



c 

331 

39700000 

660 

NTM^IO 



c 

332 

396U00UU 


GO TU 740 



c 

333 

3990000C 

69U 

NTM = B 



c 

334 

40000000 


GU TO 740 



c 

335 

40100000 

100 

IF UCAN'^) D60f69Cf710 



c 

336 

40200000 

no 

IF (JuAN-u) 720*730,730 



c 

337 

403U00O0 

720 

NTM=6 



c 

338 

*♦0400000 


GO TO 74U 



c 

339 

40500000 

730 

NTM-^4 



c 

340 

40o00000 

740 

DU 750 XCAN=1,NTM 



c 

341 

*t07000C0 


FFIRlKCAN. JCANl-F IPRUCAN)«AN( KC AN) 



c 

342 

40800000 


FFIl IKCAN, JCAN)-F'IPI iJCAN) • AN( KCAN) 



c 

343 

40900000 


JKIS= JKN^-1 



c 

344 

410000UO 


ANSR( JKNJ=hUhSTN( ID) *FFIR( KC AN , J CAN ) +AN SR( J nN J 

c 

345 

41100000 


ANSI i JKNJ =rtUhSTNI 10) «FFI I( KC AN , JC AN ) +ANS I I JKh) 

c 

346 

41200000 

750 

CUNT INUE 



c 

347 

41300000 

760 

continue 



c 

348 

41400000 


DU 770 JJ^lpJKN 



c 

349 

41500000 

770 

XM2STK(JJ)^CMPLXF ANSRUJ ) tANSUJ J)) 



c 

350 

41600000 


KMPT^NXYCNT 



c 

351 

41700000 

C 

uOMPuTE MATRIX IINM OF EQ- 28t RfcF* 

1 


c 

352 

41800000 


CALL MTaMPV IXM^STRU)) 



c 

353 

41900000 

760 

NXYCNT^NXYCNH'I 



c 

354 

42000000 

C 

LNU UF CCNTRUL-POINT DO LOOP 



c 

355 

421OO0C0 


NXY2=NXY*2 



c 

356 

42200000 


RETURN 



c 

357 

42300000 


END 



c 

358- 

42400000 
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KERNEL .- The function of subprogram KERNEL is to calculate K of the super- 
sonic kernel function of equations (3) to (6) of reference 1. The evaluation of the integral 
in K is based on the approximation of the integrand as in equation (A5) of the appendix. 
The flow chart for KERNEL follows. 



SuaftUUTlNt KtftUfcL 

c 

C COHfuTt Tht K3 aR of EQS. ^ AND <♦, R£f. I 

C THt algorithm iS OeSIGNED TO MlIxIMI/E CALLLlAT iUNi UF iiNfcS, 

c CjSimcs, anj exponentials, the integrant lf The tau integral 

L IS APPROaImATED AS IN EO. 21 Of NASA TR K-AO. CHUKUkIsE 

C INVARIANT gOANTITIES ARE CALCULATED IN SUtjROuTINE PAkTI AND 

^ PASSbU VIA CuMMON- 

C 

JlMbN:>IUiM hT(T^6) 

UlMtNSIUN X i , t ( ^8) f IDENT I 8) 
iJl-'ltNb lUiM UHX(^d«6) fH[ 48f 6) 

COMPLEX XMATAdl4t)#16) 

CaMlH0N/TRA21o2/AM, XK, T ANL E , T AN LT E • S ,NM0OE VH , KPRuo 

CUKMUi^i NXt tL>fcL,X,Y,DHX»H»^T 

CCMMuN AA,At5,XC,xN4,XN5tXN6,XN7, XN8tXN9, BtTAi> w,b2V 01, xO,yO,XXOB2, 
lXKMd2,XKAY02,x0Ti_RM,SKAY0,XKY0Sg,XKR£AL,XXjMAO,NkKUh tX«AT 2 ,KMPT 
XKArU4-2.0*AKAY02 
IF UO) 20,10,20 

10 UUMbA=^ 12. 0*XkJ «a0 

XKKfcAL=^-CUSiUUM8AJ 
xMMAG=iilNtUUMbA) 

Gu To 100 

20 IF (XKJ *fO,3G,40 

3U AMtKM-A0«*2-dtTAS0*( ) 

XXkEAL=-AO/bwRT( AKTERM) 


D 2 
U 3 
U 4 
0 b 
0 6 
D 7 
Li 8 
D 9 
0 10 
D 11 
D 12 
0 13 

0 14 

U 15 
0 16 
U 17 
0 18 
0 19 

D 20 
0 21 
D 22 
0 23 

0 24 

U 25 
U 26 
0 27 


42600000 

42700000 

42600000 

42VOOOCO 

43000000 

43100000 

45200000 

433000UO 

43400000 

43500000 

43600000 

43700000 

43600000 

439000CO 

4400000U 

4410000U 

44200000 

44300000 

44400000 

44500000 

44600000 

4470000C 

44600000 

*♦4900000 

45000000 

'♦5100000 


18 




o o cr o o n 


XK 1MAG=0 

0 

2b 

452000CO 

TU lOQ 

0 

29 

45300000 

SwKTl=SgKl tXO#«Z-BETASQ« (Y0**2> J 

U 

30 

45400000 

AUPL=AO+ aH«Su^T1 

0 

31 

45500000 

XJMZ-aO'XM*SwKT1 

0 

32 

45600000 

AXQMI=AbS(XLMl J 

0 

33 

4570C000 

TUP=b2YUl« auPL 

D 

34 

458U00U0 

ABTLu=b2Y D1 <=AaOMI 

0 

35 

45900000 

aN 10=aKu62*aO 

0 

36 

46000000 

XM l=AKM62*i> UkT I 

0 

37 

46100000 

PTUP=( i^lDV2oJ 

i) 

38 

46200000 

PTLU = (3.l4l:jS2o) ♦AbTLO 

D 

39 

4o300000 

AT U- aA*TUP 

0 

4C 

46400000 

ATl= AA*AciTtU 

u 

41 

46500000 

oTU*xa«TUP 

0 

42 

46600000 

oTL=Xb<‘AoTLu 

D 

43 

46700000 

C.TU-XC^'TUP 

D 

44 

46600000 

LTL=xC*AbTLO 

D 

4b 

46900000 

EAPL=tXPl ATUJ 

D 

46 

4/OUOOOU 

cAP2=^EXPl ATL J 

0 

47 

47100000 

E AP3 = t XP t bTU) 

i) 

48 

47200000 

tXP^-tXP ( bTL i 

0 

49 

47300000 

EXP5=EXP ( CTU) 

0 

SO 

47-^OCOOO 

t APo=t XP 1 CTL < 

0 

bl 

4750000U 

CSl^tuSl ANIO ) 

U 

52 

47600000 

CS2=tUii XNl 1 J 

0 

53 

^7/00000 

SN1 = SI Ai^l 0 i 

0 

54 

47600000 

SN2=Silsi( XNll) 

U 

5b 

4790C000 

XN12=LSI»CS2 

D 

56 

4b0UUU00 

XNl3=SNi♦S^i2 

0 

57 

46100000 

CS3=»XNU'AN13 

D 

5B 

46200000 

CS4=X/\il2 +ANL3 

D 

59 

4630C00C 

XN1A=SN1<‘CS2 

U 

60 

^6400000 

AN1S-SN2*LS1 

0 

61 

4b5UUOOO 

Sivl3-AMt*-AN liJ 

D 

62 

4bo00000 

SN<»- XIU^-XNL^ 

0 

63 

-♦6700000 

CSi)^CUSlPTUP) 

i) 

64 

4bb00000 

C:i6=C0SlPTLJ J 

[) 

o5 

48900000 

SN5=S1N1 PTJP J 

0 

66 

49000000 

SNb=SlN< PTLU) 

0 

o7 

49100000 

IF ZXUMlj bOfOOfwO 

0 

6b 

49200000 

SN4=-SN4 

D 

69 

49300000 

$No--bNb 

D 

7C 

49400000 


D 

71 

49500000 

REAL PART 

D 

72 

49600000 


0 

73 

49?00000 

kU2* l^i‘43’-bN4 i /2 *0 

D 

74 

49800000 

Rg3= UXA*CS3tAkAY'D2*>SN3J *E XP 1- ( XA*CS4*-XKA Y02<‘ *tXP2 J *xN4 

D 

75 

49900000 

Rg4= U Xb«c:>3tXNA y 0 2^SN3) ♦EXPB- (XB*CS4*-XKAVu2*bMHi#txP*>Z *Xn5 

J 

76 

50000000 

xN2a-CSb*CS3 

0 

7 7 

501000CO 

aN^I-3N3^SN3 

U 

7b 

50200000 

C S 7 = Xi'^ 2 0 - X (N*2 T 

0 

79 

5030O0Q0 

CS9- XN20f XN2 1 

D 

80 

50400000 

XN22=CSo«CS4 

0 

bl 

50500000 

XN23=i Nt)«SN4 

D 

82 

50600000 

CSb=XN22-XN23 

D 

83 

50700000 

C SI 3-AN22 ♦ XN23 

D 

84 

508U000U 

XN2^=3N:j*CS3 

0 

85 

b090000C 

XN23=Li5*SN3 

D 

86 

51000000 

SN7 = XN244'AN2j 

U 

87 

5110000C 

SN9*xN24-XN2:i 

0 

86 

51200000 

A^2^3 = SNb♦C S4 

D 

89 

5130000U 

XN27=CSb*bi^^ 

D 

90 

514000CC 

SNd = XNZ6*-AN27 

D 

91 

51500000 

SN10=aN2i>-XN2 7 

D 

92 

51600000 

U t XC=*SN /-XiVb*CS7) ♦EXPS^ l-XC ♦SNb+XNb*LSbJ XP6 i *ANb-»- 1 1AC*SN9-X 

0 

93 

5170000C 

Lr 4 7*CS9)*tAPS«'l-Xi:*SN10fXN7 ♦CSIO) *ExP6)4XNYi 

D 

94 

51600000 


0 

95 

51900000 

iMAblNAkY PmKT 

D 

96 

52000000 


D 

97 

5210C00C 

If iXJMn b0»70t70 

0 

9b 

52200000 
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non 


70 


60 


90 


100 


Aiw2=-i-CS3^C:>4)/2.0 

XIU3=( (-XA^SN34^ANAy02^CS3)*eXP l*-< XA*SN4-XAAYQ^*C St) * LXP2 J ♦ AN^ 

XIU4 = U -X6*SN3+AKAY02*CS3J ♦EXP3+(XB*SN4-XKAY02<«CS4i^tXPt)^AN3 
X103 = U ( aC*CS7<'AN6 4S^7)«EXP^-( XC#CStH-XN6«ii^<ttl*tAP6i^AM6+( ( XC<‘Ci»9 + X 
lN7»SNy) >»“£ aP^-( XC*CS10^XN7*SN10)^EXP6l^ I-XNV) i 
GJ TO 90 

Xlu2c>(-CS3-CS4^2.0)/2.0 

a1w 3=< i'XA<«SN3 + AKAY02*CS3l XP { X A«SN4-XXAY 0«i»C^4 )^EXP 2-XKAY 04) *X 
IN^ 

XlU4=t < -Ab»SN3 + XK AY02*CS3) *EXP3- ( Xb«SN4-XXAV u2«LS4 ) «LxPt*XKAY 04) *X 
1N5 

XigS = ( ( i XC^C!>7^XN6*$N7 > *EXP 5+ [ XC Sd'HXNO^SNb ) ♦tXPo-^ • AC ) ♦XNb^i iX 
IL*CS9+XN7#Sr49)*tAP5* (XC*CS10-»-XN7«SN10) ♦Expo-a • 0 *aL J “ANy ) J 
Kgi5:Ru2tKw3'»’Kgt-i-K0 5 
XIwl=XIw2^Al w3^AlQ4fXI 05 

integral in a bak finished 

XN31~A0TtKM*L:i2 
XN32 = AN31*0:>1-Algl 
XN33=-XN31<‘bNH-Rwl 
GS=2,U*XX«X0 
Cill=gOS(L>S) 

SNil=SlN(GSJ 

AKKLAL = -lXN3i+C^Il ♦XN3 3*SM1 ) 

XKIMAij = XtJ32»SNli“XN33*CSll 

RhTUHN 

END 


0 

99 

52300000 

0 

100 

>2400000 

D 

XOl 

52500000 

0 

102 

52o00000 

D 

103 

52 /OOOOO 

0 

104 

52600000 

D 

105 

529COOOO 

D 

106 

53000000 

D 

107 

531000C0 

U 

108 

53200000 

0 

109 

53300000 

0 

110 

5>4C00C0 

0 

111 

535UOOCO 

D 

112 

53600000 

i) 

113 

53/OOOOC 

D 

114 

53800000 

0 

115 

53900000 

D 

116 

54000000 

U 

117 

54100000 

0 

iia 

5420000U 

0 

119 

5430000C 

D 

120 

54400000 

D 

121 

54500000 

0 

122 

54O0000C 

D 

123 

54700UO0 

0 

12<» 

54800000 

D 

125 

5490000U 

0 

126- 

5500000C 


MTXMPY .- The function of subprogram MTXMPY is to compute the elements of the 
matrix III of table I of reference 1 and carry out the indicated matrix multiplication of 
equation (28) of reference 1. The flow chart for subprogram MTXMPY follows. 


Cmtxhpy"") 



XKAT3 = COLUMNS 5 TO 15 
OF III MATRIX 

EQ. (28) AND TABLE I, REF. I 


XMAT5 = COLUMNS I TO 5 




WRITE ON TAPE9 
RECORD It COLUMNS I TO 4 
RECORD 2, COLUMNS 5 TO 16, 
OF MATRIX HI 


( 100 V- 

X 


ONE ROW OF IINM MATRIX, 
EQ. (28), REF. 1 


(return ^ 


20 






SUbRUUTIuE HTXMi>¥ (XM2STP) 


E 

2 

55200000 

CDMPiiTt ThE iiT MATRIX OF EUS. 28 TO 3i»ANU TAbLt i OF KEF, 

i- 

fc 

3 

5530000C 

ANU HOLTIKLX as in EC, 26 TO GET IINH ^aTKIa Of EuS, 18r2bf 

32 

E 

4 

55400000 

KEAL «<;.^b*KAC«K40»K4 


E 

5 

53500000 

REAL KE,KF,k2»*I^2C*K20,K2 


E 

6 

55600000 

DlMtNSiUW wT I'ib) 


E 

7 

55/00000 

Dir^tN:^lDN bl9l *F125>tKEI SJ «KF( 


£ 

8 

55800000 

OlHtMSlUn X(<»o>»V(4BI«TOEKTIB> 


E 

9 

55900000 

UiNENSiQN UHX44b,6)»H(48.6) 


E 

10 

56000000 

COMPLEX XMAT314btl2l»XHAT2I48t I4>) ,XM2STR(<i6J 


E 

11 

56100000 

complex XMATh(%0»4> 


£ 

12 

56200000 

CaMHaii/TKA2ib2/AM»XK,TANLF,TAf^LTE.SfNMQDt ,NSVM,XPKGG 


E 

13 

5O300000 

CUMf^QM f4AT*UEL*A.Y*OHX.M,bT 


E 

14 

56400000 

COhMOU XA,Xb ,XCf XN4vXM5»XI^6»XN7» XN8tXN9«BE7A6U»82Y01rX0« V0|AK0d2 

f 

£ 

15 

56500000 

IXKMd2«XKAY 02 .XOTER M.SKAYO » XKYOSU » XKREAL * XxlhAG t NKKOx t XMAT2 ,KMHT 


e 

16 

56600000 

x2=2,0*xx 


E 

17 

56700000 

JUMP TQ lOCiO FUK OTHER THAN FIRST CONTROL PUINT 


t 

18 

:>6bOOOOO 

IF 4KrtPT-U lOpiOflOO 


£ 

19 

56900000 

COMPOTE ELEMENTS OF 11 I MATRIX, E0S,29 TU 3i ANU TAdLfc 1» RfcF. 

1 

E 

20 

57000000 

811 = 0, !> 


E 

21 

57100000 

tti2=0.2113248b 


£ 

22 

57200000 

822=0. 788^7314 


£ 

23 

57300000 

613=0.11270180 


E 

24 

57400000 

82 3=0. D 


£ 

25 

57500000 

833=0. U8/29634 


E 

26 

57600000 

Cl 1 = 0. /5^11.Q>SI»ITANL6^TANLTE i 


E 

27 

37700000 

C12 = CU* 1 0.453012/C/S) TANLE-TANLTE) 


E 

2B 

578000CC 

C22=2.0'^til-C12 


e 

29 

57900000 

C13=Cll4^*0.3o09^7Sn/S) »I TANLE-TANLTEJ 


E 

30 

58000000 

C23=Cll 


£ 

31 

58100000 

C33=2.0*L11-l13 


£ 

32 

58200000 

CNST=4-l.0/l 5.0*3*^2)) 


h 

33 

58300000 

iUi = CNST^Cii 


E 

34 

58400000 

012=CNST«C12 


t 

35 

:>850000C 

022=LNST«C22 


H 

36 

58600000 

013=CNST*C13 


£ 

37 

58700000 

023=CNST»C23 


b 

38 

SboOOOOO 

033=CNST«C53 


£ 

39 

58900000 

E4 IMo.O 


E 

40 

5900 0000 

t < 2l=<-6.Q^822-6.Q*B!2) 


E 

41 

59100000 

E4 3) =4-o.0^C22-o."^*C12) 


E 

42 

59200000 

E 14) =4-8.0*022-0. '>*012) 


E 

43 

59300000 

1 4 !>) =1-0.0*812 )*1 -622) 


t 

44 

59400000 

t{8) = -O.0*Bl^*l-C?2)-6.''*C 12 *4-822) 


E 

45 

5950 OOCC 

1 4 7)=-8.U*bl2*4-o2 2 )-6.0*C12*(-C2 2 )-6. 0*012*4-822) 


E 

46 

59600000 

E46)=-8.0*C12*4-022)-6.?*D12*4-C221 


E 

47 

59700000 

t 1 9)=-6. 0*01^*1-022) 


E 

46 

59800000 

FllJ TO F428) HERE APE G(l) TO G425) IN REF. 1 *** 


t 

49 

59900000 

f 4 U=-20.0 


E 

50 

60000000 

FI 2) =20. 0*825*20. 0*613 


E 

51 

60100000 

F 4 3)=20.0*C2:s*20.'^*Cl3 


E 

52 

6020000C 

F 1 4) =20. 0*023*20.0*01? 


E 

53 

60300000 

F 4 5)=20- 0*bl3*4-o23 ) 


E 

54 

O0400G00 

F 4o)=20.0*bl3*4-C2? )*2 0.C*C1 3*1-B2 3) 


E 

55 

60500000 

F 1 7)=20.0^813*l-O23)*20.0*C13*4-C23>*20-0*Oi3*(-825i 


L 

56 

60600000 

F18)=2 0.0*C13*l-J23)*20.'"*l313* (-C23) 


e 

57 

607000CO 

F 19) =20. 0*013*4-023) 


E 

58 

bOdOOOCO 

F 110)=F4 I > 


E 

59 

60900000 

F4 ll) = F42)-b33*F41) 


E 

60 

610000C0 

F4 i2)=F43>-C33*Fll) 


E 

6i 

ollOOOOU 

F (lj) = Fl4)-033*F4 1 ) 


E 

62 

6120000C 

F 4 l4)=F4 >)“8 :>j*F12) 


E 

63 

61300000 

F 4 l3) = F4o)-633*F 4 3 )-C3 3*F4 2) 


t 

64 

61400000 

Fl 18) = F4 7>-833*F 4 4)-C??*Fi3)-033*F4 2l 


E 

65 

61500000 

F 4 1 n=F 4 d i-C33*F 4 4 1-033* f 4 31 


E 

66 

616UQ0Q0 

F 4 Id ) = F4 9)-0 33»F 4 4 1 


E 

67 

61700000 

F< 19) = -833*F 18) 


E 

68 

618000C0 

F(2aJ=-833*F(b)-C3 3*F4 5) 


E 

69 

61900000 

F 421> = -633*F4 7)-C??*F{ 6)-D??*F( 3) 


£ 

70 

62000000 

F 4 22)=-b33*F \ d)-C3 3*FI 7)-033*F 4 6) 


E 

71 

62100000 

F 4 23)=-83 3*F4 9)-C33*F4 8)-D33*F(7) 


E 

72 

62200000 
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Fl24> = -C33<‘Fty>-u33*FJ 8) 

UU 20 J=lf^6 
□0 20 L=lfl2 
20 XMAT31 JfL> = U**0*) 

K2B=2*0*AK*l^bUi 

K2C=2^0*XK*(-C11> 

K20-2.0*XK*l-Dlli 
6STAk=2*0^bl 1 
CSTAK=2*0«Cii 
OSTAR=2.0*D11 
K^tt=— ♦.0*AK* l-bLi ) 
kaC=-a* j«xK^ i-C.il i 
k^O--4pO*XK^ i-OIl * 

L COMPaTE C0L«JMI^:> 5 TO 3 CF III MATRIX* 

JO 30 

XHAT3t Jt JJ =CNPLXC8STAR ,0.1 
L = J>i 

XMAT3I L * Ji=LMPLX< CSTAR»0. J 
L = L+1 

AMAT3Ut Ji-LMPLXiOSTAR *0.) 

L = L^8 

XMAT31L* J* =CMPtXi-2 .'^,0.1 
L = L*-b 

XMAT3( L ,JI=CMPLa4-2.0*K46> 

L = L + 1 

XHAT3< L ,0*=LKPLa(0. ,KAC> 

L = L>1 

XMATSiL, JJ^LMPLAi 0. *KAO> 

L = L*-6 

30 XMAT3lL»b>^CMPLXi 0-tKA ) 

JO ^0 J-l,y 

^0 Kti JI=2.0«XK^E<ji 

E T AM==t 4 1 >«TANLfc**24E C6 I 
C CJMPUlt COLUHfViS 9 TO 12 OF III MATRIX 

liO 

XMAT31L, J>=CMPLXi Ei 5»,0. J 
L = L4l 

XHAT3i L# Ji^CHPLXlETAN, I 
L=L4l 

XMAT34L, j4-CkPLXi 6(7), 0. I 
L=L4l 

XMATii L, J)=CMPLX4E I a) ,0. ) 

L=L4l 

XMAT34 L# J) »CFPLAiE4 9) .0. ) 

L=L46 

XMAiai L, Ji^CWPLXiE (?),0. ) 

L=L4l 

XMAT34 L, J4 = lMPLa 4 El 3i,0. ) 

L-L + l 

XMAT34L * J) =CFPLa4E4AJ,0. ) 
i- = L4b 

XHAT31L, J)=tMPLX4E(2) ,KE(^) ) 

L = L4l 

XMAT34L, J)^CMPLA( E43I ,KE (61) 

L = L+l 

XMATSa, J)^CMPLXIE( A) tKE 17) I 
L^L4i 

XMAT34L, J)=cMPLa4^.,KE 4 9) ) 

L=^L41 

XMAr3(Lf J)=CHPLXl0.pKEI9)> 

C=L46 

XMAT34 Lp J i-CHPLAi 1 3 .?♦ E 1 1 ) ) , KE 12 > ) 

L = L4l 

XMATjiLp J)-C MPLA4 0.pKE43l) 

L^L+l 

XMAT34LP J)=CMPLXIO.pKElA)) 

L=L46 

bO XMAT34L, J) = CMPLa4^.,KE ID) 

FTA^41^*F 4 15>>4Plia)»TANtE“4»2 




t 

73 

62300000 



E 

7A 

6240000J 



E 

7b 

62bOOOOO 



E 

76 

6^600000 



e 

77 

62700000 



E 

78 

62800000 



E 

79 

62900000 



E 

80 

63000000 



E 

81 

63100000 



E 

82 

63200000 



£ 

83 

63300000 



E 

84 

634000CO 



E 

8b 

63300000 
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SSSSS .- The function of subprogram SSSSS is to store in sequence the span stations 
7] and the associated spanwise integrating factors for subregions I, n, HI, and IV of fig- 
ure 3 of reference 1. No flow chart is needed. 

GAUSS .- The function of subprogram GAUSS is to calculate the span stations r} 
and select spanwise integrating factors for a 6-, 10-, or 16-point Gaussian quadrature in 
regions I, II, and IV of figure 3 of reference 1. No flow chart is needed. 

LCOMP .- The function of subprogram LCOMP is to calculate the values of ^^*(^, 77 ) 
(n = 0,1, 2, 3,4, 5) of equation (13c) of reference 1. The flow chart for subprogram LCOMP 
follows. 
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30 WbHFLTl i )=bMA^Kli )/VALUE 

RETURN 
END 

SDBRUUTiNc LCuMP 1 DP$1 , T AN ,DPS I LE » XLNOt OTht ) 

C 

C COMPUTE Iht uNSTAR (N=C TQ S> CF EU* liCt REF. 1 

C 

DlMbNSlUN XLNUloJ 
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GO TO 80 
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RETURN 
END 
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PRT2.- The function of subprogram PRT2 is to solve equation (32) of reference 1 
for the weighting factors a^. A heading is written. If HMLTPLY(l) and HMLTPLY(2) 
are both nonzero, the real part of the downwash input (-3h^/ax) is multiplied by 
HMLTPLY(l), and the imaginary part (-2hj) by HMLTPLY(2). The imaginary part is 
also multiplied by k. When there are more than 16 control points, the simultaneous 
equation set is solved by a least-squares subroutine GLSP. The resulting a[^^ matrix 
is printed, written on TAPE9 for communication to D2182, and can be optionally punched. 
The flow chart for PRT2 follows. 
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GLSP .- The function of subprogram GLSP is to solve for the weighting factors 
a^^^ of equation (32) of reference 1 with the aid of subprogram SIMEQ. A weighted 
least- square-error solution is made if the number of equations M exceeds the number 
(N = 16) of unknowns. The matrix of weights is supplied by the user. For WT(1) = -1,0, 
all weights are made 1.0. After a reduction to a normal form, solution of the simulta- 
neous equations is done by SIMEQ. If the number of equations equals the number of 
unknowns, SIMEQ is called directly. If the number of equations is less than the number 
of unknowns, solution is not made. Instead a RETURN is made with lER = 1. The list- 
ings for subprograms GLSP and SIMEQ follow. 



SUBkOUTlNt GLSP ( A , M, N , 6 , I P , X, U i SUM , rtT t 1 EK , I 0, t.A , H AXM , MA*N J 


J 

2 

93/00000 

C 

DuCuMLNT UATc 08-^1-68 SJbROLTINE RtivI3tu 07-09-71 

J 

i 

96800000 

c 

This KUUTINE iS ADAPTED FCR USE AT LANULEY i-kuM SMlhSUNlAN 


J 

4 

93900000 

c 

ASTRUPHYSICAL OBSERVATORY SUBROUTINE GL6P* 


J 

3 

94000000 

c 

A4MAXh,M«xN) = GIVEN RECTANGULAR MATRIX 


J 

6 

94100000 

c 

K = NJ. kOwS of a 


J 

7 

942COOOO 

c 

N = Nu, Cuts OF A 


J 

8 

993000CO 

c 

bJMAxM,iP> - MATRIX CF RIGHT SIDE 


J 

9 

94400000 

c 

IP = NU* COLS OF B 


J 

10 

943000C0 

c 

X|MaaN*IPJ = SCLUTION MATRIX 


J 

11 

94600000 

c 

U4MAXM,IP) = RESIDUAL MATRIX 


J 

12 

94 700000 

c 

SUMUP) == SUM OF WEIGHTED SQUARES OF RtSiiiUALS 


J 

13 

94SOOOOO 

c 

WT4MAXM4 = WEIGHTS, IF WTm='l.Q* ALL wElbhTS AKL SET = 

1. 0 

J 

14 

94900000 

c 

lER = ERROR RETURN, = 0, NORMAL* = Ir M LT N. - DE T 

= 0. 

J 

13 

93000QC0 

c 

TbIMAXN) TEMP, STORAGE IN SIMEQ 


J 

16 

93100000 

c 

WA4MAXN, MAxn) = MATRIX OF NORMAL EQUA7 luNS 


J 

1 7 

93200000 

c 

MAXM = MAXIMUM ROWS OF A 


J 

18 

9i>300000 

c 

MAXN = MAaIMUM cols of A 


J 

19 

93^00000 

c 

SlMhW = ROUTINE TO SOLVE NORMAL EQUATIONS 


J 

2C 

93300000 

c 

UlMENSIuN A(M^AXM,MAXN) , B(MAXM,IP>, XtMAXritlPJt uIMAXM,IP), 


J 

21 

96600000 



J 

22 

95700000 


1 SUM{ IP J , wTIMAXMJ , TBIMAXN > , WA 4 MAXNtMAXNI 


J 

23 

958000CO 

c 



J 

24 

959UOO00 


DtT=0*0 


J 

2b 

96000000 


I£K=0 


J 

26 

96100000 

c 

SET ALL wEIGTS = 1-0 FOR NO WEIGHTING IwTlU - -i,0J 


J 

27 

96200000 


IF I WTI ll.Gt-0*) GO TO 20 


J 

26 

96300000 

10 

DU 10 I ^^1 * M 


J 

29 

964UOOU0 

wTU J^l-Q 


J 

30 

96500000 

2J 

IF 4M-NI 30.40,70 


J 

31 

96600000 

c 

IF M LESS THAN N, NO SOLUTION, I Ek = 1 


J 

32 

96700000 

30 

lER^l 


J 

33 

9o800000 


return 


J 

34 

96900000 


28 



c 

it- = Nf solution with normal EQUATIOnSi ax ^ tt 

J 

3^ 

97000000 

40 

uu 60 U1,M 

J 

36 

97100000 


00 ^0 J=1,N 

J 

37 

97200000 

t>0 

WAU ,J j=A( i. JJ 

J 

30 

97300000 


00 ijQ K==l t IP 

J 

39 

97400000 

60 

X« l»KJ=bt I,K) 

J 

40 

97500000 


bO TO 100 

J 

41 

97600000 

C 

Ik- M uT Nt rtA ^ A<T) ♦ A ♦ WT , X = A,(T) ^ 0 ♦ WT 

J 

42 

97700000 

c 


J 

43 

97800000 

70 

00 dO K=i *N 

J 

44 

9790 UOOO 


UO dO 

J 

43 

96000000 


rtAi J.KJ=0.0 

J 

46 

98100000 


UG 60 1=1, M 

J 

47 

982000C0 

60 

WA( J,Ki=:rtAlJ,XJtAn iK) ♦A(I,J J+WTU) 

J 

4b 

98300000 

C 

WEIGHT B 

J 

49 

98400000 


JU K= 1, IP 

J 

30 

98500000 


GO 'JO J = 1,N 

J 

51 

90600000 


X U, Ki 30.0 

J 

52 

96700000 


00 90 1=1 ,M 

J 

53 

9860U00U 

90 

XU,K)=Xl J#Ki tAi I* J)*Bn ,K)^WT(I ) 

J 

54 

98900000 

C 

SOLVE twUATIuNS BY SIMEQ 

J 

55 

99000000 

100 

CALl SIMty ( ifc A f N f X , I P f OE T f T B f M AX N f I S ) 

J 

56 

99100000 


(DET-Nb.O-i GO TO 110 

J 

57 

99200000 


iek^2' 

J 

5b 

99300000 


bO TO l<tO 

J 

59 

99400000 

C 

SULVt KfcSlOUALS AND SUM OF SgUARES OF RtSiuUES 

J 

60 

99500000 

110 

OU 130 K«i # 1 P 

J 

6l 

99600000 


SUM(XJ=Q.0 

J 

62 

99700000 


00 130 1=1, M 

J 

63 

99800000 


UiUKi =-611, K) 

J 

64 

99900000 


UO 120 J^1,N 

J 

65 

lOOOOOOOO 

120 

Ui I ,K> =un ,KJ At I, J)^X1 J,K) 1 

J 

66 

lOOlOOOOO 

130 

SUMlK)=SUMlK>+40t ItK) l♦*2♦WT^ ) 

J 

67 

100200000 

140 

RETURN 

J 

68 

100300000 


ENU 

J 

69- 

100400000 



SUBKOuTlNE SIMEU (A*N,0tMf 

OETeRM,lPlVOT,NMAX»lSLALEJ 

K 

2 

100600000 

c 

♦ UUCuMLNT UATt Ofl-''l-68 

SUfaROLTINE RtViStU 08-01-66 

K 

3 

100700000 

c 

iULUTIUN OF SIMULTANEOUS LINEAR EQUATIONS 

K 

4 

LOO 0 OOOOO 

c 



K 

5 

1U09000UO 


UlMtNSIUN IPIVUUNJ,AiNMAX 

,N),0iNMAX,M) 

K 

6 

lOlCOOCOO 


t JUI VALfcWCt 1 ikUWt JROW ) f ( I COLUM, JCCLUM) t { AMAXt T , SWAP J 

K 

7 

101100000 

c 



K 

0 

101200000 

c 

initialization 


K 

9 

101300000 

c 



K 

10 

1U140Q000 

10 

ISCALt=0 


K 

11 

101500000 


Ri=iO. 0*»100 


K 

12 

lOIbUOGOO 


R2=1.0/R1 


K 

13 

LOl TOOOOO 


UETcRM = l*Ci 


K 

14 

101600000 


UQ 20 J = 1 1 N 


K 

15 

101900000 

20 

IPI VGT t J) -0 


K 

l6 

1020000CO 


UU 3b0 I“lfN 


K 

17 

1U210D000 

c 



K 

IB 

1O220OC00 

c 

6LARCH FOR PIVUT ELEMENT 


K 

19 

102300000 

c 



K 

20 

102400000 


AMAX=0-0 


K 

21 

102500000 


UO 7C J=lfN 


K 

22 

102600000 


IF liPIVOnJi-1) 30,7:»,3C 


K 

23 

. 102 7000CO 

30 

UU bO R=1,N 


K 

24 

102800000 


IF ( IP IVUT (K J-l> 40,60,39'" 


K 

25 

102900000 

40 

IF lABSt AMAXJ-Abcl A(J,K> ) ) 

SOi 60,60 

K 

26 

103000000 


2 RQw = J 


K 

27 

103100000 


1CULUM=K 


K 

28 

i032COOCO 


AMAX = A 1 J , X i 


K 

29 

iU3300U00 

60 

CuNTINUE 


K 

30 

IOBhOOOCO 

70 

CCNTINUt 


K 

31 

103300000 


IF tAMAXJ 90,60,9^ 


K 

32 

103600000 

80 

UhT£RM=0.0 


K 

33 

103700000 


I SCALL^O 


K 

34 

IO30OOUOO 


UU TO BVO 


X 

35 

103900000 

90 

1 PIVOT UCGLuMi = I PI VOT( 1 COLUM 

K 

36 

10400C000 


29 



c 

L iiSiTliKCHAiVbt TO PUT PIVOT ELEMENT UN uImuun^^L 

C 

If t iKOw-lbULUi^i 100rl40,100 
100 LJbTt PM--L>t TLKM 
UU 110 L=1,N 
!> =A ( 1 RUw • L 1 

AUKOh*L J =Al lCuLuy,U 
110 A( Il(jlUN,L ) = S«AP 
If (MJ 

120 OU 130 L=ltM 

S«AP = ti URUvJt U 
6 i lKO«,Li = b{ ICuLUMtL J 
130 8( ICULUMtL J = SrtAP 

l<rO PI ICULUM, ICOLUM) 

If IPlVUn lt>0»b0.1!>^ 

c 

c SCALE THh UfcTfckMlNANT 

C 

1:>U PIvuTi=Pli/UT 

If ( AbSl JET tKMi“Kl ) 180,lbCtl60 
UjQ JtTEKM=UcTERM/Kl 
lSCALt=lSCALLfl 

If ( AbSUibT tKN J-Kl J 2l'^#lTC»17^ 

170 UtTbKM=Ut TEKM/ kI 
ISCALt^iSCALb^-l 
OU TO 210 

ItiO If ( AoSt UbTtkM I-R2 ) l90tl9C,210 
IVO OETtRM^UtT 

1 SCALfc-i SC ALb-l 

If I AdSluETtkMi-KZ) 200»2^C,210 
200 UfcTERM=UcTEKh*Rl 
I SCalL^I SLAL t-1 

210 If UbSlPlVOin-KU 24C, 220^220 

220 PI VUTl =P WOf I/Rl 
lSLALt*I$CALL+l 

If lAbS(f I VuTi)-Kl} 270*230,230 
230 PlVUTi*Pl VOTI/Kl 
ISCALb = lSLALL^‘l 
bO TU 270 

Z<*0 If lAtiSI PIVuTn-k2) 250, 25^,270 
25U PIVUTi-PUUTI’^Rl 
ISCALE*lSCALb-l 

If lAbSi P I VDT 1 )-k2) 260t260,270 
260 PI VUTI=PJ VOT l*Kl 
I SCALt=*ISCALb-l 
27.0 OE TriKM=uETtki^»Pl VOTI 
C 

C Ui^iQE PIVOT KUd BY PIVOT ELEMENT 

C 

00 290 L=1,N 

IF (lfJvdT(Li-n 280,290,390 
280 A( ICUlUM,L)=A( iCDLUM,L )/ pivot 
290 CLNTlNUt 

IF iM) 320.320,30'' 

300 DU 310 L-1,M 

310 bllCGLUMtU=bI ICGLUMtL J/ PIVOT 
C 

C KtUUCb NGh-PIVOT RCWS 

C 

320 OU 3o0 Li^l.N 

If (Ll-ICULUMI 330,380,330 
330 T=Ai LI , 1 LULUM> 

UU 350 L=1 ,N 

If UPIvGTIlJ-U 3AO,350,39'> 

3^U AIL1,U=A(L1,LJ-A(ICGLUM,L)«T 
3^0 CLNTINUE 

If IMJ 3b0,3b0,3u0 
360 DO 370 L=l,M 

370 b Ui ,L I =bl L I ♦ L)'b< ICnLUM,L)<‘T 
3b0 CCNTlNUt 

390 RETURN 

END 


K 

37 

104100000 

K 

3b 

10^200000 

K 

39 

i 0*1-30 00 

K 


1 UA^OOUOO 

K 

<♦1 

1045C00UQ 

K 

A2 

IOhoOOOOO 

K 

A3 

104700000 

K 

AA 

1 UAbUOOug 

k 

A5 

104900000 

K 

46 

1 Q:>OOOOGO 

K 

4 7 

105100000 

K 

4d 

10:^200000 

K 

49 

1U53U00U0 

K 

50 

10:>40000C 

K 

51 

105500000 

K 

52 

105600COO 

K 

53 

105700000 

K 

54 

1 05800000 

K 

55 

105900003 

k 

56 

106000000 

k 

bJ 

106100000 

k 

58 

106200000 

K 

59 

106300000 

K 

60 

ICc* **00000 

k 

6l 

1 06500000 

k 

62 

i C6600000 

k 

63 

iOo 700000 

K 

64 

106800000 

K 

o5 

10690000C 

K 

66 

107000000 

K 

67 

107100000 

k 

6b 

107200000 

k 

69 

107300000 

k 

70 

107400000 

k 

71 

107^00000 

k 

72 

107o000G0 

k 

73 

107 700000 

K 

74 

10 7o0000C 

k 

75 

107900000 

k 

76 

lObOUOOOO 

k 

77 

1O81OODO0 

k 

78 

108200000 

k 

79 

I083000CO 

k 

80 

108400000 

k 

81 

lObSOOOCO 

k 

82 

LQdbOOOOO 

K 

83 

108700000 

K, 

84 

lObbOOOOO 

k 

85 

1089000C0 

k 

86 

109000DOC 

K 

87 

109100000 

k 

86 

109200000 

k 

89 

109300000 

k 

90 

I094DOOOO 

k 

91 

i 09500000 

k 

92 

109600000 

k 

93 

1O97OOUO0 

k 

94 

109800000 

K 

95 

IU9900000 

K 

96 

liOOOCOOO 

k 

97 

IIOIOOOOO 

K 

9b 

110200000 

k 

99 

110300000 

K 

100 

11U400000 

k 

101 

1L0500000 

K 

102 

11O6D00O0 

k 

103 

110700000 

K 

104 

liOBOOCOO 

K 

105 

1109000G0 

K 

106 

lllOOOOCO 

K 

107 

llllOOOOO 

K 

108 

m^ooooo 

k 

109- 

111300000 


30 



OVERLAY (JMF,2,0); D2182 

Subprogram D2182 is the controlling subprogram in overlay 2,0. Its function is to 
direct the calculation of the generalized aerodynamic forces and, optionally, of the 

flutter solutions as described in reference 2. A NAMELIST input is read. If it includes 
previously calculated Qfj values, their calculation is bypassed. Otherwise, the III 
and matrices are obtained from TAPE9. During the calculation of Qy, most of 

the intermediate matrices can be optionally printed. Finally, flutter solutions are made 
for up to 25 values of air density, and the output is printed. The flow chart for D2182 
follows. 


FLOW CHART FOR FROGRAM 02ISZ 
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( eoo;) 


/if CWPR6.NE.0 PRINT MATRICES t6l AND(7l/ 
/if CHPF8«NE.0 print ITM3 AND MATRIX 18)/ 


/If CWPRgX,NEaQ PRINT NATRIX IBX)/ 


ITH4 - M ATR IX 18X1 T IMES MATRIX III 
IC0MPUT6 QIJSTAR MATRiXt EQ. 




PUNCH QIJSTAR MATRIX 

AND MACHtKfTANLF»TANLTF,S,NMOOSt AND NSYM 
IFORMATTED FOR SUBSEQUENT NAMELIST RE API/ 


PRINT MATRIX 0*^ OEFtFCTkONg Hit 

PRINT SPAN STATIONS 
PRINT ANM MATRIX 

PRINT GENERALIZFD MASS MU MATRIX 
PRCNT MODAL FREOUCNCIES OMEGAIII 
PRINT HDOAL OAMPlNO COEFFICIENTS GlI )> 




/ PRINT QIJSTAR MATRIX / 
YES 


j 


CPHPUTE elements OF FLUTTER ' 
DETERMINANT, EQS. A12 AND A13 





/print value of alpha/ 


SPRINT - 

FLUTTER^ SOLUTION/ 

C ^310 ) 

Y 

( RETURN ) 


/print flutter^ determinant/ 

A L L EECMf COMPUTE EIGENVALUES A ND EIGENVECTORS^ 

I 

/mNT eigenvectors / 

/ NORMALIZED TO LARGEST ELEMENt/ 


(TIz*o ) 


COMPUTE DAMPING COEFFICIENT G 
AND STIFFNESS PARAMETER /V 


I 


r PRINT EICENVALUESt EIGENFREOUENC I ES Wf 
OAMPING COEFFICIENTS 0 » AND STIFFNESS 
PARAMETERS /V 

~ ^ 

r RETUKW TO MMiQ 
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oveRLAmwht 2.0 j 1 

PKUoKAM J21b2 L Z 

L)lMt^SIU^ XbTtil3(lC).XBTB2C(20>*TaiCl0( l0J*ToiC20(20JtYSP(llJ. L 3 

lXL(2CJtALt(lU.AlEIU) pM'fp 116.20). f^TK2 ( 20J .MTK3 4 20 .6 > , HTk^ it> » 6 > , L 4 

2MTR<fA(o. 11 ) ,nLtlll,6J,MTR5<6,6,H),MTR6(lli.MTK71il,LJi, MTKiilo L 6 

3.0. 10) .MTKaXto,4Li J . ITM 1( 6. 20) . ITM2( 6,6) , 1 TM3io,u .1 1) .lUtNT (6J , L b 

4ACbARl20) ,nd( 6) , YSPS( 1 II ,H IC U 1 , 20 .6 ) . HT At) UO J . Iyl»X (6) ,U^L4o) . L 7 

sTeMPU) . ol^i) .uhboA [6 I .MASS <6.6 J , AlFHAi 25),Aiol20J . ATAO120) .SMG( 6) L b 

0, fcGMto(oJ .lSPi>ail.XTft(2C),INTHACT{6.2J L S 

oLMPLbX T-^tn46.16) 16 ,6 > . GA^^MA l6 J .CAPoIopOJ ,i^iJ4o.oJ,tlt>l6) L 10 

1 , lTM4i6.*t£> J , U-fc, AD[6), STOP (6 ,6).HACT (b,6), VtLJ 6, O) ,VtC6AV L 11 

KbAL ^^TRl,MTK2,MTP^.MTR^t.MTR^A,MTR3TMTR6f?lTKr,MTKotrtTK6A. ITMIf lTh2 L 12 

lfirM3fK*MASifK,dfMACHtTSPS L 1^^ 

CJMinUN/TKA21o2/Art, XK.TANLE .TAf^lLTt ,S,NMUD£,Ni1fN L 14 

INIbGLK Cv 6P6T,C«SPI ,C WCDr.CV(PK2 .CWPR3tCWPK^,C«PKl>, L 16 

iCwPKo, OwPKb , LwPKtiX L 1<> 

C UATa TAdLL^ PU'’ integration STATIONS Aiw PACToPS PJk 10 ANu L 17 

C 20-PuiNl oAuiSTAN QUACRATURt PUk OIJ.S L lb 

oATA I XbT a 1 0 1 J ) . J = 1 1 1 *' I / ♦ 1 204o73 oE- 1 , • 67Aub3 1 71;;- 1 . #1 o02Sp22 • L 19 

U2B33U230, * '+4i!>:>o28 3 , . 5 74 ^ ?7 1 7, . 7 1 6697 70 . • o39 ? 047 d , ,9;J2 631od. L 2C 

2.9bo963d.b/ 1- 21 

JATA 1 XDTa20U J , J = 1 . I / . 343 57Q04fc- 2 , • 1601 403oc- 1 , -43ao27 6bfc-l , L 22 

1 .dD'.^lbWh-i 12 06 3405 1869731b, • 244566 bO . - 3151 ^oVo . *3 oul070 7 , L 23 

2«'»ol7jio7‘*i «:>3o2ob26. ♦6 1389293, ♦6J6fib3^4t* / :>0. # di ^02o64 , L 24 

3-a731ob9i, * 9 19 bbo49 ♦ - 956 1 172 I , . S 81 98590 , ♦99obo43C/ L 2b 

JATA (TuiClOlJ) = 1/ .333 3567 2b“l,- 74 725o7:jE-if 1095431 6, L 26 

1. l3‘*6333o. . U7 7 o 21 1.- 1 4 7 762 1 1 , . 1 34 633 36 , ♦ 10'V5^31 o, . ?*♦ 72bo 7bt-i > L 27 

2 .3333bo72b-l/ ^ 28 

JATA UblG23 U) . j = l,2"' )/ .B8''7003 5L-2,^2030C7l5t-if . 3i 3oo024L- i , L 29 

1 .41o3o3 7l- 1 , •:^09o5060F-l . ♦ 5909 72 66E-1 , • 65ot4 3 l9b-l , - 71 C4d0 5b t- i . L 30 

2. 7456o493l-1 , . 7o 376694E- I , . 7 63 7o694E - 1 , • 7Abdo493 L- 1 , ♦ /i04b055t-l, L 31 

3.tDd4‘t3i9t-if*b909726 6E-lf*5C9650o0E-l,«^lo3oi7t-if.3iJuo02‘*t-lf L 32 

4.2D300/l:^c-l,.65J7C^35E-2/ L 33 

NAMbUbf / IP2io2/NSPST ,NCOSTt«ASS. CMtOA ,0M£b 6,0 , male hA.AoPHAtVbPS, L 34 

1 ISPS.HIU.ottSPbT .Ck^SPT ,CV.CDl ,XTAB,C>yPR2,CrtPK3, L 3b 

CrtPR^ ,L«Pkd ,l^.Pk 6.CWPR8 .CwPRdXtXC eAK,Xl L,N)iP L 3b 

3fINTUr\tNPNCHw,NkjlJ,0IJ,^ACH,KTTANLE.TANLTt,S,NNuLiS,NbyM L 37 

NMUuS=NMuJb iPALn=XNSK=XK L 36 

PkInT 10 i- 

IJ POKi^AT 4 Ihl , ///22X , 87HPR OGPAf^ J2162 *♦. CbNPolt oC Nt kAL i IcJ AtKJUY L 40 

InAMIC Pu^CbS ANU SCLYE FLUTTER J bT EKMI N AHf 5 ) L 41 

C JtFAuLT VALJiS Ot CONTROL PARAMbTbRS FUR .<UAuh.AT vjKb . PKiNTiNo, L 42 

C PUNOtilNGt ANu INPUT L 43 

O^^CJI=ChPK2^LwPKj = CWPP4 = CWPP 5-CwPR6 = CwPk8^trt^'^^8X =0 L 44 

LV.SP i = ChSPsT=NPNCHt = NRQIJ = C ^ 4b 

C dtoIN INPUT L 46 

RtAu (bf2DI IoLimT L 47 

20 FURMAT UaIOI L 

PRl.a 30, iJtNT )- 

30 FUKMAT ( ///lUA.dAlO/ZI L 50 

JJ tO 11=1, 5o L 5l 

jijnn=^o.o L 51.1 

43 MASS(in=J*0 L ^2 

jO bC I 1, 1^23 £ i>3 

bJ HlUii > = 1 /7/00000000007777778 L 54 

REAJ 15, lP2id2) L 55 

L»0 P0R;-1AT (//lx.32n ECHO Of NAMELIST IP2L82 FOLLUV^S ) L do 

rtKiTt io.oO) L S7 

hKll fc (6, IP21B2) L 56 

IF ( J. bw . 1 > oG TO 1030 L 59 

rtbwiNJ 9 L 60 

RLAJ (9) i ITKEri JK.KJ) , JK=1 ,43),KJ=L,4) t 6l 

kbAj (9i I ITKi.lf4 JK.KJI ,JK==1,4 c)tKJ = 5.16J L 62 

JU 7 0 xJ = I.NMUiJS L o3 

7 3 RtAJ (9J UNHUa,KJJ,JK=1,16I L 64 

RurtlNu 9 L cb 

aKITL Io,dO) ('^ACm, K.TANLE iTANLTE , S ,NMUl>StN5Yh L 6o 

dD FUKMAl <//oH MACri=El5.8, 5X,2HK = El5.8,5Xt6HrANLb = tlb.o,5X, 7hTANLTt= L 67 

It I 5- ti, 5X, 2HS=t 1>. 8//7H Nf^00S = l 2, 5X,5 hNSYM=I 2/// j . L 6d 

C ^ L 69 

C*** lNlTlALl2t:i‘iATKiX ELEMENTS TO ZEkU L 70 

JU 9u J=l.i'<CUST L 71 


111500000 
lllbOOOOO 
111700000 
ill bOOGOC 
1119000UO 
112000000 
112100000 
U22000Q0 
112300000 
112400000 
112500000 
Il2b00000 
il27a0UU0 
112600000 
112900000 
11300000C 
1I31O03QQ 
i 13200000 
113300000 
113400000 
U3500000 
113O00000 
113700000 
1 13800000 
113900000 
LI4000000 
114100000 
114200000 
li‘t300U00 
114400000 
1L45JOOOO 

ii4ooocoo 
1 147UUOOO 
114600000 
114900CC0 
iibOOOOOU 
lUiOOGOO 
11520U00U 
11530000:' 
115400000 
libbOOUOO 
liboOOOOO 
1157U00U0 
115600000 
115900000 
116000000 
iiolOCCDO 
llo200000 
116300000 
116400000 
116500000 
liboOOOOO 

116 /ooooo 

Ll6dOOOOO 
11O900000 
117CC0C03 
117100000 
1 17200000 
11/300000 
11/400000 
11 7500000 
llVbOOOOO 
11770000U 
U780000C 
11 7900000 
1 1 dOOU jOC 
lldiOOOCO 
116200002 
118300000 
1164U0000 
1 IdbOOOCO 
lidbOOOOO 
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o o 


^0 (Mf<2 IJ ) = 

iJO 1 lU J-= 1 1 o 
Ou I LQ J J = L i i'^MUL> S 
lMJ LCJ Ji>=lthSHi>l 
iOO ITM3 iixS J *NJ J * JS J - ^ 

.110 CLNTI^Ut 

UQ I dO J = 1 1 a 
Uu 12U NJ = i«i%CUSI 
120 ITMKJtKJi^U.O 
C 

JO 130 J^= l.fv/SPi>] 

130 MT^oi jPJ -0. 0 
1^0 CLi\TiNUt 

DO 1 /U 1, o 
JU I o 0 Mk = 1 t o 
jO 1:^0 Ml=1* 10 
IjJ MTh.d<rt jf MK ,Ml) =0*0 
1^0 CLWTIfviUt 
17j CbN7l<\Ui: 

Do 190 MJ=lfM^Uub 
JO IbO WL=l.hMUJi 

IbJ 

190 CLMINUt 

o ^PA*4rtlSc: bTATiuiM^ FSTABLISHEJ 

C 

IF C LrtbPST *bw.0J GP tg ?in 
C 

0 N^.\-iT A.NoaKU 5PAN STATIONS 

C 

JO 2C0 J = 1 .NSPST 
200 ysp* j) = ySPSi 
oL TU 230 
C 

C STANuARu SPAi^ stations I?»3t7t9,0KU CtjUALLY SPAuuu) 

C 

210 Lnr-i/^LUAni\SP3T-l ) 

VSPil)=D 

uO 220 J=2 ,NSPSr 
220 YSPl j|=TiPlj-l J>uY 
230 Lui>lT I N J t 

JO :>70 1 $PAfs = l tN:>PST 
UU ‘♦20 Mjufc^ l*NMODS 
JU 290 J=lto 
UO 2 40 r.J=i,KMUiij 
243 iTM2I JtKjj=0.0 

CCMPUl t A;MU SrUKt CHOFD INTEGRALS UF LNSTAR - MTK4A 

ALL ( liPAiMJ = AjS I r ;,P ( ISPANJ J *T ANLE 
XTb(l3PANj = l,*-A[ib( YSP{ ISPAM I*TANLTk 
AVAL=aLL ( I SPANJ 
JO 233 J = ltlN,riP 

XTijI Ji^ALbi 13PAN) + XTA,<^[ J 1 S PAN ) -XLE U ;>PAN J J 

2:>0 HT A j I J J shl C 11 SPAN, J ,M00C ) 

C INTfcRPuLATt Lt JcfLECTIDN 

CALL FTluP ( AVAL ,MVAL ♦ I NTOR , NHP , X Tb ( U t hTAt>l i i i 

HLEI ISPAN,Mujl J =HVAL 

ALT^g=xT tl ISPANJ «*2-XLE{ ISPANJ **2 

XL TSk^SUK T i alt Sui 

ALTS3^AL rsk«=>3 

XLTS3- ALTSkfl-^D 

alL2-XLL I I SPAN) ♦ ^2 

IF txLLI ISPAN J GO TO 2^3 

xTUAL-AdS< XTbl I 5PA M /x LF U SP AN) ) 

xCGSh=ALUij(aTUXL^‘$CKT( XTCXL*'0^2-i . J ) 

^U TO 270 
2oJ xCu:>H=0.J 
270 CCNTlNJL 

MTK*»A( 1 , i.iPAN) =XL TSR 

i' 1TK4AU, 1 sPwNi ^.u^XTE I IS PAM ♦xLTSRF-5«XLb2»ACu:>h 
MTR‘fAl3t ISPANI=,:>*XTE( I S PAN ) ★XlT SR-- :>♦ XL b2* XCu3M 
MTk9A14,ISPmN)=*33333333^XLTS3 


L 

72 

1 IdZOOOCO 

L 

73 

11 cjtiOOOOO 

L 

74 

116900000 

L 

7S 

i 1 sQooooa 

L 

70 

llSiOOOOO 

L 

77 

119200000 

L 

7d 

i 1 9300000 

L 

79 

119400000 

L 

‘ tJO 

i 19SOOOOO 

L 

til 

1 1 SdOOGu J 

L 

b2 

119700000 

L 

o3 

1 ISbCOOOO 

L 

B4 

119400000 

L 

as 

120000000 

L 

tit 

120100000 

L 

o7 

120200000 

L 

as 

120300000 

L 

as 

U0400COO 

L 

9C 

120SOOOOC' 

L 

91 

12U6000UU 

L 

92 

1^.0700000 

L 

9? 

UOBOOOOO 

L 

94 

1^.0900000 

L 

93 

121000000 

L 

96 

121LU0000 

L 

97 

12U00C00 

L 

9b 

121300000 

L 

99 

1214000CC 

L 

100 

121S0000C 

L 

101 

121jOOuuO 

L 

102 

121700000 

L 

103 

UibOOOjU 

L 

104 

12190C0C0 

L 

lOS 

122000000 

L 

106 

1221000U0 

L 

107 

122200000 

L 

108 

122300000 

L 

109 

12ii40UOOO 

L 

11? 

l2nSOCOOC 

L 

Ill 

122o00000 

L 

112 

1227000CO 

L 

113 

122600000 

L 

114 

122900000 

L 

US 

1230COOOO 

L 

lie 

12310000C 

L 

117 

12320000J 

L 

lib 

123300000 

L 

119 

123400000 

L 

120 

123S00000 

L 

121 

12:>oOCCOO 

L 

i22 

123 700000 

L 

123 

I23o00000 

L 

124 

I239000CO 

L 

12S 

124COOOOJ 

L 

L2o 

124100000 

L 

127 

1^:4200000 

L 

126 

124300000 

L 

129 

12440COOO 

L 

130 

1241JOOOOO 

L 

13i 

1<:^600000 

L 

132 

124700000 

L 

133 

124600000 

L 

134 

124900000 

L 

13S 

i2SOOOOOO 

L 

13j 

i2bl(JOOOO 

L 

137 

12:^200000 

L 

136 

lc>300000 

L 

139 

12^400000 

L 

140 

12:>dCOOOO 

L 

141 

12S60000 J 

L 

142 

l^bVOOOOO 

L 

143 

12SoOOOOO 
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MTK4A4 3. ISFAN)=.2 5*XTE M SP AN) ♦XL TS3*. 1 2 5»ALfc2«AT £4 ISPAN J ♦XLTSK- . 12 L 

144 

125900000 


It»^ALE2^*2*XCUi»M 

L 

143 

12oOOOOOO 


MTR4AI f>,l SPAN J*. 2 ♦XLTS5^- 3333333 3 + XLE2*XLTi3 

L 

146 

12bI00000 


UO 280 NN^I.C 

L 

147 

126200000 

Z60 

HTk4(NN*MjDb J = MTR4AINN, I SPAN J *HL E llSPANt MUu£ i 

L 

146 

126300000 

c 


L 

149 

126400000 

c 

COMPUTt ANO STOKt CHORD INTEGRATING MATRIX - MTK2 

L 

150 

126300000 

c 


L 

151 

12CO00000 


IF ICk^COI *EiJ.O) 1.0 TO 300 

L 

152 

126700000 

c 


L 

153 

1Z6600000 

c 

NON-STANOARU LHUKDhISE INTEGRATION 

L 

154 

126900000 

c 


L 

155 

127000000 


DO 290 J=1.NCUST 

L 

156 

L27100000 


MTK2U J = XlCIJ>*(XTeiISPAN)-XLE nSPANiJ 

L 

157 

127200000 

2^0 

XC U)=^AL£t ISPANi 4XCRAR IJ <XT£( I S P ANi -XL E UsPAN Ji 

L 

156 

127300000 


GO Tb 340 

L 

159 

127400000 

C 


L 

160 

127500000 

c 

bTANDAku CHOKJHisF INTEGRATION (10 OR 20-PuiM GAUSSi 

L 

161 

127600000 

c 


L 

162 

127700000 

300 

IF ( NLUST.Ew-lOI GC TO 320 

L 

16 3 

127600000 


DO 310 J=1,nCDST 

L 

164 

127900000 


MTR2{j)^rblC20U)*(XTE( ISPAN)-XLE( I SPAN I > 

L 1 

L65 

12800 OOCO 

3iU 

XLi JJ^XLt UbPANJ +XBT62 0( J)* I XTEI I S F AN)-XLt ( iSPANi J 

L 

166 

126100000 


GO TO 340 

L 

lo7 

126200000 

320 

DU 330 J-1,NCDST 

L 

168 

126300000 


MTR2U>^TdIC10l jJ*{XTE(I 5PAN)-XLE i ISPANI > 

L : 

l69 

126400000 

33Q 

XC ( J) = XLt ( ISPANI ♦xaiBl^ U)*< XT£( I SPAN ) -XLt U bPAN JI 

L 

170 

126500000 

34U 

CONTINUE 

L 

171 

1266QOOOO 


DO 410 XChU=l.NCDST 

L 

172 

126700000 

C 


L 

173 

126600000 

c 

compute and bTOKc DEFLECTIONS H-Hie(MD0E U “ MlK3 

L 

174 

126900000 

c 


L 

175 

i29UOOOUO 


IF (ChOuUNt-Oi ^0 TO 36^^ 

L 

176 

129100000 

c 


L 

177 

129200000 

c 

standard iNTbkpLiLATION OF DEFLECTILNS 

L 

178 

129300000 

c 


L 

179 

129400000 


DO 3!?0 J = 1*NHP 

L 

160 

129500000 


XTb( J) =XLL ( ISPAi'^J +XT A8 ( J ) ♦ < XTE < 1 SPAN l-XLE ( ISPANl J 

L 

I8I 

129600000 

3f>0 

HTAdUi=Hibl ISPAN, J .MODE J 

L 

162 

129700000 


XVAL=xCt KCHD) 

L 

163 

1296000UO 

C 

INTERPOLATE lEFLcCTIONS 

L 

184 

129900000 


call FTLJP i aval fHVALf INTCR ,NHP,XTfi m , HTAbd i i 

L 

165 

130000000 


MIRilKCHU. ML JU = hV AL-HLE ( ISP AN, MODE J 

L 

166 

130100000 


GO TO 370 

L 

167 

130200000 

c 


L 

166 

1303UUU00 

c 

NGN'ST ANuAkD, Uit CEFLECTIONS NIC h/Q INTEKPucAT I JN 

L 

169 

130400000 

c 


L 

190 

130500000 

3b0 

MTk 3 IKCHJ .NUuU=Hl C ( IS PAN# KCHD,M JOE) -HLE( iSPANfMuuU 

L 191 

13U60000CI 

370 

CONTINUE 

L 

192 

1307000C0 


SgkTkAO^SuKT I XC(xCHD)*«2-XLE n SPAM«*2) 

L 192.1 

1308CD000 


xCN=1.0 

L 192.2 

130900000 


DU 400 N- 1 1 cj 

L 

193 

I3LOOOOOO 


IF IN.LU.IJ uU Tu 380 

L 

194 

131100000 


IF IN.tw.2J bO Tu 390 

L 

195 

131200000 


MTRl iN.KChJ J-XCN*SORTR AD 

L 19o 

i;>I30C000 


XLN=XCN*XL INLHDJ 

L 196.1 

131400000 


uO TC hOO 

L 

197 

131500000 

3d0 

IF (bUKTKAD.NE. O.^J CO TO 382 

L 197.1 

131600000 


MTRl IN.XCHDJ =0.0 

L 197.2 

131700000 


GO Tb ^30 

L 197.3 

131600000 

3d2 

MTRl IN# kCHDJ=aGI kCFDJ /SOPTRAD 

L 197.4 

131900000 


iyO TU 400 

L 

199 

132000000 

390 

MTRl (N .^CHDJ=XC^ xCFDJ^MTR !< N-l,KChD) 

L 

20Q 

132100000 

400 

CGNTINUi: 

L 

201 

132200000 

410 

CGNTlNUt 

L 

202 

132300000 

420 

buNTINUt 

L 

203 

132400000 

C 

end UF MUUAl JG LOOP **« 

L 

204 

132500000 

C 


' L 

205 

132O00000 

C 

MATkIX uPtKMTIUi\ FOR CHGRDVvISE INTEGRATION 

L 

206 

1:»27000CC 

C 


L 

207 

132800000 

C 


L 

206 

1329C0000 
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L 

CLMPUTt AND 6TUKt: CHORO INTEGRALS - MTR6 

L 209 

133000000 

C 


L 210 

133100000 


DU 430 J=I»o 

L 2L1 

133200000 


UU 43u JJ^ltNCuST 

L 212 

133300000 

430 

I TMl i J , JJ J =MTK1 IJt JJ) ^MTR2UJI 

L 213 

133400QCO 


DO 440 J - 1 f o 

L 214 

13350000J 


UU ^40 Jj=ltNMUUS 

L 215 

133O0C000 


DG 440 JJJ=ltNCUST 

L 216 

133 700000 

440 

I TM21 J * J J> = I TMil J , JJJ) ♦MTR3( J J J, J J ) 4 1 TM2 1 J , JO ) 

L 217 

133800000 


DU 430 J=l,o 

L 216 

133900000 


DO ‘♦DO JJ=l.NMODS 

L 219 

134000000 

4t>0 

rtTRDiJ J* Jf I dPAN) -1 TM2t J , JJ)<-MTR4 « J , JJ> 

L 220 

134100000 


Ih iL^PK2.tO-0i GO TG 480 

L 221 

134200000 


k IT £ 1 tj p 460 ) 

L 222 

134300000 

4b0 

FUKMAT ( //10A*21nMATRIX PTR! - BV RCwS ) 

L 223 

134400000 


rtRITfc (o*iOQO> UMTRK JK pKJ) ,KJ= ItNCDSTJ pJk^lfuj 

L 22 4 

134500000 


wRlTt lo»470) 

L 22 5 

134600000 

470 

FORMAT W/lOA* 3DHMATPIX ^*TR2* DIAGONAL fcLtMLNTS i 

L 226 

134/00000 


wRITt UflOOC) (MTR?UK), JK = 1,NC0ST} 

L 227 

134800000 

460 

IF ILwPK3.£o-0J TO 500 

L 226 

134900000 


wRlTfc (o,490J 

L 229 

135000000 

490 

FORMAT l//10Af34hDEFLECTICNS H-HL£ , MTR3 t dY KuwS J 

L 2 3C 

135100000 


WRITE lo.lOCG) i iMTP3UK,KJ) ♦ J 1 1 NCDST J t K J- 1 ♦ NMLUSi 

L 231 

1352OD00U 

■jOO 

IF lUPK4.fcw.J) jO TO 52"' 

L 232 

13 5300000 


WRITE loibLO) 

L 233 

13D400000 


FORMAT ( //10a»33hCHORU INTEGRALS LF LNSTAk - MTK4A/i 

L 234 

135500000 


WkITfc (utiOOO) IMTR4AIJK,ISPAN)t JK=LtO) 

L 235 

135oOQOGO 

i>2Q 

IF ICwPR4.LU*0) go TO 560 

L 236 

135/OQOOO 


WRITE (o.630i 

L 237 

135600000 

p30 

FORMAT (//lOXt^lhWATRlX MTR4 - BY ROWS ) 

L 236 

135900000 


«kUt (otluOOJ UMTR4I JKpKJ) ,KJ=l,NMGD$J,jK=i,6> 

L 239 

136000000 


wRITt (6,D40) 

L 240 

136100000 

5^*0 

FUKMAT 1 //10x*33hMaTRi X ITMl * MTR1«MTR2 , uY RUwS J 

L 241 

136200000 


WRITE I OtiOOCJ UTTMl ( JK ,KJJ ,KJ- l,NCDST ) » JK=l,oJ 

L 242 

136300000 


wRlTE lo»550i 

L 243 

136400000 


FORMAT I //1CX»33HMATRI X ITR2 = irMl*MTR3 t oV RUwS J 

L 2 44 

13O500000 


wRITk lotlCOCJ i i I tm?{ JK »KJ ) ,K J= ItNI'^OOSI f JK=Lf u) 

L 245 

136600000 


CUNT INUh 

L 246 

136700000 

■>10 

COwT iNUb 

L 247 

136600000 

c 

END uF $FAN-STAT1CN OC LOOP *** 

L 248 

136900000 


IF (CwPK4.Ey.0J GO TO 590 

L 249 

137000000 


WRITE (o*b60i 

L 250 

137100000 

>b0 

FORMAT t//10X»6lHJ£FLECTI0NS HLE OF MTR4, bY 6PA»t* STA., MOuE 1 FIkS L 25L 

137200000 

IT ) 

L 251.1 

13/300000 


WRITE (bt lOOC J UHLFUK.KJ I » JK==1 , NSPST ) ♦ K J = 1 , NMUU6 ) 

L 252 

137400000 

t»90 

CQNT INUE 

L 253 

137500000 

C 


L 254 

137600000 

C 

COMPUTE AND STOKE SPAN INTEGRATING MATRIX ^ hTRb 

L 255 

137700000 

c 


L 256 

i37o00000 


IF ICwSPi .t j-OJ jO TO 61C 

L 257 

137900000 

c 


L 25b 

138000000 

L 

nln-jtandakj spanwise integration 

L 259 

136100000 

C 


L 26C 

1362000C0 


DO oOC JS=1.NSPST 

L 261 

136300000 

600 

MTkoI JSJ -i SPSUiJ 

L 262 

136400000 


GO TO o30 

L 263 

138500000 

C 


L 264 

13bu00000 

C 

standard spanwIse integration 

L 265 

136700000 

C 


L 266 

136800000 

ol 0 

yuANT=S/FLuATI B^^iNSPST-l J ) 

L 267 

136900000 


MTKoU ) = UJANT 

L 2 6b 

139000000 


NSM=NSPST-l 

L 269 

139100000 


Ou 62U JS^2.NSM.2 

L 270 

139200000 


MTRt>USl ^‘♦.O^JUAnT 

L 271 

139300000 

o20 

MTR6US^I)=2.0*UUANT 

L 272 

1394O0C0D 


MTRblNSPSTJ-gjANT 

L 272.1 

1395000UU 

C 


L 2 73 

139600000 

630 

CONTINUE 

L 274 

139700000 


DO 64C JS=1 pNsPST 

L 275 

L39buOOOO 


«TK7(JSf IJ^I.O 

L 276 

139900000 


DO o80 JS=1*NSPST 

L 277 

140000000 


IF INSVM.Eu.U GO TO 660 

L 278 

I4C100000 
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00 650 JP-I, 10 

L 

279 

140200000 


JSPs2«i jp-1 j i-1 

L 

260 

140300000 

ob3 

MTR7 (JS* JPi =VSP( JS )**JSP 

L 

261 

14040000U 


00 TO 6b0 

L 

282 

140500000 

663 

OU o70 JP-2tlO 

L 

283 

140600000 


J6P==2» ( jP-21 

L 

264 

140700000 

o70 

MTR7ijS. Jpj = JS)**JSP 

L 

265 

140600000 

6b3 

CLNTlNUt 

L 

286 

140900000 


DO 700 NSJ=1»6 

L 

287 

141000000 


DO 700 WjJ-l*NrtOL)S 

L 

288 

141100000 


DU o90 J6=1,N6PST 

L 

289 

141200000 

t>yo 

ITMBiNiiJ. NJj , JSi=f^^TR5( NJ JtNSJ. JS) ♦^^TR6t JS> 

,L 290 

141300000 

703 

CL^NTlNUt 

L 

291 

141400000 

C 


L 

292 

141500000 

C 

-■MATRIX OPfckATION FOR SPANWISE INTEGRATION 

L 

293 

141600000 

c 


L 

294 

141700000 

c 


L 

295 

I4ld00000 

c 

COMPUTE ANu STOAt OIJ-STAP MATRIX - MTR8 

L 

296 

141900000 

c 


L 

297 

142000000 


UU 74C MJ=l,o 

L 

298 

142100000 


DO 730 MK=1,NMuU6 

L 

299 

142200000 


DO 720 ML=1* 10 

L 

300 

142300000 


DO 710 MM=l»N6P$T 

L 

301 

142400000 

no 

MTRolMJ*MX*MLJ =1TM3 (MJ,MK,PM)#MTR7 (MMt ML) ♦'Mf Ko ( Mj » MK , ML i 

L 

302 

1425000UU 

720 

CLNTlNUt 

L 

303 

142600000 

730 

CUNTiNUt 

L 

304 

142700000 

743 

CQNTINUc 

L 

305 

142800000 

C 


L 

30 6 

142900000 

C 

COMPUTE AND STuKt CIJ-STAP. REARkANGEO MTK6X 

L 

307 

14300000U 

c 


L 

308 

143100000 


L)u 770 Mj = l»NMUDS 

L 

309 

143200000 


NN-0 

1 

310 

143300000 


DO 7oO MX^l.o 

L 

311 

143400000 


ih lMK,twn»LK*MK*EO-?l NTM=13 

L 

312 

143500000 


IE lMx.EW-2*0K-Mx.EQ-4) NTP-8 

L 

313 

143600000 


IE lMK.tg.5J 

L 

314 

143700000 


if iMK.tU.ol NTM=4 

L 

315 

143800000 


UU 750 ML=LtXTM 

L 

316 

143900000 



L 

317 

144000000 

7I>0 

MTRdXi MJ »NN) =MTK6(MK*M JtVL 1 

L 

318 

144100000 

/60 

CUNT INUt 

L 

319 

144200000 

770 

CUNTINUt 

L 

32C 

144300000 


IF (CrtPKi-Eg.O) uO TO 800 

L 

321 

i-,4400000 


rtrt I Tfc ( d T 760 J 

L 

322 

144500000 

7oO 

FORMAT i //iUXt22rlCHORO INTEGRALS - MTR5/ 6A, 15hN=0 UN lNSTaKI, L 323 

144O0000U 


IIU. 3HN-U l/At3HN-2, 17Xf3HN = i» 17X,3HN = 4, 17X,3hN=3 ) 

L 323.1 

144700000 


DO 790 MJ = l*NM(JJ6 

L 

324 

144600000 


print 7b 2 ♦ Mj 

L324.1 

144900000 

7a2 

FORMAT (5x# :>hMuJE ♦!?) 

L324.2 

145000000 

793 

WRITE (o.lOOCJ ii MTR5< MJ , JK »KJ) » JK- 1 , 6 ) , K J- 1 , F 5T J 

L 

325 

1451UOOOO 

600 

IF iCriPKo.Eg.O) bO TO 8?C 

L 

326 

145200000 


WRITE iotDiOi 

L 

32 7 

145300000 

blO 

FORMAT < //10a*47mSPAN lNT£GfiATlNG MATRIX - MTKo lUiAGONAL 

TERMS)/) L 

328 

145400000 


white iOflOOC) UiTR6UK)t JK=L,NSPST) 

L 3 29 

145500000 


wKiTt (dt620J 

L 

330 

145600000 

620 

FORMAT J// lOXt 23HMATRIX - MTR? - BY ROWS /I 

L 331 

1457COOOO 


white (btlOOCJ UMTR7UKtKJ)«Kj=l»10).JK = i,N6PiT J 

L 

332 

145800000 

633 

IF U.wPKo.tU.0) GO TO 

L 

333 

145900000 


WRITE Ib,64JJ 

L 

334 

146000000 

6^0 

FORMAT i//10A*36MMATPIX ITM3 - MTK5*MTK6 t or COLUMNS ) 

L 335 

146100000 


JU a50 JK=1»NmOU6 

L 

336 

l*t6200000 

6dO 

WRITE I6*100CJ 1 l!TM3( JO , JK , JM ) , J J=l ,6 ) t JM- 1 , NSP ST J 

L 

337 

i463000u0 


WRITE 1 6 1 ooO i 

L 

338 

146400000 

660 

format l//10A,50HSPAN INTFGRALSr MTRBtBY ROwS.oO tLEMtMS 

EACH MODE L 

339 

146500000 


1 / J 

L 

340 

146600000 


00 670 Mj^ltKMOuS 

L 

341 

146700000 

673 

WRITE I 0,1300) ( iMTR8( JK , M J , KJ ) , K J= 1 , 10 ) # JK= 1 , 6 ) 

L 

342 

146800000 

ddO 

IF ICwPRoa. tO-Oi GP TO 90C 

L 

343 

I4o900000 


WRITE J6,690) 

L 

344 

147000000 

690 

FORMAT i //IO a, oZrtSPAN INTEGRALS RE ARRANGED, Ml K6A t KOwS OF 

•♦6 ELEME L 

345 

1471000CCI 


ints each mode/) 

L 

346 

147200000 


WRITE (o,ICJC) 1 l«TR8XUK,KJ l,KJ = l,46) , J K=^i , NMUD $) 

L 

347 

W7300000 
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9(J0 

CU(\iTlNUL 


L 

34b 

147400000 


DO 920 KJ=l*NrtODS 


L 

349 

147500000 


DO 910 HL=lt4o 


L 

350 

147600000 

910 

I TM4 (Mj 1 MLi =0.0 


L 

351 

1477O0UU0 

920 

CUiMT INUc 


L 

352 

147600000 


OU 940 MJ = l,NMUJi> 


L 

353 

147900000 


DO 9^0 ML= !♦ lo 


L 

354 

146000000 


OD 930 MM=1*4o 


L 

355 

148100000 

9i0 

lTM44MJ.MLJ = MTHoxtf'J»MMJ =>T RE Y ( ^ ) ♦! T M4 IMJ , ML J 


L 

356 

148200000 

940 

LQNT INUt 


L 

357 

146300000 

C 



L 

358 

148400000 

C 

MATkIa uPtRATlCu FOP. GENERALIZED AERQOYNAhlL I'QKL.hS 


L 

359 

1^6500000 


DO 960 MJ=l,NMUDi 


L 

360 

1 '^8600000 


DU 960 Ml = l,<\MUDS 


L 

361 

148700000 


DU 950 Mi. = X*l6 


L 

362 

148800000 

950 

Ol 0 IMJiMl ) = ITM4(MJ»ML ) ♦ANMlMLi MI I J 


L 

363 

148900000 

9o0 

CLNT INUE 


L 

364 

149000000 

C 

gIJ-6TAR KATKiX IS COMPLETE 


L 

365 

149100000 


^RITE lo,970J 


L 

3o6 

149200000 

9^0 

FuKHaT 1 //ivJA*77riI NPUT deflections H SUB i» IN UKUtM 

L.E* 

TO T-t-f L 

367 

149300000 


1 INdUAKD FIRST# MODE 1 FIPST/J 


L 

36B 

1494000DO 


DO 990 Mj=l,NMJDS 


L 

369 

149500000 


DU ^60 M I =1 f NSPSr 


L 

370 

14960UOOO 

9dO 

WHITE (6»100lI inlCtW! ,J=i,NHP) 


L 

371 

14970000C 

990 

LLNTINUL 


L 

372 

149800000 

1000 

FURMAT (ol5X,El5.8n 


L 

373 

U99000CO 


WRITE lUflQlQJ 


L 

374 

150000000 

1010 

FuKMAT * //10 a,13HSPAN STATI0N$/| 


L 

375 

15O1OUOO0 


wKITb tOfiOCC) t YSP(J>TJ=ltNSPSTl 


L 

376 

1 50eOOOOC 


WRITE ( 6 ♦ 1 020 J 


L 

377 

15J300000 

1020 

FURMAT t//10A#4oHANM MATRIX* BY MODAL COLUMNS UF LOMPlEX ELEMENTS/ E 

378 

1504D000O 


15X#4HREAL# loA#4rilMAG* l6X#4HREALt l6X , 4H I M Ay » luA , 4fiK EAL 

# 1 6X f 

4HIMaG) L 

379 

150500000 


white Id* lOOOi j lANMIJKfKOif JK— 1 #16) tKJ = l« i^Mlj U S i 


L 

380 

150600000 

1030 

Cunt INU b 


L 

381 

15070000C 


AkITE (o#l04G) 


L 

382 

150800000 

1040 

FORMAT 1//10a#10mMIJ MATRIX/) 


L 

383 

150900000 


WRITE <6tlU00) MMASSIOK ,K0) ,KJ-l,NMOUS) # JK=l»NHUDSi 


L 

384 

151000000 


WHITE (6,1050) 


L 

385 

15UOOOOC 

1050 

FORMAT l//10A,17riM0DAL FPEOUENCl ES/J 


L 

386 

151400000 


WRITE (bflOOC) (UMEGA( J) ,J = 1 #NMUUS J 


L 

387 

151300000 


WRITE (6,10o0J 


L 

388 

151400000 

lOoO 

FORMAT 1 //L0A,31HM00AL OAKPING COEFFICIENIS 0(1)/) 


L 

389 

151500000 


WRITE (o.lDOOJ ( u( J ) # J=1 ,NMOOS ) 


L 

390 

151600000 


WRITE ( 6, 1070) 


L 

391 

151700000 

10/0 

FORMAT (//10a# 26H GIJ MATRIX, By CCMPLEX kUw;>/ ) 


L 

392 

l5l800UOU 


WRITE (Of 1000 J 1 ( UlJ( JK,KJ ) ,KJ = 1 ,RMGUSJ • Jh = 1 iNMUDS) 


L 

393 

151900000 


IF (NPNCHy.NE.l) GO TO 1102 


L 394 

L52000000 


PUNCH 10o0» MACHfK ,TANLE,TANLTE, S ,NMODSiNSYM 


L 

395 

152100000 

lOoO 

FURMAT ilA,5HHACil = E 15 • 8 1 , K = E 1 f>* B , 7H , TANLE = E 1 5. d , IH , 

/I Af 7MTANLTt= L 

396 

152200000 


1 £15* b, 3H, s=El5.o, 7H,NMODS= !2 ,6H, NSYM=12, IH, ) 


L 

397 

152300000 


DU 1100 Mi=l,NMUUS 


E : 

398 

154400000 


DU 1 100 MJ=1 ,NMUDS ,2 


L 399 

15^500000 


kWU-REA4. 4wiJ( HJ,MI )) $ RQ UP 1= REAL ( Q I d ( MJf 1 , Mi i ) 


L 400 

152600000 


AIUlJ = AlMAG(wl JlHJfMl) ) $ AIQI J Pl = A1 MAG ( wi J ( M J ♦•1 # MI i J 

L 401 

152700000 


PUNCH 10^0, Mj#MltRQlJ,MQ!J,KQIJPl,AIOIJPi 


L 402 

152800000 

1090 

FORMAT ( 5h gU( 1 1 1 » IH, 1 1 ,2H) = , 2( IH ( ,E15. b,ih, ft l5.d#2H) •) ) 

L 403 

152900000 

IlOO 

CUsiTiNUE 


L 

404 

133000000 

c 



L 

4C5 

153100000 

c 

btblN FLUTTER SQ.LUTICN 


L 

406 

153200000 

c 



L 

407 

1533QOOOC 

1102 

IF (K.NE.O.J GO TO 1120 


L 408 

153400000 


PRINT 1110 


L 

409 

153500000 

1110 

FORMAT (//33h K=^', NO FLUTTER SOLUTICN ♦♦»//) 


L 

410 

153600000 


bO TU 131U 


L 

411 

153700000 

1120 

JQ 13D0 NaF=1,NALPH4 


L 

412 

15380000C 


UO 1140 M1=1,NMJdS 


L 

413 

153900000 
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UiO 

L143 

llbO 


U70 

lu:) 

ii^a 

c 

uoo 

U13 

1233 

1240 

Ui>D 

12o0 

12 70 

I2o0 

12:^0 

1300 

1310 


KPT= UOMt:bb/CrthGAtMIM*^2/tMASSlMl,Mli*4 l.+utrtl l**2J > ) 

X1PT=-G(rtl i^KPT 
GAM**1AAMli =CMPLAl KPTtXiPT ) 

UiJ 1130 MJ = l,NMuOS 

CAPG(hUrtJJ = bAMMA( MiMiMASSAMI , j ) +ALPHA ( N AP J ♦ t w IJ AM 1 1 M J ) /K=>^2 J i 
CLNTINUt 

V^RITt ALPHAANAF) 

FClK*MAT I/// / 10a, oHALPHA= Elb*8) 
riRlTb AotllbO) 

PUKrtAT A//10A,3c»rtFLUTTER DETERMINANT, BT DUMFlEa KUwi/ J 
WRITE lotiOOCj UC APGl JK ,K J ) ,K J= 1 *NMCDS ) • JN=i * NM0Ui>) 

NHAA=^6 

LsiTHAbTI l,li^lNTHACTl2,l )=NMnDS 

CALL fcECM KAPG*ElG,VFC,STaR,HACT,CFE,AOf lNThACT,NMAAJ 
IF UNTHACTlLtli .GE.NMODS) GO TO IIBC 
WRITE (6,lI7Ui 

FORMAT (/oOH NCN-CCNVE RGENCE RETURN FRUM i>UokuUT INE **** 

DO TO 130J 

PRINT 1190 , . ^ 

FORMAT ( / bA,9brtElGENVECTDPS IN SAME ORDER A3 IHfc E 1 uEN V mlUE 3* • • fcAC 
LH tlGENVcCTUK NJRMALI2ED TO ITS LARGEST EtEMENTi 
UU 1230 LVEC^l,i'*MDOS 

KbNUKMALUE H^EnVECTHP, MAKE LARGEST ELEMENT ♦ 0-Ul 

JSAVE=i 

JU 1200 J-2tNM0US , . . 

it nV£C( J5 AvE,LVECJ *C0NJG (VBC( JSAVE,LVECn J*lT. IVECtJrLVtCl^CGNJG 

u veciJ,L»/LCn i J jsave= j 

CLNMNUt 

VbCSAV-VtU J3AVL, LVEC) 

DU 1210 J=l,NM0U3 

VEC IJ, L VELi =VtC( Jt IVEC J / VfcCS AV 

COnT INUE 

print 1220, LVbCflVECt I,LVFC),I=l,NMOOS) 

FORMAT I 14H 1 1 ocNVECTOR ( , 1 2 , 1 H) / 1 6E 21 . 8 ) / ) 

CUjmT INUt 

OU 12t>0 J1 TL^l ,NmODS 
eiGcNK = KLALlbIU( J ITEM 
t IGbNl =A1 MAG i b 10 IJ TTE n 
IF (bIbENK.GT,Oi GO 1240 

;,MG{ JlTbi = bUhEblJ I TE}^KB(J!TE) =1 777000000000077 77/ /8 
GU TU 1^30 

SMbi JiTb)=tlGENl/EIGCNR 
EJMEGi J 1 T bJ=LME^D/SURT ( E IGENR) 

K8I JiTEi =K^uMi:Gb/ECMEG (J ITF ) 

CONTINUE 
WkITb I6,12 oOJ 

FORMAT 1 //13A.41HEIGENVALUFS, RE m , IM (1 ) ♦ kE I , IM 42 ) t tTC ,/ ) 
rtKiTb 4o,100CJ IcIGl J ) ,J=ltNMJDS) 

WRITE ib,L27C) 

FORMAT (// IOa.IoHEIGENFREQUENCIES / ) 

WRITE 4o,lQ0Cl Ib'lMEGUi tJ = l tNMOOSI 
wRIT£(o*12o0i 

FURMAT 4 / / IOa ♦ 22HDAMP ING fCEFFIClENTS G/> 
write iwilOOOl 4 SMG( J I , J=l ►NMUOS ) 

WRITE 4ofl29Gi 

format i // lUXt 20hSTIFFNESS PARAMETERS/) 

WRITE lo.lOOOJ iKBi Jl , J=I tNMGOSJ 

CUNT INUE 

KETuRi^ 

END 


L 

414 

134000000 

L 

4li> 

154100000 

L 

41 0 

134200000 

L 

417 

134300000 

L 

41B 

154400000 

L 

419 

134500000 

L 

420 

154600000 

L 

421 

154700000 

L 

42 2 

13480OO0C 

L 

423 

134900000 

L 

424 

153000000 

L 

423 

155100000 

L 

42o 

153200000 

L 

427 

133300000 

L 

42d 

133400000 

L 

429 

135300000 

L 

430 

133600000 

L 

431 

135700000 

L 

432 

135800000 

L 

43 3 

1339000C0 

L 

434 

i5o000000 

L 

433 

156100000 

L 

43o 

156200000 

L 

43/ 

136300000 

L 

438 

136400000 

L 

439 

136300000 

L 

440 

13O600000 

L 

441 

136700000 

L 

442 

136800000 

L 

443 

136900000 

L 

444 

137000000 

L 

443 

157100000 

L 

446 

137200000 

L 

447 

137300000 

L 

448 

1374COOCO 

L 

449 

157500000 

L 

45 0 

13760000U 

L 

431 

157700000 

L 

432 

13/800000 

L 

453 

157900000 

L 

434 

15O000000 

L 

435 

158100000 

L 

45 6 

15820QO0U 

L 

457 

158300000 

L 

438 

138400000 

L 

439 

138300000 

L 

4toO 

13bo000CC 

L 

461 

138 700000 

L 

462 

138800COO 

L 

463 

138900000 

L 

464 

139000000 

L 

4o3 

139100000 

L 

466 

139200000 

L 

467 

159300000 

L 

466 

159400000 

L 

469 

159300000 

L 

470 

I59o00000 

L 

471 

139700000 

L 

472 

159o0C000 

L 

47B 

139900C00 

L 

474- 

160000000 
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Description of External File TAPE 9 

The file TAPE9 is written in Overlay 1,0 and read in Overlay 2,0. As shown on 
the flow charts, it is handled as follows, for KPROG = 0: 

(1) Rewound in subprogram D2181. 

(2) In subprogram MTXMPY the first four columns of the matrix in are written 
as record no. 1, and columns 5 to 16 are written as record no. 2. 

(3) In subprogram PRT2 the matrix of a^^ is written as record no. 3. 

(4) In subprogram D2182, TAPE9 is rewound, the matrix m is read as the first 
two records, and the matrix ajj^i is read as record no. 3. TAPE 9 is then rewound. 

For KPROG 0 TAPE 9 is not used. 


USAGE 

Program Information and Input Deck Arrangement 

The program D2180 is written in the FORTRAN IV language for the Control Data 
6000 series digital computer under the SCOPE 3.0 operating system. To compile, load, 
and run requires a field length of 63000 (octal) locations with variables dimensioned as 
follows: 48 control points, 6 vibration modes, 25 air diensity values, and 220 points on 
the half- span wing for modal deflections. An example of a central processor time to 
compile and run is 166 seconds on a 6400 for a case with 48 control points, 3 modes, 

21 density values, and 110 deflection points. The following sketch shows the input deck 
arrangement 


i 


42 











Input Description 


The input quantities are loaded by using the FORTRAN IV NAMELIST, with the 
exception of the IDENT cards that are in 8A10 format. The input symbols are as fol- 
lows: ^Recall that all NAMELIST input begins in columns 2 to 80. Also, see descrip- 
tion below of "standard integration" for Qy values^. 


$IP2180 

KPROG 

$END 

IDENT 

$IP2181 

DEL 

DHX 

H 

HMLTPLY 

NMODE 

NSYM 


NXY 

'S 

TANLE 

TANLTE 

WT 


XK 

XM 

X 

Y 

$END 

IDENT 


control word: if less than 1, execute both programs D2181 and D2182; 
if equal to 1, execute only D2182. 

identification information, one card in 8A10 format 
(equation numbers below refer to the number given in ref. 1). 
nominal value of width 6 of subregion m. Default value = 0.04, 
real part of downwash ratio on left-hand side of equation (18), 

-ahjy'Sx, one value per control point 
imaginary part of downwash ratio on left-hand side of equation (18) 
(ref. 1) not including k, -2h, one value per control point 
optional two-word array that if nonzero multiplies: HMLTPLY(l) 
times DHX, HMLTPLY(2) times H. Default values are 0.0. 
number of modes of downwash (maximum of 6) 
control word for spanwise symmetry in lift series A„* of equa- 
tion (13b) of reference 1 and A^ of equation (14) of reference 1: 

0 for spanwise antisymmetry (odd powers on 77*^), not 0 for span- 
wise symmetry (even powers on 77*^). Default value = 1. 
number of downwash control or collocation points, maximum of 48 
ratio of wing semispan to root semichord 
tangent of leading- edge sweep angle, positive for sweepback 
tangent of trailing- edge sweep angle, positive for sweepback 
array of weights used in the least- square -error solution of the simul- 
taneous equations (32). WT(1) = -1.0 causes all equal weights. 
Default value = -1.0. 
reduced frequency k 
Mach number of undisturbed flow 
array of x-coordinates of control points 
array of y-coordinates of control points 


identification information, one card in 8A10 format, may differ from 
preceding IDENT 
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$IP2182 (equation numbers below refer to the numbers given in ref. 2) 

ALPHA array of values of a in flutter equation (A 12), maximum of 25. 

ALPHA and MASS must be in consistent units. 

CWCDI control word for chordwise integration in equation (lib) of reference 2: 

0 for standard, #0 for nonstandard, integration. Default value = 0. 

control words for printing matrices (2), (3), (5), (6), (8), 
and (8x), respectively of equations (14), (16), and (18): 
*0, do not print; ^0, print. Default value = 0. 

CWPR4 control word for printing matrix (4), its integrals of 

hj^x^g,y^^: sO, do not print; ^^0, print. Default value = 0. 

CWSPI control word for spanwise integration of equations (lib) and (16): =0 

for standard Integration, #0 for nonstandard. Default value = 0. 

CWSPST control word for span stations at which chordwise integrations are done: 

=0 for standard, for nonstandard. Default value = 0. 

G modal damping coefficients of structural, solid-friction damping, gj. 

(See eq. (A7).) (Damping force proportional to amplitude, but in 
phase with velocity). 

HIC three-dimensional array of input values of modal deflections from which 

the values of hj[(x,y) of equation (14) are interpolated or used. In 
HIC (JSPAN, KCHORD, IMODE) subscripts JSPAN, KCHORD, and 
IMODE indicate the index of the span station, the chord station, and 
the deflection mode, respectively, of the value of HIC. JSPAN 
increases inboard to outboard, KCHORD increases from leading 
edge toward trailing edge. 

ISPS input array of spanwise- integration factors, which the program multi- 

plies by s to produce the Ig. of matrix (6) of equation (16) for non- 
standard spanwise quadrature. Omit if CWSPI a 0. 

INTOR order of chordwise interpolation (1 if linear, 2 if parabolic) used in 

obtaining the deflections hj of equation (14) from the input data HIC . 

MASS a two-dimensional array of generalized mass elements m^j of equa- 

tions (A3) and (A 13). 

NALPHA number of a values in input (see eq. (A12)), maximum of 25. 

NCDST number of chord stations used in every chordwise quadrature, =10 or 

20 for standard quadrature (with CWCDI = 0); a maximum of 20 for 
nonstandard. 

NHP number of h (deflection) points input in HIC for each chord, maximum 

of 20 


CWPR2,CWPR3,CWPR5, 
CWPR6, CWPR8,CWPR8X J 
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NPNCHQ 

NRQU 

NSPST 

OMEGA 

OMEGB 

QIJ 

XCBAR 

XIC 

XTAB 

YSPS 

$END 


control word for punching matrix of (QU) and seven other parame- 

ters for subsequent use in NAMELIST input: =1, punch; do not 
punch. Default value = 0. 

control word: =1 when QIJ is included with IP2182 and its calculation 
is not needed; #1 when QU is to be calculated. Default value = 0. 
number of spanwise quadrature stations: 3, 5, 7, 9, or 11 for standard 
quadrature; a maximum of 11 for nonstandard 
array of input modal frequencies 
base or reference value ct>„ (see eq. (3)) 

nondimensional generalized-force-matrix elements Q^j of equa- 
tions (11), (18), and (A12) 

array of nondimensional local chordwise coordinates x^^. at which 
hi(Xc,y<jj of equation (14) of reference 2 are input and used in non- 
standard chordwise quadrature. (See eq. (B2).) Same for all 
chords. Omit for standard quadrature, 
array of integrating factors for nonstandard chordwise quadrature, 
which the program multiplies by the local chord length to produce 
the of matrix (2) (eq. (14)). Omit for standard integration, 

chordwise array of x (see eq. (B2)) at which the deflections HIC are 
input; the same for all span stations 
array of nondimensional span stations y/s for nonstandard spanwise 
quadrature. (See ISPS and NSPANST.) 


The '’standard integration" (for obtaining the generalized aerodynamic forces) 
referred to with certain of the input variables is that described in reference 2 and is 
summarized as follows: The chordwise integration (see eq. (13) of ref. 2) is done by 
Gaussian quadrature, the values of modal deflection h^ being interpolated (linearly or 
parabolically according to INTOR) at each span station from the deflection data HIC sup- 
plied at local chord locations XTAB. The span stations (at which the chordwise integra- 
tions are done) are equally spaced for spanwise integration by Simpson's rule for second- 
degree parabolas. That is, the wing semispan (from y = 0 to y = s) is divided into an 
even number of equal- sized parts, and this division results in an odd number of equally 
spaced span stations. 

For nonstandard spanwise integration the user assigns CWSPI 0 and inputs the 
arrays ISPS of spanwise integrating factors and YSPS of span stations. The values of 
ISPS must be those for a unit semispan; the program multiplies them by s to obtain 
matrix 6. The values of YSPS are y/s. The array of input deflections HIC must, of 
course, be at the user-selected span stations YSPS. 


46 



For nonstandard chordwise integration, the user assigns CWCDI ^ 0 and inputs 
the arrays XCBAR and XIC. Furthermore, the input deflections HIC are not interpolated, 
except for the leading-edge values, for use in the chordwise integration and must there- 
fore be supplied for the chordwise locations of the user-specified quadrature scheme. 

The values of XIC must be those for a unit chord length; the program multiplies them by 
the local chord length x^.g ■ 

Listing of Input Data for Sample Case 


»iP2iao 

KPRQG=0. 

«E^0 

HT -7 FLUTTER MOOEL , CPTSi 3 MCDESr 8 / 22/73 
SIP 2181 

XK=.4, NMUUE = 3, TANLE=1.216, TAKLTE= «3:>2>& f S = .8l» NXX = 48t 
X(1) = .03..2,.^,.6..8,1. , X{7) = .1915» .5£35, .5355, 

Xt 101 =. 365 , . 541 ,. 7 i 3 l 8851 , 1 . 0571 , XI 151 =. 3745 ,. 532 5 , . 8485 , 1 . 0066 , 

XI 191 = . 538 ,. o 821 , .8 2 o 1 , . 97 '> 2 , 1 . 1142 , Xt 24 1 = . 5575 , . 687 o , . 8176 , . 94 / 7 , 1 . 0778 , 

XI 291 = . 707 , . 8231 ,- 9392 , 1 . 0553 , 1 . 1713 , X ( 341 =. 740 s , . 6426,. 944 7, 1 - 0468 , 1 . 1488 , 

X ( 39 ) =. 8761 ,. 9642 , 1 . 0523 , 1 . 1404 , 1 . 2284 , X ( 44 1 =. 9235 ,. 9976 , 1 . 0717 , 1 . 1458 , 1 . 22 , 
ym = 6 * 0 . 0 , 3*. 0 81 , 5 *. 16 2 , 4 *. 243 , 5 *. 324 , 5 *. 4 C 5 , 5 *. 466 , 5 *. 5 o 7 , 5 *. 648 , 5 * . 729 , 


HMLTPLyil 
0 HX ( 1,1 
-. 627 , - 
-. 605 , - 
-. 573 , -. 
■“ • 66 3 ♦ 

OHXl 49 flJ 


=1.0fHMLTPL7{2)=2 

hlti, -.455 

. 60 if 

bbbt '~«6b7i 
7 34 * ” # 696 f 
d0b9 -.b41t 


503 f "'• 267 f ”* 430 f ■^• 67 bi 
-• 553 f '• 6 l 9 f ^♦ o 64 f -“♦* fr 24 
■~#587f “*6369 **629f •*•664 
-# 657 f -♦ 604 » ** 667 # -•699 
-• 739 t -• 733 , -, 772 , -• bl 2 
-• 873 , -• bl 6 , -• bb 7 ,- K 049 
1 . 7 bC , l . b 30 , 1 . 620 , 1 ^ 17 Q , l 4 l 70 , I - 3 i 0 , I 


-• t >18 

-•648 

-•7<^5 

-.881 


070 , 

• 627 , 

• 782 , 

1 . 910 , 

1 . 440 , 

. 696 , 

1 . 130 , 

.223 

290 » 

• 650 , 

. 290 , 

.552 , 

2 . 600 , 

1 . 530 , 

. 104 , 

.790 

130 , 

1 . 740 , 

2 . 610 , 

. 370 , 

. 220 , 

. 370 , 

1 . 940 , 

1.730 

990 , 

• 530 , 

• 290 , 

. 240 , 

. 620 , 

. 960 , 

. 420 , 

.220 

071 , 

• 120 , 

. 170 , 

. 420 , 

. 350 , 

. 270 , 

. 200 , 

.150 


0 HX ( 97 , . 33 U , . 440 , . 440 , 


»140j 


i038. 


. 130 , 


. 150i 
-.018* 
-.500. 
-2.170t- 
-2. 990t- 
H( l.lj 

.330, .260, . 

.003, -.120, 
-•52Q, -.410, 
l.olOff -.660,-1. 
2.160,-1.530,-1. 
= .21i>0, .l69C 

.058, .100, 

190, -.200, -.170, 
770, -.910, -.610, 
680,-2.760,-1.940, 
790,-1.670,-2.040, 
♦ .072'^, -.0045, 

.300, 
— • loO, 
-.6o0, 
-1.610, 
-2.170, 
-.0670, 

.140, 
-.530, 
-1.650, 
-<: • 3 50 t 
-1.720, 

— • 1 690 f 

,lo30. 

-.•3090, 

-.1690. 

•07b6, —.0040, 

-. 0630, -.1680, - 

•2710, 

• 05:>3 f 

-.0310 t 


-.1755, 

-.2730, -.040^, 

-.117^, -.1955, - 

• 2650, 

-.3760, 

-.0790, 


-.1500, 

-.2330, -.3145, 

-.4000, -.2180, - 

.2900, 

3o45 , 

-.‘,966, 


-.5250, 

-.2900, -.3o3':', 

-.434^, -.5080, - 

•5800. 

-.4500, 

-.5200, 


-.5900, 

— .ooOQf — -7350, 

-. 5550, -.6200, - 

.6850, 

7450 , 

-.8200, 


H ( 49,1) 

= -.9500, —.6400 

, -.2875, -.0115, 

.2130, 

• 4660, 

-.5850, 

-.1310, 

.3870, 

-.3360, -.2280, 

.0080, .3250, 

.4870, 

-.3550, 

-.2610, 


. 1920, 

.4300, — .26:?0, 

-.2450, -.0055, 

.3370, 

.4lo0f 

-.2880, 


-.2610, 

-.1620, .l4:>0f 

.3425, -.2575, - 

.2460, 

1090, 

.0770, 


.2710, 

-.3000, -.25/0, 

-.2340 , -.1920, - 

• 1 060 , 

-.34^0, 

— . 3 loO t 


-.3050, 

-.2975, -.2040, 

-.4880, -.4590, - 

.4350, 

-. 4l65 f 

-. 4065 , 


HI 97,11 

^ -.1130, -.0490 

, .0460, .1055, 

.IlbO, 

.1400, 

-.0075, 

•1320, 

.1645, 

.1075, .loi:), 

• 189^^, . 1980, 

.2130, 

.1910, 

.2275 , 


.2325, 

.2320, .3205, 

.2960, .2690, 

.2425, 

.2190 , 

.4200, 



334 I — • 4 l b 1 


400 , 1 . 370 , 


310 , * 230 , 


47 



-37o0t .29tj0* • 1775, •39c5, ,3060, .2-tOOf .1100, 

-.1240, .409t>, .29d0, .1950, -.0^85, -.2950, .2110, .03^0, 

-.2090, -.4*<-o5, -.oOO?, -.066^, -.3120# -.4740, -♦59:>0f -. 7250, 

SEND 

HT-7 FLUTTfck MODEL, 48 CRTS, 3 MOCES, 8/22/73 
$IP2i62 

INTQR = 2, NPNCHO=L, kii>PST=ll , NCOST^IO, 

H1C( 1 , 1 , 1)=— .218, — • 192,— .158, — . 110, — . 047 , .036, .1^7 , .ZJ2f« 385,. 525, .o75, 

H1C( I ,2,1 } = ”.l59, — e lo3, —.119, — .0557 , — .0026, .079, .179, .29, .417 , .558, . /06, 
HlC(i,3,il = -.ic>9,-.127,-.'^7 9,-.<'n73,,04,. 12,.2l8,.327,.^51,.591,.73^, 

HIC( lf4flJ=-.124f — •Q83f-.037,.031 ,.07 9, .159, .255, .363,. 48 o ,.623,.7b7, 

HiCl 1 t5,n=-.Q7 2,-.Q34, .??4 1 , .r' 62, . 1 1 7, .196, . 29, . 396 , . 5 21 , . 6 55 , .796, 
HICIit6flJ“‘“.027, .0 09,. 044 , *C93,. 155, .233, .326,. 434 , .55/,. 686, .825, 

HiC(l ,7 ,1 )-.D046, .044 , •'1329, *121, .196, .273, .365, .471, .592,. /lo, .856, 
H1C1I,8, 1J=.U31, .077, .123,. 176,. 2 39 315, .405, .50b, .627 ,.746,.8b9, 
HlC(l,9fli^.067,.ll6f.l66,.225,.2 85 ,.358, .446 , . 546 , . fao2 , . 7b3, .925, 

NICi 1* 10, li^.l22ff *l6^f . 218, .273, .331, .4CC,.486,.5b4, .698,.b^l ,.962 , 

H1C( 1 f 11 , li 1 o8,.230f .271 f.3 20 ,.376, .442, .525, .b21 , .735, .860, 1 .00, 
HIC(l,l,2i=1.0,.72lt.545f .4 55 , .383, . 3 5 , . 328 , . 32b, . 385 ,. 504, . o39 , 
HXLU,2,2>=.b22,. 5b 7 , . 4 29, . 355 , . 309 , . 28 8 , .284 , . 30 • .364 , .4 88, .621, 

HiL ( 1 f 3,21=.u 4, .4o*t, .336, .286 ,.2b5, .257, *258, .275,«345,.473f .605, 

HiCi 1,4,2J = .459, .355, .282,. 261 ,.2 59 ,« 261, .254, .257 , .328 ,.459, .591, 

HlC(l,5,2 )*.28b, . 2^9 , . 228,. 233, .245,. 254, .246, .245 , .31o,. 447 ,.576, 

HIC( 1,6, 2)^.1 3 8, .131,. 136,. 149,. 165,. 182, .202, .234, .3 09 ,.436,.558t 
NlCCl ,7 ,2i“.Di4i, — .0046,— .308, — .006, .006, .034, .109,. 21b, .305, .427,. 53d, 

H1C( l,8,2J — — .101, — .l47, — .178t — • 193, — . 18o,— . 145,—. 013, .192 , .3 02 , .419, .519, 
HILll ,9, 2i=-. 2 13,-. 2/b,-. 32 5, -.34 7, -.3 37,-. 283,-. 1^7, .155, .298, .412,. 506, 
HIGH, 10, 2 J=-. 3 36, -.3 8 7, -.423,-. 431,-. 4C4,-. 343 ,-.2i , .108 , .^:91 , .40o , . 5 , 
HIC(l,ll,2i=-.4 6o,-.‘+b2 ,-.487,-.468 ,-.4l8,-.3 51,-.271f.05b,.264,.4,.49 7f 
H1C(I,I,3)=.I22, •03 5f-.C53,-. 166,— . 335,-. 454, -.505 #-. 49 8,— . 3 75, -.075, .53 4, 
H[C(l,2,3)=«0b8,.C07/, — .08,— • 191,-. 32 8, — . 42 ,-.451 ,— .41,— . 292 ,.0 bo , .587, 

3,3J = .049,— .025, — .1*^8 , — .213, — .^2 1, — .385,— .397, — .341,— . 211 ,.198, .636, 
i 1 , 4, 3X — .UOlOf — .Oo*tf— . 136, — *227,— .311, .351#— .347 .29b, — .129, .312 , .68, 

HICX 1 ,5,3)-— .04Of — .103, -.161, — .235,— .298, — .321, — .30b,— * 2b, — .034, .403 , .719, 
HXCXlfb, 3)-— .084,— .132,— . 179,— .236,— • 284,— .298,— .^74,— . 195,. Ob, .4 7 4, .757 , 
HICi I f 7,3)=— .lC5f — .l4o , — • 189 , — .234, — . 2b9, — .278,— •24,-.0db, .2 09 , .53 5, .796, 
HXC(l,8f3i=-.li4,-.l5i,-.194,-.233,-,256,-.256,-.i9,.04br.337,.595,.843, 
HIC(1,9,3J=-. llb,-.l54f-.198,-.232f-.243,-.223,-.ll, .183, .446, .659, .894, 

HI C a, 10, 3) =-.127 ,-.165,-. 2 05, -.2 32, -.23 I,-. 17b, C.O, .295,. 53,. 727,. 94b, 
HIwll ♦ 11, 3) — 14, — . lb, — .213, — .232, — .2 19, — .126,. 124, .3b8,.o0, .798, 1.0, 

NHP = 1 l,xrAb(U = 0.0f.if.2,.3#.4,.5,.6,.7,.8,.9,1.0, 

MASS( l,ii=. 00035704, MA S S( 2 , 2 ) =. 00051 78 , MASS < 3 , 3i =. 00026352 , 

QMEGAil J=l62.5, 391., 725., OMFGB=725., G41)=3«0.0, 

ALP HA X U 6.0004, .0006, .0007, . 0C08 OOC 9 ,. 0010 , .0011 0012 0013 , .0014, .0015, 
. 03 lb,. C017,. 00 lb,. 0019,. 0020, .0022, .0024,. 002b,. 002b,. 0030 , NAlPhA=21, 

$ENO 


Discussion of Printed Output 

The only printing from D2180 is one line that states the value of KPROG. Then 
follows from D2181 a heading, the contents of the IDENT card, and an echo of the NAME- 
LIST IP2181* Then comes a heading identifying PART I, and the case data K, TANLTE, 
TANLA, S, MACH, DELNOM, NSYMMETRY, and number of control points NXY. Then 
follows, for each control point in sequence, 
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X 


X 

Y y 

DSL DIL plus D2L (fig. 3, ref. 1) 

DIL width of subregion I, if any 

D2L width of subregion n 

DSR width of subregion IV 

PTSDlL,PTSD2Lri one half the number of Gaussian quadrature stations in sub- 
PTSDSR j regions I, n, and IV, respectively 

the half-width 6/2 of subregion III (the control-point strip), and a row of the Ilnjn 
matrix (16 complex elements). 

After the Hnm n^a-trix for the last control point is a heading that identifies the 
beginning of PRT2. Then, if the input includes a multiplier for the downwash input, a 
notice is printed regarding that multiplication. The complex downwash elements follow 
for each downwash mode. Then for each mode comes the 16-element complex 3-nm 
matrix and the matrix of residuals from the least-squares solution. ^The sequence of 
subscripts nm in a^^ is 00, 02, 04, 06, 10, 12, 14, 16, 20, 22, 24, 26, 30, 32, 34, 36 
for spanwise symmetry. For spanwise antisymmetry, m is one unit higher: 01, 03, 

. . etc., . . . 37^. 

The output from D2182 begins with a heading, followed by the contents of the second 
IDENT card input, and by an echo of the NAMELIST IP2182 input. Then the input values 
follow for the case data MACH, K, TANLE, TANLTE, S, NMODS, and NSYM, the matrix 
of deflections hi (x,y), the matrix of generalized masses my, the modal frequencies coj, 
the modal structural damping coefficients g^, and the matrix of generalized aerodynamic 
force elements Qy. Finally, for each value of the air mass parameter a, there is 
listed a, the flutter determinant, the eigenvector columns the eigenvalues O, the 

eigenfrequencies w, the structural damping coefficients g, and the stiffness parameters 
bQCOg^V = kRe(n) . Where the eigenvalue D, has a negative real part, u), g, and 
bo^B/V are printed as indefinite quantities. 

Discussion of Punched Output 

The punched output is obtained when NPNCHQ = 1 and consists of the QIJ matrix 
and the quantities MACH, K, TANLE, TANLTE, S, NMODS, and NSYM in the format for 
subsequent NAMELIST/IP2182/ input when flutter solutions are desired for additional air- 
density parameters ALPHA. For the latter usage, KPROG = 1 and NRQU = 1. 
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OUTPUT FOR SAMPLE CASE 


HT- 
ECHO OF 
SIP2161 
XK 
XH 

TANLE 

TANLTE 

S 

NXV 

NSYM 

□EL 

X 


i 


□HX 


NASA - LANGLEY RLSEAKCH CENTER - HAHRION, Va. 

JEAN FOSTER FOR HERB CUNNINGHAM- 021 BOi lNCLUDtS02 Ibl AND 02182^" ’* “ 

UNSTEADY LIFTING-SURFACE THEORY BY THE SUPERSONIC KERNEL FUNCflON' HETHQD ANo A 
GALERKIN MODAL FLUTTER ANALYSIS UF ARRQWHEAU PtANFORMS WITH SUBSONIC LEADING 
EDGES AND SUPERSONIC TRAIL INu EDGES 

***REF^ 1 = AIAA JOUR.r N0V«1966t P«19ol~19oB 
♦♦♦REF« 2 = NASA TN D-A012 

3 * NASA TM X-29l3t THE DESCRIPTION OF THIS PROGRAM 

FLUTTER MUUELf 48 CPTS, 3 MOOES» 8/22/73 

NAMELIST 1P2L81 FOLLOWS 


0*4E+00t 
0 * L 2 E^ 0 l 9 
0.1216E+0L9 
0 . 352 S 3 E +009 
O.BlE+ 00 ^ 

4 d* 


I. 


0,4£-01, 


0.3E-01, 
0 . 19 l!»t +009 
0 * 7131 E +009 
0 « d^doE^OO V 
0.9 702E + 00» 
0.9477E+00* 
0 . 1 O!? 53 E> 0 i* 
0. I046db*'0i , 
0 . 11 V 04 E ^019 
0«ll9f>bE>0i , 


284 - 00 , 0 p 4 E 4-0 

0.5635E+0Ct 
0.8851E*C0, 
^• 10066 E* 01 , 
0.11142E + 0U 
0.10778E+01, 
0.11713E+01* 
0^1 1488E4-01 t 
0.12284E4-01 9 
□ • 122 E +01 • 


0, 0.6E«00, 

0.935l>£>00* 
0.lC371E*-01f 
0.538E+O0f 
0.5f>75E40U» 
O.7O7E+0Of 
C.7A05E+00* 
0-d761E+009 
0 . 9235 E 4-009 


0.8E4-0D» U.iE+01, 

O.Jb 9 E 4 -OOt 0 . 5 ^ 1 £ 4 - 00 » 
D.374i>E4-00f 0.i>32>E4-00, 

0.bd21E4-00t 0,d26LE*00, 

OM6t>76E^OOt O«dl76E>0Oi 
0.B231E4-00, 0.9392E + 00* 

0*d<^2oE4-00, 0.944 7 £4*00 f 

0.9d42E+00i 0-10523E4-01 

0.997dEt0Cf 0.10717E+01 


0 *0 « 0 .0 V 

"'-Of 0-0 

0. 162E+U0, 

0-162E4-00 

0*243 E ♦'00 t 

'^.243E4-'^0 

U* 324t-f00f 

0. 32464^00 

0.40:jE*00 f 

'^.405E4-^r 

0.4doE4-QQ^ 

0-486E4-CO 

O.So7£+CO* 

0*S67E^OC 

0 •tj4dt4'00* 

0*64dE4-CC 

0*729E*00» 

0*729E+C0 

-O.ld^E+OO# 

-'^.381E4-^r, 

'0* o 7dE4-G0« 

-0*334E4-0C. 

‘0.473E+00f 

-''*455E^rA, 

‘'0*424E4-uO» 

-0*605E4-00, 

-0.b87£+00t 

-O-GBSEi-OOi 

■0.3ddE4-00» 

-0.657E4-^C, 

■0* 66 7E4'CQ, 

“^*699E4-G0, 

■0* 70>t+00, 

-0.739E4.CCi 

0- 763E4-00* 

-C*a08E4-C0, 

■O.diot+OO, 

-0*887E4-cCi 

0*l62E4-0i, 

0-117E4.01, 


O.Of O.Of 

O-dlE-Oi, 

0*162E+O0f 

0-162E4-U0, 

')*243£4.00» 

0.243E*00t 

0.324E+O0t 

0*324E4'OOt 

0.405E4-00, 

0.40f>E*O0t 

0*486E4-00f 

0*4b6£4’00 9 

0.567E+00, 

0.367E4-00, 

0-64eE+00» 

0.64dE4'00 V 

0.729E+00, 

0.729E+009 

-O.S03E+00, 

-0.2o7E4-00f 

“0*41dE4-00, 

*0. 62764-009 

-O.&53E+C09 

-0.6l9E*00f 

-0.601E4-00, 

-0-:io9E4'009 

-0.629E+00, 

-0.6i>4E+D0f 

-Oq640E4'OOi 

-0.6S7t4-009 

-0.663E+00* 

-0.734E*00, 

-0.733E4-00* 

-0-7726+00. 

-0*d41E*00. 

-O.ddlE^OOf 

-0.1049E4-019 

0. i7dE+0i, 

3.117E+01, 

a.i3iE+01f 


O.dlE-01, G.dlE-Oli 
U*162E4-00, 
0*324E4>00, 
0 . 403 E 4 -O 0 f 
0.4d6E+00» 
0.367E4>00i 
0 . 6 ^d£ 4 - 00 i 
0.729E+00» 


-O.^DEf-OO, 
'“0 . ^76 E+OO r 
“ 0 • ot>4£4-00, 
-0.&l6E4-00» 
-0.373E+00, 
- 0 .o 04 E 4 ' 00 f 
- 0 . 696 E^ 00 1 
B12E+00, 
-D.d 73 E 4 ’O 0 f 
0.1diE4-01, 
U.14E4-01r 
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0*7£j4tE+00t □•L9lE^0if 
0*229E+01# 

0<*2oL*-01» 0.I53E + 01* O.lO^E + 00* 

0*l74fc<-01» 0.26it + 01, 0.3/t + 00» 0.22E+00» 

'D.LTSE^plt O.^Vt^-OQ, O*53E*00t 0-29E4-00r 

0.96E+00, 0.42fc+00i 0-2^t+00f 0.71t-0lt 

^.42E<-00, 0.3:?fc*-00, 0*27E*00, 0.2E*00t 

0.44E + CC* O.^^t^OOi 0*l4E^00, 0*3dfc''0lf 

0.23E + 00, 0.13E*-0Uf 0-33£^00, O*2bE*0Cf 

O.lE^OOt 0-3E+00, 0.14fc+00t -O^ltiE-01# 

0-3£“02f '0*12E4-??, -0,W6*-0C« -0.2E*>00» -0.1/t^00» -0.16E + 00, 
-D.53c*OOt -0-5E+0C* -0*32Ei-CC* -0-4lt+00, -0*77E+00f -0.91E + 00i 
-O.oLE + OC. -0,6aE+OOt -0*l65E4-01, -0.2I7t+01» -0-i6iE + 0i, -0.o8E+00t 
“O.lodh+01, -0.276e+0lt -P.194£+01» -O.lbLE^Ol# -0,235E+01t 
-0*299ti*Cl, -0-218Ef01, -0^153£ + 01, 


0.137t^01 
0, l44Ef QL 
0 E+OG * 
0. 79c+00t 
0*37t^90« 
0.24t*-0G* 
0- IZE^OO, 
D.L3E+0Gt 
0, 13E + 0C, 


* 0.107E+^1, 0*827E+00i 

► 0.698E4-00, 0*LL3E+0l« 

0*29E«-T"'t 3.352E+C0» 

-U.13£-»'00t 

0-194E4-O1 , 

\Jm 62E ♦*00 j 

0-1 TE^OO, 

0*33E*00, 

0*31E+O0f 
3.58E-01. 


~0*i79£4-Ul, -0.1d7£^-01» 


0 .2Q+E+01 ♦ 

-0 

.217E+01, 

-0 

.172E+01, 

If 

1 t 

1 f 

1 f 

If 

If 

If 

1 1 

I f 

1 f 

1 f 

t 1, 

1 « 

1 , 

It 

1, 

I, 

I, 

It 

I » 

I f 

If 

It 

1, 

It 

I f 

if 

1. 

If 

It 

It 

If 

It 

* It 

It 

It 

it 

I. 

It 

I. 

It 

It 

I, 

It 

If 

if 

If 

If 

It 

If 

If 

If 

If 

if 

If 

t I. 

It 

1, 

i f 

It 

I. 

It 

It 

I , 

I t 

It 

If 

1 f 

It 

1 1 

it 

if 

It 

If 

if 

If 

If 

, It 

It 

if 

i f 

If 

It 

I, 

It 

it 

it 

If 

If 

if 

If 

If 

If 

It 

If 

if 

if 

If 

If 

t it 

It 

It 

If 

I, 

It 

It 

It 

I t 

It 

1 f 

It 

1 . 

If 

If 

If 

If 

If 

if 

If 

If 

If 

t I . 

i . 

i . 

I * 

It 

It 

1 f 

If 

I t 

It 

It 

If 

It 

If 

if 

1. 

If 







H = 0-2Ii^£^0a, •lb9Ei'C'^# Op72E-Clf -0^43£-02f -0.b7fc-0If -0.lB9E+00t 

Q.lbit^DO, -0.9E-02t -0pl69£+00t 0.7d5L-01, -0.4t-02# -D-aSE-Ol, 

“Opibdb+OOt -0.271E+00» 0p3&5E-01f -0*31£-01» -0 • 1 733E+00t 

-Op2 73£+UOi -Op^E-Olt -OpUTE^-OO, -0.1 9ii3£*00, -0^2d3£+00» 

-0*37t>£+00. -'>.79E-01f -0.156 + QC» -0.2j3t«*00# -0* Bl^bE^-OO, -0.4£^*00# 
-0.21d£ + 00t -0.29E^00f -0.3645E + 00t -0.4466t ♦OOt -0.:>23H+00f 
-0.29£*00t -0 .363E>'^^< -0.434E+C0, -O.^Odt+OOt -0.3dt + 00» -0.45£ + 00» 
-0.32t + 0o* -0.59E + O0p -0.66E<'CC# *0.73bt + 00f -0*33&£+00f -0.62E+00» 
-0.oa3£+C0» -0.745E+00, -0.82E+00, -0-9p£*00, -0*d4£*00, -0. 2«75£^*00f 
-0.113E-01* 0.213£ + C0p 0.466E+O0t -U. :>bPt ♦* 00 t -0-131E+00f 

0 .3d7£+CG» -''.33BE + ^0i -0.22aE+00f 0^d£-02t 0.325£ + 00i 0.487E + 00* 

-0.333t<-00, -0.26ie+00t 0.192E+00, 0*4^£ + 00f -0.2b3t+00f 

-0.243E^00* -0. 55E-02f 0.33 7E*0C» 0*4lb£i-00p -0.2aa£+00t 

-0.2oie+00t 182E<*00, 0.L43E + 00t 0*->t4:3£ +001 -0*23 7i>b + 00t 

-U.246E+00, -0-109E+00, 0.77E-0Lt 0.271t+00f -0.3£4*00, -0.237£^00t 

-0.234£t00i -0.1926+00, -0.1086+00, -0-3^3fc+00, -0.316E+00, 

-0.3Q3E + 00, -0.2975E + CC, -0-264E + 00, -0.4od£ + 00t -0.4!>9£+00, 
-0.433E+C0, -0.4185E+0C, -0.4055E+00, -0.113£+00, -0.49E-01, 

0.4bt-01» 0.1055E+00, 0.118E+00, 0.l4c+00, -0.73E-02, 0.132£*00, 

0.1645E+00, 0.1075E+CO, 0.161E+00, 0.id9£+00, 0.l9dE+00, 

0.213b + 00f n.l9lE+00, 0.2275E+00, 0.2323t«*00, 0.232E + 00, 

O.32O3E+0Ot 0. 2986+00, 0.269E+00, O.2423E+D0, 0.219E+0U, 

0. + 2t+00f 0.376E+0^t 0.29bt+00, 0-2!>ot+00, 0.1773E+00, 

0.39o3£+00, O. 3 O 6 E+OC 1 0.24E+00, 0 -ilE+OO, -U.124E+00, 

0.4093E+00, 0.298E+OC, 0.195E+00, -0.4d3t-0l, -0.293E+00, 

O.211E + G0, 0.34E-01, -0.209E-+00, -0.44o5E+00, -0.t»E + 00* -0.66E-0I, 



-D 

.3l2t+00. 

-0. 

474E+00, 

-0. 

3936+OOf 

-D 

.723£+00f 
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HMLTPL7 = O.lE+01, 0.2E+01, 


SEND 



PART I Qf 02181 GENERATE IINH HATRIX 


K TANLTE 

4.00000000E-01 3*i?2550000E-01 

NO, Of- C«PTS.(=^NXirf 
48 


TANLA 

i.2it>oooooe*oo 


S HACM 

6.IOOOOOOOE-OI l«2000QQaOE*^00 


DEL MOM 

4,OOOOOOOOE-02 


3.00000000E-02 


PTSDIL 

0 


PTSD2L 

5 


HALF-hlOTh 01= CCNTRQL POTNT STRIP 


OSL OIL 

L.4366a60BE-02 0« 

PTSDSR 

5 

L.59o3l7&6k-03 


ROW 1 Of- CQEPPiClbNT MATRIX ! !KM 


RE 

4.7400‘*833E*00 
-l,0200l030e-li 
-6. 3720l494£-0b 
-7* 3990684BE-D4 
2. 7321952 7E+-00 
-7.600784l5fc-l2 


IMAG 

1-036^918BE-01 

1.91S04522E-13 

l 4 r' 20307596-09 

6417071446 E -08 

643213S190E-02 

1-43001278E-13 


Rt 

-8.74>460l0ofc-04 

4.371992718^00 

- 9473 V 907 T 3 fc -12 

^5.93tt72373E-0B 

-64l43tt0432£-04 


02 L OSR 

1*4366B608E-02 1« 43668608E-02 


IHAG 

7.2946 5038E-0& 
9.62B53525E-02 
1.83163797E-13 
9422390544 E -10 
5413573241 E -06 


RE 

-6.9012731 3E-08 
-B.29653527E-04 
3«67640126E4 00 
-B.85359850E-12 
-5.04648391E-08 


X 

Y 

OSL 


2.000000QDE-01 

0. 

9. 5779071 7E-02 

0. 

PTSDIL 

PTSD2L 

PTSDSR 


0 

5 

5 



D2L OSR 

9-577907176-02 9 . 5779071 7E-02 


HALF-WIUTH Of control PQENT STRIP 


l.0642ll9lt-02 


RU- 2 OF COEfFlCltJ 
RE 

4. 7 b 74644 OE 4-00 
-8.861488118-07 
-9.43559432E-05 
-5.642517 o 7E-03 
-2.48155710E4UO 
l*0642:>Il3E-a 7 


MATRIX TINP 
IMAG 

6.B6743974E-01 

1.115600346-07 

9.84505588E-06 

3.64^172306-04 

-3.313000796-01 

-1.26297M46E-08 


KE 

-3-d7969774£-02 
2.20694221E400 
-6.36799452E-07 
-3.2739 7969fc-Cl6 
1.09093482E-02 


IMAG 

2.147U711E-03 
3.675611786-01 
8.043431196-08 
3- 54256769E-06 
-5.32759607E-04 


RE 

-U351922776-Q4 
-2.53123758E-02 
-1. 0052042 5E400 
-2.55U3560E-07 
2.30910575E-05 


X 

V 

DSL 


4.000000006-01 

0. 

1.928423826-01 

Q. 

PTSOU 

PTS02L 

PTSDSR 


0 

5 

a 



HALF-WIDTH Of CONTROL POINT STRIP 2.00000000E-02 


D2L OSR 

1. 928423826-01 1. 92842 382E-01 


ROW 3 OF COEFFlCIfcj'jT 
RE 

4-84447078E40C 
-5.491766826-05 
-9 . 10714940E-04 
4. 75762074L-02 
4-063l06o2t-0i 
1 .S48 86565t- (J D 


MATRIX IINM 
IMAG 

1.349B94l7Et00 
1.4M10C733E-05 
1.972400726-04 
-4.38747998E-03 
-2.93899561E-01 
-4. 878646096-06 


Kt 

-1.5332ol7ie-01 

-5.>9t>4lVl2E-Dl 

-2.041972586-05 

3.O0364929E-04 

5.^33o0544E-02 


IMAG 

1.68377327E-02 
1.09921166E-01 
6. 870446636—06 
-6.33655145E-05 
-6.20152507E-03 


R 6 

-2.10897208E-03 

-4.96076225E-02 

-3.19926747E400 

5.39‘130707£-06 

7.78o78>64E-04 


X 

6. OOOOOOOOE-01 

Y 

0- 

OSL 

2.99263572E-QI 

0. 

OIL 

PTSOIL 

0 

PTSD21 

a 

PTSDSR 

8 




02L OSR 

2.99263572E-01 2. 99263 572E-01 


half-width uf CuMKUL POINT STRIP 


ROW 4 OF COEffICiENT 
RE 

4.95201113t400 
-5.928304? It-OA 
-2.018389551-03 
1.40584350E-01 
4.59089168E+0G 
1.4iab745ut:-04 


MATRIX IINM 
IMAG 

1.973152786400 
2.31969971E-04 
7.55374178E-04 
-2.07047544E-02 
9.694403981-0 1 
-5.695784806-05 


2.000O0C00E-02 


kE 

-3. 3 644 809 IE -01 
-3. 6080508 6E *00 
-1.073 741186 - 04 
3.31815905E-03 
-6.3735o790E-03 


IMAG 

5. 503982 83E-02 
-7.25121470E-01 
7.03O71205E-05 
-9.39446653E-04 
-3.738708236-03 


RE 

-1.024l90s>5E-02 

-3.51690474E-C3 

-1.1235237^6400 

1.51661883E-04 

1.7V795043E-03 


X 

8. OOOCOOOOE-01 

Y 

0. 

DSL 

4.04400525E-01 

0. 

PTSOIL 

0 

PTSD2L 

6 

PTSDSR 

B 



OZL DSR 

4. 0440 05 256-01 4. 04400525E-01 


HALF-WIDTH Of CLNTROL POINT STRIP 

RUW 5 OF COEFfiCIENT MATRIX t INM 
RE IMAG 

5-O6B94820L4OO 2.551674936400 

-3.088856O4E-03 1.664215896-03 


2. 12842362E-02 


KE 

-5.8p907b26b-0l 

-6-95580079E*00 


IMAG KE 

1.25079596E-01 -3 , 0 5792844E-02 

-2.07654832E400 I .707433036-01 


NSVNHETRT 

I 


IHAG 

1.07124020E-09 
6.90942694E-06 
8.236 72533E-02 
1. 6652495 lE-13 
7.84296966E-10 


IMAG 

1. 4016688 5E-05 
1-4247762 5E-03 
-6.0533U13E-02 
3.27819066E-0B 
-2.23929356E-06 


I HAG 

4.398732346-04 
6.089434676-03 
'7. 467010826-01 
1.14210452E-06 
1.60235851E-04 


IMAG 

3.235555636-03 

4.929264516-03 

■8-39407132E-01 

•5.36833647E-05 

■6.494201116-04 


IHAG 

1.30571603E-02 
■1. 98450 203E-02 


52 



4.06-^369766-04 

6.59964653E-04 

-3.5S137244E-04 

2.53319436E-04 

5- 99001315£*-00 

7.350658O5E-01 

l.2OS250l6t-0l 

-3.019821416-02 

8. 61 >440636-03 

-3.43469930E-03 

8.02524663E-04 

-3.939556196-04 

-i,24lb6023t+00 

1.0457C932E«-00 

-1. 631751266-01 

2.23833921E-02 

-1.56370745E-03 

3.12195010E-05 

l,T7l940lst-3^ 

-1.36466489E-04 





X 

Y 

DSL DIL 

02L 

DSR 


1 ,00000Q00E*Q0 o. 

S.CS50Q6S6E-01 0- 

5.055006566- 

-01 5.0550D656E-01 


PTSDIL 

PTS02L 1 

31SUSR 




0 

8 

e 




MALF-MiOTH UF CuNTRQl POINT STRIP 

2.660529?7E-02 




6 OF CUEFFlLltNT 

MATRIX IINK 





RE 

IHAG 

RE 

IMAG 

RE 

IHAG 

5. 16766783E+0U 

3.0961808B6*00 

-8.9Q775585E-01 

2.32365561E-01 

-6.953880446-02 

3.77732259E-02 

-1.068783b>t-02 

7.515192216-03 

- 1.0 5856 8666 *-01 

-3.b9439l43£^00 

5.208342936-01 

-9.06709 4276-02 

1 .44747555E-02 

-4.141589036-03 

3.O0876713E-04 

1.99284872E-04 

1.888137166+01 

4.926659166*^00 

-2.646299336-01 

2.303480766-02 

6.080117456-03 

-4.42160746E-03 

1.72663320E-03 

-1.165966896-03 

-3.0800361b ti'Oi 

-5. 311255366 *00 

1.26781383E-01 

8.03020836E-03 

-5.63951917E-03 

2.62042586E-03 

-4.331408806-04 

1.315542916-04 






Output for downwash points 7 to 47 not shown 



X Y 

1 .24i000a00t-»00 


OSL OIL 02L DSR 

1- IC211439E*00 3.i9A»^7od:^fc-Ol 7- 836!>4»703fc-01 i*!>9?28642£-01 


PTSOIL 

8 


P 7 SDPL 

a 


PTSDStt 

i* 


MAlF-WIOTH of CuNTKUL point strip 


l*77^7*4»92£-02 


Ra*» UF COtFFlCltNT MATRIX 1 TNM 


Rfc 

3.078733l4t-0l 
-1.714<*0132E*CC 
6 .18j966V2L-01 
2.3072a7«jok-Cl 
-1.3to77009ob-Cl 


1HAG 
3.62648213EOO 
4.339574C0£-0i 
-3*53337521E-01 
4,6 3570668E-01 
-4* 30066946E-01 
-2.41L81865E-02 


4.91739907t^00 

-5.:>'4970799E-0i 

>L.a^tt 69 i 9 iE *00 

6 * 49 lfe 23 lab- 0 l 

-l*3394o!>36b-01 


iMAG 

1*3342 2 6*i76^00 
-2*V6973966E*‘00 
-1.69023457E-01 
l-7208980BE-’01 
-1.720 S8246E-0I 


KE 

2.44403309E4-00 
-2.29 79li>31t+00 
-e*38903682E-01 
4.660D7O49E-01 
-1.7616291OE-01 


IHAG 

7.I9471056E-01 
-8.299S8B50E-01 
Ue 4 B 927 fi 8 E «00 
T*2449eil7E-02 
-6. 0341322 9t-02 


PART n OF U2lt)l *** SOLVE SIMULTANEOUS EwUATIClNS FOR ANM 


♦ •***NUnLt- IRPUT OQWNVJASH QUANTITIES HAVE bfctN KULTIPLitU AS FULLOhS - BY 1.000 FOR THE SLOPES OHXt ANO 
BY 2.0UU Fu^ THE DEFLECTlLiNS H 


UJiiNwASH m/V* MUUt I lit NXY REAL PARTS FlhSTf liM LuNTRuL FUlNl SEUUENCE 


-1 .8S03C0COL-C 1 
-3.340000oat-01 
-3.53000000L-CL 

-s.t> 90 oooaot-ci 
-S./30000JOL-U1 
-b.d7000000E-C 1 
- 7.390000UUE-U1 
-8- 4 1.0000 30 t-C I 
1.72000000r-Ol 
I .30*»COOOjc-C1 
-1.34400001C-01 
- 3 , 2 O 00 iJ 00 aL- 0 i: 

-1.2COOOOOOL-C1 
-2 .91bQQOOOt-Oi 
-^.Co400000fc-0l 
-S . 8dOOOOJOd-..31 


-4.180COOOOE-01 
-6.l90COOOOe-'}l 
-5.180'rOCOCE-Ol 
-5.88C0CC00E-D1 
-6.9 9?'^'>'?40e-0 1 
-7.33000000E-01 
- 8 . 8 X^'?r 0 ' 3 'i £-01 
l.?52C?C00E-01 
-7.2OOCCOO0E-O2 
-2. 166Crrt0<'£~01 
-9.3 600000ab-C}2 
- 1 . 864 '^OO'^OE-Ol 
-3.5728C000E-01 
“4,64CO0000e-01 


-3. 0300000 OE-01 
-6.270000UOE-01 
-6.o^000000E-01 
-5.d7000000t-Ci 
-6-D/tiOUOOj£-01 
-c».o30C0000E-Ci 
-7*7200000gE-Ol 
-8.73O00DO0fc-Ol 
5.7b000000t-02- 
-1.332'JCOOOt-Oi 
<»-‘»*iOOOOOOb-02 
-1.3O900300E-01 
-2-ilbOOOOO£-01 
-9.ZOOOCOOOE-CI 
-B.oOOuOOOOE-Ol 
-<i* 9 QOOCfOOOE “01 


-2-b70000COE-Cl 
-4 *.VuOOOOOOE- 01 
-4.24000000E-01 
-toH.3bOOQOOOE-Ol 
-o. 460000006-01 
-7.3<»DaOOOO£-OL 
-b.l2u000006-01 
-8.16000000E-01 
-3,bOOOOOOOE-03 
b.2bOOOOOOE-02 
-2.‘»d00O0C0£-02 
-2. 2 6 0000 00£ -01 
-3. 200000 OOE-01 
-2.32000000E-01 
-4-ltoUOOOOOE-Ol 
-S.^oOOOOOOE-Ol 


-«t.30000000E-ai 
-^.730000006-01 
-6.03000000E-01 
-6.Z9000D0OE-01 
-o. 370000006-01 
-6.96000300E-01 
-7.830000006-01 
-8.67000000E-01 
-3.J<»000000£-02 
-3.200000006-03 
-1.^04000006-01 
-3-JOoO 00006-01 
- 1.74<»OOOOOE-Ol 
-2.90^003006-01 
-4.72000000E-01 
-3.9O000000E-01 


-6.780000006-01 

-4.S5CCOOOCE-01 

-6.OIO0Q000E-OI 

-6.3^0000006-01 

-6.0400000CE-C1 

-7.0SOOOOOOE-01 

-a. 080000006-01 

-1.049COOOOE4-00 
-1.31200000E-01 
-6.6400000C6-C2 
-2.18400000E-C1 
-6.320000006-02 
-2.32 0000G06-01 
-3. ^720000 06-ni 
-5.2d00Q00CE-01 
-6.56000000t-0l 


53 



W/V, Mj,>t t 


L , /b0000:^i^t*00 

I * ^cacoooot+oc 

1 •44000000t<^00 

2 .'JODODOOCt-Cl 
-I -3O0030C0L-C1 

I , 9<*00uQU Jt*OU 

b^zoQooooje-ci 

I . /OOOOOOOe-Oi 
-7.6Ci200000t-01 
-^.o600'3000h<‘i 

Z .bOOOOUOiit-Cl 
-2, UOOOOOCrc-oX 
-2 .OdbOOOOOt-Ol 
-d . 72u00000t-C2 
-i .?3oOOCDdc-Cl 
-2.27200i3jOt-Ui 


UJn'irtftiiH aH* noubi 3i 


3-3C000300L--X 
3 -lOOOL^OOOc-Ci 
i»,dC0QC03Cc-02 
-1 .2C000000L-01 
-3. COOOJOOOE-01 
-t>.BC300000L-Cl 
-2 . fbOOOOQOt*OQ 
-l-!^3000033t*0':' 
-•9. 0^003 3 03 t-C^ 
-6 -OOOUDO jJL~Oj 
i. bb^OOQ jOt-Cl 

2 .i»6^00000c-Cl 

3 .OOoOs^OOyc-Ul 
1 .9200C073L-C 1 

-i.aoULH>OOUL-ui; 

aO300QO0..-:i 


rtblbHTlN;, bACTOiiS fllsd* 


— i • cjO Oooo 3 3 L - 0 2 
lx I9^7i>ti3 1 r ♦Ci 7 
I . 3<* /Jf7lo7c-C-- 
-9«993o9^o2t-02 
-1 .320a900ie-03 
-3- 932929v7c-02 


LtA^T SwUAKL nL^iULAL^ 


“1 *403oci L9 7t-" 2 
-d.S<>251023l- 02 
2.235Si2jjQb-i:2 
2-33997^30^0^ 
-I .ooil>l)2o/t“G£ 
L * v3o5k^3 1 3t-C2 
1 -2 lo3V‘»d3t'G2 
-I - pdb 7t.i»oaL-0Z 
7-93i>2o353t-0* 
2- ll2959dHL“u3 
l>2oo39 92ob''02 
-1 - a 7V207‘iaL-Lf3 
'*t> « t>2 923 I J'Jb - 33 
-3-399 720O n-ip 
3. l3bt>t*0 3yc-09 
1 .9l2377o:>c«03 


NXr RbflL PARTS fiRST, [K i^UNThOL PuIhT ituUcNCt 


l-0?OCOCOOEtOO 
1 .3 7''''C000£*30 
6.9aorocno£-oi 
5.52CCC000E-01 
1 ♦ 7^^?00'!*OE*00 
! .730C0000E<^a0 
9.60''^'^000£-01 
A.200r0000t-01 
-5.!2aOCOOOE-Ol 

3,696CCC0CE'-6l 
-1 .SbCOOOOOt-Ol 
-1 -4b6rrcooE-Oi 
6-I6000000E-02 
-0,64'>^"''’if^E-O2 
-^.9C4C^000E-01 


l*s>20C0000t*00 
L.07000000t*00 
t.i3000000t>OG 
2 -bOUUUUUUk^LU) 

2- O1000000&+00 

9*^0000000t^01 

4.2i}OOOOOa£-Cl 

3- >oooaoook-oi 

-2.300G0k)UG£-0i 
3.09pCOOOCfc-Cl 
-2.d4000000t-OL 
-4.40a00000t-03 
UltxOOOOOOt-Cl 
2-io60ooaut-ui 
-2-7i*030000fc-CI 
~3*p 72 OC OOOc~C 1 


l-17000000E*00 

tt-27000000E'01 

2.2^0000COt-01 

1- L>iOOOOOOE+00 
3- /OOOOOOOt-01 
!».30000000c-01 
2.200000C0e-01 
2. 700000 OOfc-01 

-9. 200000 00£-03 
-2. 704000C0E-C1 
-2.0od00000L-01 

2- o9o00000e-0L 
2. 74OO0OCOfc-Cl 

-2--*OOOOOOOt-Ol 

-2-S2dOOOGOE-Ol 

-3-4dOOOOOOe-01 


1- lVOOOOOOt^OO 
7*d20GO000£-01 

2- 2S300000E+00 

1- O400O000E-01 

2- 20000000E-P1 

2- 90000000E-01 
7.10000000E-0S 
k-OOOOODOOfc-01 

1-704000006-01 

-l.d240000CE-0l 
I -bioOCOOOE-Ol 

3- 344000006-01 
--x-aoDC3D00fc-Cl 

- 2 *o;>bcoooofc-oi 

-2.4t000000E-01 

-3-34OC00O0E-01 


NXr REAL PARTS FIPST, ti> COMKOL PulM SfcwUtKtfc 


2-^ ro'^cccot-oi 

1 , _r^ 1 

-3.2000C000E-01 
-l-65?''nCOO£*0O 
-1 *940O^C00c*-00 
-l-790?rf'?0E*00 
“"'-92^''^^''''E-C'2 
1.056PCC00fc-91 
’ . 7 *^4 70E-7 1 

2-3640 '‘n A r^t-n ^ 

2.3a4^c000t-01 

-2-'‘'6OCr0C0£-01 
-5 .2 8'''^<^077E-Q2 


4.4000000Qt-Ol 

l.:>OOOOOOGt-Ql 

3-OOOOOOOOt-Ol 

- 2 -OOOOoO vOt-Ol 
“4-1 OOtjoOOuL-Ol 
-2-17JOOOOOt+OG 
-l.dlOOOOOOfcFOD 

— 1-o/OOLrOuOfc+OO 

3-ot>000C00t-02 
1.3luOLOOOi:-Cl 

1- >2tfOOOQOb-OL 

2- i>200000b-Cl 
2.04600000t-Ul 

-9.92000000t-02 
l.booOOu>OGL-Gl 
-2 •'•4i>u000C£-0 I 


1-400000COE-01 

3-300000006-01 
l.HUUOOOOOt-Ol 
-1-700900006-01 
-7.700000006-01 
-1-61000000E+00 
-2-3>OOOOOOb+00 
-2.O4OOOOOOEFO0 
O.44000000E-02 
d-oOOOOOOOt-02 
i-b2000000E-01 
1-94000000E-01 
1-42000000E-C1 

3-27o000CX-0l 
2. 720000 OOh-02 
-3-79200000E-01 


5-dOOOOOOC£-02 
2.b0000000£-0l 
-l.b0000000E-C2 
-i.bOOOODDOE-01 
-V-IOOOOOOOE-Cl 
-b-dOOCOGCCE-Ol 
-2.990C0DCCEF00 
-2-1700 ooooe <-00 
9.440000006- 02 
i.2oa00000t-01 
l.aoOOOOOOE- 01 
i.7:.20C090E-01 

3-17200000fc-0l 

2-3<J4000C0E-01 
-U67200000E-01 
-4. 7bOOOOOOE-Ol 


MCOF( n, 16 real parts mksi 


-2-913B19476-01 
-r-ir9?2995E-01 
-2-3 77!4r97t-3 I 
1.1 0799670E-01 
“3.92313979E-03 
5-‘97?Ob9lE-02 


6- ol2 >4ol oc-01 
1 -3b479baoc-0t 
a- >03c d7->* at — 01 

7- l407o992c-03 
4# »4di:3<i2t-C2 


-9-a0O36329fc-0l 
-4.5>704lllie-01 
- 4 . 900274 74E-C1 
4.910737336-02 
o.lal6l900E-02 


l-7obl4927E-02 
1.4223374bE*0C 
3 ♦3e0boob36-02 
-2-d402l909E-03 
-9.0347190SE-C4 


MODE! IJ, Nxy K6AL PARTS WRSF, IK CUnTRUL P□I^T SEwUEAfCE 


5.14616231E-33 
- 1 - 3629r6 7'>b-02 
-7- 60Cr6S79£-02 
-6.?33'>^51b6-93 
'•1,?0496169£'C2 
1. 9367^61 7E-G3 
-1. 566972046-02 
-6. 1 B429765E-03 
-6.225420816-03 
2. 793300216-03 
-2.20?a4''25E-03 
-1-04731 7daE-0 3 
-3.021149d3E-03 
-7-34514273E-C4 
-1.205 67609E-03 
1 -37'T954616-9B 


d>-a7t»a /uG^b— C 4 ; 
7.4iJ<»36l49b-02 
I.^o23o01 2b-01 
o.>l>b24o9t-G J 
5.073iix:4 7 Oc— 03 

- X -a424u24 1 fc-03 
3*2lb0i424t-02 
L *43.^9 / dbOE-02 

- 7-io3x.77x.ok— 0^ 
7-ob04u>aok-u3 

&-iab7907St-C3 

-4-349o 73 7kc-04 

- 3 « 094 oJdulk— 03 
2* >2 j>C0 9:it-04 
a*UjiioUl oE-04 

-2.00b97232L-05 


-1.5ibi447ae-0l 
2.560293296-02 
-4.o32o1613E-02 
7. 6S05Q9C0E-04 
-9.4d553550E-0? 
-2.5/230413E-C2 

7-lod22647£-02 
-5.2 74o0862E-02 
-7.7d:»2T79lE-03 
-4.223ad3ClE-03 
1.037711046-02 
3.091724106-03 
-2.5i?903aUt-03 
-l.d73bb8996-03 
2.402994l5t-03 
-3. 290099136-04 


1 . vO^ojobbE- C2 

k. dlloa3b4t-02 
4.72127654E-C2 
3.222o3b4b£-02 

-3 •4 3>a006dc-'02 
-1.09o2l469E-0Z 

l. 7d052774£-02 
-3.o792572dE-02 
-l.2a3579t>7E-02 
-2.04lop4l 96- 03 

7.045211016-03 

-6.39Gil310E-03 

2.0oi90dibE-03 

2.31225977E-C3 

3.OO313035E-03 

-3.02201741E-03 


l-3iCC00CCt*CC 
1.91 OOOOOQE+00 
6.500000GCE-C1 
7.90CC0000fc-0l 
3.70000Q0CE-O1 
2.40000000E-01 
1.2000000CE-01 
1.50000000E-G1 

3- 726CQ0C0b-01 
6-4O00D0C0E-03 

3 -44 0000 CCk- 01 
-2.3040000Qb-01 
-1.9bOOOOOOE-CL 
-Ufl72COOCCE-OL 
-2.3800000CE-01 
-3-2440000CE-DI 


1. 30000 OCOE-01 
8.000000006-02 
3.00000000t-05 
-5.3000000CE-01 
-6. lODOOOOCE-Cl 
-l-6bOOOOOCE40C 

-z.iaooocoofetoo 

-1.720000006+00 
l.UOOQOCOE-01 
1.51200CCCE-01 
1.8i>bOOOOOE-Ql 
3.3600GOOCE-C I 
2.44bC00CCfc-0l 
1 -56Q000C0E-01 
-3.572C000CE-O1 
-5.8OC0000CE-O1 


-3.d40il078E-Ql 
-1.493946056+00 
-4.7db44660fc-C2 
-1. 26622790E-O1 
4.27l5l4bab-C4 


3.9ttl22db?fc-03 
-I.U233311L-02 
-L.bd96473 Lb-02 
-3.474120C76-02 
-4.38ll6075b-C2 
2.38ai«L90e-O3 
-3.097957396-03 
4.202812o2E-02 
-7.0320O410E-03 
6.56585o63fc-03 
7.982368056-03 
-5.5435C7llfc-C3 
9.357121136-05 
7.009033C5E-04 
3.711d2105fc-O3 
-I .b684l6'5b-G3 


54 



i»i£: I (jhT i Ntj FACTCj.-<S 


DD&l 2>» 16 REAL PARTS Hk^T 


2. 61 707 6iki:-01 
-2*6l7<«okj:i7t + CC 
-1, i9‘f0l6t):>L-J2 
-2 . 0?0 yul VI t-0 I 
-y . iVttUyuJcL- J3 
-ij . 2o 7 V2 6ou t -0 1 


-6, Vd76b^66c-D2 

2. 51091 1316*^00 
2.2'^A252AlE-02 
l.l9l03aAlL-Ol 
7.S7^» aA?2 9L-0L 


-1 « y ytit^Ou 

720t*6T^>t-Ci. 
-6-o6i^/tt21t+0u 
- 5.ui2 Vt>ol oc-0^ 
“0*o>yl’»Vi 7E“01 


2. 6101^2 6 94E <-00 
^-6^1200 S^E*-C0 
2-743D29tt^e+00 
2-301 d^o 71E-01 
7*/7i26568E-'01 


1. 13301 Ob t?E- 01 
-9.G01l4202fc*0? 
-I.l4<»lv5l0t-01 
-2*6303/34CE-Ol 
4,6443941 lfc-03 


L. 5215951 3t-0l 
6*70 B6S273E + 0C 
7. 9144102 8E-02 
4.34366D0OE-01 
-3.226749026-02 


LEAST SgoAftt kcS luUAl.;», MHOF 1 ?) r N>tT REAL PARTS 6i»vSlt 11^ LUNTkOL POINT SEOOENCfc 


o* i3T99-*24L- 
1 .3a73'*d60c-31 
O.0l44Jly*tL-O2 
4.224439O2L-0- 
4. 59237344b-Cl 
-o.b27s52j2E-01 
2 .420653^76-01 
1 . 40 i 52 62 5L -Ui 
-6 . 507bo32 1 t-9 ' 
-6 . 19 7002aot -0 
-2.l65i34l3t-Ci 
-3.62 u 7052U-02 
-2.29i>2o5l 76-02 
-9 .60313 Obo c-03 
2.2b011324t-02 
5 . VS 7t04> oL-03 


-2. Hd536555b-01 
-4.324691966-02 
-7.016255726-02 
3.9497B103E-01- 
-4.69215!l4e-01 
-2.07^761086-02 
2.9b60e«98£-C2 
I, 70454800E-03 
4.53138240E-32 

l,0a005265b-02 

-2.^1''I4374F-02 

7.*‘63486a2c-03 

1.56475751E-02 

-1.353023986-03 

3.2289C939E-J2 

2.31815899E-04 


- 1. 94 654 46-01 

3 .»i>3»*>67t-Ci 
-2.4910692/t-Ol 

- 5 .>5 759 Is it- 01 
-4* o3364t22 46-01 
-4.422S44l9fc-01 

— 1. SSu6645ot — 01 
1.0?45493 3fc-Ul 
I.i27j0970t-0t 
1. 01o4g9Ci96-C2 
6- 422V465lt”03 
1, 6l4bb23 6C-0^ 

-3 • JE-Li 

-3.401 Vss 2 >c -02 
-6.29s9s 93 1 h-03 
6.4b4o0026t-03 


l.lt>o/4966E-01 
1.3JV04155B-C1 
4.77t>ld468E-0l 
s.95l4l354b-Cl 
1.0600330864-00 
6. Z446S922E-02 
-b. 751878596-02 
3.419377676-02 
-3.052715656-02 
-l.o73l4537E-02 
-2.2934L076E-02 
-2.42871026E-03 
4.626640246-03 
7.383315816-03 
3.58O04214E-03 
3.661044376-03 


-2.059<.3o526-0l 

3.74l2Vol3E-01 

-2.045346O7E-01 

-5.92450ci92E-01 

3.7402i<056-C2 

-1.209668396-01 

1.240350266-01 

-6.94513130E-02 

-1.O2177879E-02 

1.607868756-02 

5.6lo520o0t-03 

1.190203406-02 

-b.o9o47798E-03 

-9.13 ^o 0560E-03 

7.O09686V5E-03 

-5.76a9l345E-02 


1-4 1 73o2066-02 
-3. 40540 2276-01 
5.272540966-01 
-1.901933V2b-Dl 
7.081445416-02 
1.55 739920E-01 
1.878578566-01 
-1. 71956651L-01 
4.499282036-03 
1.230264946-02 
2.134798736-02 
2.493630466-02 
-0.O734O85OE-O4 
1 -98964190E-C4 
S.492107826-C5 
-2.346950186-02 


WElbMTlNs FACTuks aj<w* ‘*0061 3lr 16 REAL PARTS FIRST 


5. 52o305vac-02 
4.5O921Z4OC-01 
-1.474lLoo2t-0t 
-2.2939055OC-01 
5 • a 7d59 1 776-03 
2 .4098103^6-05 


-1.53C62741E-01 
-1. J9^91256E+0C 
"..93971364E-01 
-Z,622665''9E-0L 
1.537371 30E-02 
5.66136587E-01 


-2.502i>oli«t*00 
- 1 .80^45 o 03E-0^ 
-1.00 1405136 -01 
-2 *8i 7>707«6-06 
8 . 9 524 ^9546-02 


3-787400236+00 
4.20232532C-01 
-5.069871 796+00 
7. 209365626-05 
-4.145197736-01 


3. 7508lo9l£-02 
— 5.1 8bo3o3 4 E +00 
-2.43199V456-04 
2.80471049E-01 
1.20o01il6E-03 


-5.686O3904E-02 
6. l30olCl4£+0G 
2.41287364E-01 
1.03055026E-01 
1.912499C9E-03 


LEAST SUUAkE KtSluoALSf 'ICOLI *> ) ♦ NXT REAL PARTS hlKSI* IN LuNTkOL POINT SEQUENCE 


-8.o7d7531dE-0^ 
8 . 3£.a770326-0^ 
o.6220Vl99t-02 
9. 572 7oO/5t-02 
2.12037 87ifc-0l 
-4. 7002797ob-0l 
o.o726u33aL-0i 
-i.42 /O0712L-02 
-2.333B9947L-05 
- 7.44o43122b-u:; 
-5.1io4925lt-03 
-1 . 24l8o31it-L2 
-8. 7433043{»t-C3 
-1 .8080452/6-02 
o*4234ol d9b -03 
2 . 348910576-03 


-2.O018297OE-O2 
-1.2n67a25l£-<'2 
2.516647666-01 
1.76277243E-01 
-6.370276176-02 
1.293761966-01 
-2.1 599^445t-0l 
1. 649430356-03 
1.483339156-02 
-1. 136126996-02 
-3. 51'' 347926-04 
4.22952793E-03 
1 .66967563t-»02 
-4.680185436-02 
!.5457a8''0e-02 
-8.78343193E-04 


-7.7o5X377ot-y2 

-I.s742ti21>t-Ci 

- 1 , 432 7 7 1-01 

I.3567o5to7b-01 

- 2. JoU7 Js54t-01 
2.43ol742Q£-0i 

- 1.5522>233t-02 
7.733364/56-02 
1.3531i2o> 96-02 

-1.225l<>22b£-02 
1 . 189529196-03 
l*o90424£>5t-02 

-7.5044 74356-03 
5.4312277/6-03 

-2.7307 77oot-02 
1.754177C9 c-02 


3.80751301E-02 

-4.151655186-02 

5.815684316-03 

9.459731006-03 

-2.123155296-01 

-2.952553576-03 

3.402oo010E-01 

-1.194456486-01 

1.048188506-04 

4.397442556-03 

-4.31643512E-03 

1.359481556-02 

-2.320316C76-02 

2.042232116-02 

-2.14971466E-02 

1.220^21376-02 


-9.0o0834l 16-02 
-4.50320906E-02 
1.3176568 66-01 
-1.995057046-01 
1.387598576-01 
-0. 040352316-01 
4.74iol502E-01 
-2.64187855E-01 
-5.392bU326-04 
-1.9228B58 9E-03 
1.01182209E-02 
1.058310846-02 
1.3733704 DE-02 
-2.720209516-03 
2.759963996-02 
-2.002517356-02 


-1.589522806-01 

5.079735036-02 

1.82658499E-G1 

1.191079606-01 

-2.208548456-01 

-1.34336148E-02 

-3.73290O52E-01 

1.78009211E-01 

2.214567806-03 

-7.45544516E-03 

1.11714225E-02 

-1.22944136E-02 

2.04O08347E-0Z 

-3.173030186-02 

4.003801656-02 

-1.79669131E-02 


PKDGBAM D2182 ... CCMPUTE GENERAL 126U AEKUUYNAMIC FORCES AND SOLVE FLUTTER DETERMINANTS 


hT-7 FLUTTER «OOEL» 4d CPT5» 3 MUUtb, 6/22/ /3 


ECHO UF NAMElIsI lP21o2 FOLLOWS 
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NSP S T 
NCl>ST 
MASS 

LiMEGrt 

G 

.'^ALHHA 

ALPHA 

YSP5 
i SPS 
MIC 


= 11 , 


= 10 , 


0.3 d/ 0^L-03» 0, 

O.Q, O.Gt ^-0, 
0.0, 0-0, 0.?, 
0.0, C • C t , 


I C.O, 0*C, 0,Ct O.Of O.G, C*S17faE“'03r O.Ot 
0,0, 0.0, C ,2635<it“03 f 0.0, U.U, O.U, 0.0, 

o.r, o.o» 0 . 0 , 0 . 0 , 0 . 0 , o.c, o.o, o.c, 

O.Off 0*0, 0.0, 


:.io^3L + 03, 0.725E^03, i, 1, 1, 

0.7xi Dh<-C3, 

O.'Jt o*u, 0*0, T, I, I, 

2U 


= 1 , 


= 1, 


1 lfc-u<_, 
0 * 2 oL'' C 2 , 

1 , 1 • 

I , i , 


0.6E~'^3, '".7E-03, 0*9t-03, L.it-02, 

u.l2E“02, 0.13E-02, 0.14t-*02, 0.i:>E--02, 0.ibE-*-02r 
"i.lOE-V, ^.19E-02t 0*2t-02, 0.22t-02, 0*2^b-02, 
J.20E-"'?, 0.3E-02, I, I, I, I, 


i 


t I, If I, Ip I, If 
i» 1* I, I, I, I, I, I, 


-0.2ldL+CUf 

O.IS 

'•0 . i 9 y b +00 , 
J* 1 ?9b+ Cu , 
-:'.io9t+ro, 

0.L^Ob + Luf 

• i 6 + C 0 , 
0.2:>L»£ + C0, 
~'v ♦ / d b i , 
0.29C+0 L t 

0 .93^t+C0, 
0* b2 vt'ol, 
2.^7lii + c:, 
0*123fc^^;U, 
0* li^Obb^’L 1 , 
0. llbL + OOf 
0."t*+bt+CCf 
0. 122t + C0, 
0. ^ouc + CO, 
0. ldOE + 00, 
0* L>2Si: + Co, 
1, i, I, I, 


-0.192E + OC, -0.l3tib + C0, '-O.llt + 00, -l.97l-01, 0.36E-01, 

■^.25?E+rC, 0.385E + 00, u*b2i»E + 0Ut Q.b7Sb + 00, 

Ll9fc + 00, -0.:>j7t-Ci, -0.2c»b'02, 0.79b-01, 

0-A17E+0C, 0.i:>Gb + UU, 0.70bt + 00, 

-?.79E->^1, -^.73c-02r O.+t-Oi, 0.1^£ + 00, 
0.9i>lE + 00, 0.:>9LE + OG, Q.73 7b+00, 

'^*37E-01, 0.3ib-01, 0. /9E-CI, O.lS9E + 0Cf 


0* 29E+C0, 

127E + '''', 

O*3'^7E + 0C, 

.83b-01, 

^.36"'E + 0C, 0*^,866 + 00, 0.o*::3£ + 00, 0*7o7t+00, 

-C*3AF-^1, ^.4lb'02, D*6^fc-0-i, O.U7b+00, 0*l96E+C0, 

0.^98E + ‘^r, 0.921E + 00, O.obi>t+Ou, 0.79oE+00, -C.27b-01, 

O.AAF-01 , 0. 936-31, 0.15^6,^00, J.233E+0U, 0.32bE+00, 


I, 
1 . 

i, 
1 . 

J 


i , 

i . 

I , 


I f 

i , 

i t 

0, 


i , 

I , 

1 , 

i * It 

3 b C C • 
O.o39b+C3, 
0*:;09b + C3, 
0. o2 It + LU , 
0 * 2 o ^ t + O t 
0*<,7j?6 + 00, 
0-2ol6+C0, 
0 * 3 o L C 0 , 
0.22dt+CO, 
3 *^^^t + 00, 
3* 131b+j3» 
O*202b+C0, 


^.5576+ 

0. 131E+CG, 
*593E + "'0, 
3*1 76E+C0, 
''•637E + '"?, 
0. 168E + <^C, 
0.5A6E+^^t 
3.169E+^r, 

*58AF+'"0, 
3.23E+00, 
3.621E+0C, 
I, I, I, 
I, I, I, 
I, I, If 
If I, I , 


0.68b£+00, 
J-1966+00, 
0*7 I6fc+00, 
3*239t+00, 
0.748E+00, 
3*225E+00, 
0.6o2b+G0, 
0.2 186+00, 
3*o98E+00, 
-^•aTlE+OO, 
0*735E+00, 
I , I , 1 , 

T , i » 1 , 

I, It 1, 

It 1 , 1 , I 


IFi-?!, ^.721E+00, 

0.328t+’^^, 0.3286 + CC, 


B22F + 0C, 
*33aE+<^^, 
0*64E+00, 
.2>7E+ "C, 
0.605E+0C, 
*259E + ^<", 
O.A59E+^C, 
^*272 E + "'0 , 
'^*^16E + ^C, 

0. 3346 + 00, 


0*Sb7E+00, 
3.2a^E+?0, 
0*4646+00, 
3*2686 + '^'^, 
■3*45SE + 00, 
0.26lt+00, 
0* 591E + 00, 
3*2456+30, 
0*44 7E + 0C, 
0. 149E+0C, 
0*3C9fc+00, 


G.O4.i#t + 00, 0.406-02, 0.446-31, 

U*27j»c+ 00, O*i>o5t+00, 

0.o>o6+00, 0.3It-0i, C.776'01, 

0.31:>b + 00, 0.^05t + 00i 

0*oo9b+jCt O.o7t-Ol, 

0.2o6t+00, 0.3 60 6 +00, 

0.7o;j 6+00, 0.9ii:>t + G0, 

0.2736+00, 0. 33it + 0 3, n.<,£+00, 

0.o^.lb + uU, 0.9o26+00, 

^•32b+0o, C*37oc: + 00, 0.4^,26 + 30, 

u.ooL+00, J.lb+Jl, 1, 1, I, i, 
It 1, 1, I, if i, i, 1, 

I* 1, 1, 1, i, 1, I, I, 

J, -i, I, 1, i, i, i, 1, 

I, ii If I, 1, 1, 1, 1, 

545L+0G, 0.^:ii»t + 0Ct 0.3odE + 00t 
C.3o:^c + Gu, 0*:>04b + CO, 

U.429E +Uu , 0.355E+0 0, 

0.3L+00, 0.3O4L+00, 0.^88E+0C, 

0.33o£,-00, U.2dbfc+00, 

0.<:736 + C0, 0.j4:>b + 00, 

O.355b + 00, 0.<id2 £ + 00 , 

a.2S<,b + G0, 0*25 7t + 00, 

0.2bt>6+00 , 0. ^49t+O0 , 

0.^:p*tb+00, 0.246t+00, 

0.:>7bL+u0f 0.13bt + 00, 

0*lobt+0J, 0*18c6+00, 

3.^j>ob + G0, 0.:j3db + QCi 


if U 
if If 
if I , 

it If 


It 

It 

If 

If 
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'0.^6E-02t -0.dG-02» O.ob-U2* 0.34E-Qit 

+ "'-218E^"''f 3.303E + OOi 0*^27t^-00f 0. :>3bt <-00 1 

-*0*iuL£ + Ou» -O.T47F^CCf -J.178fc + 00f -O.iVJt + uO, -O.lobt^OQf 

.l‘t:>b + CO» C.lS2t*OOf 0,30--fc + 00# Cl* ♦OO t 

O.^l^L^LO. -0.21’?*E + CCf -'0*278£+00f -0*323t+C0t -0.347E+00# 

-0* 33 /E + C ■:, -:',28? F+'^Ot --).L27E+?0» 0 . 1 i> 00 1 O.^^aE+'OOf 

0*ti2c + C0, ^•^j'^6E + 00t -0*?3oE+0Ct -U.3b/E+OOr *0.H23b<-00f 

“3*^3ib<-s:0. -?.4^4E+"«Ct *-0.343b+00* -u.2it:t00f O.lCoE+00* 

0*2'7'iL^OOt ?.4?6E^'^Ct O*5b + C0# - 0 . 00 r -O.46^L-^00» -0.^b7t+00* 

- -:»,418E^^C, -0.351b+00t -0.27it + 00, O.bob-Uif 
O.^b^tL + COf '".4E + 00, 0.‘t97E-t-CCf It it I» I» it i» 1» I* 

it If If if If tf If I» It If If I* It It -It It !• I* I* I' 

i, 1. If If A, If I, Tf If 1# If If It If If if It It 1# If it It It 

i« if if If if If If If If It If If If If It if if if It If if If it 

If If If i, i, If If If If It It if If it it At It It 0.l22fcA-0Ci 

0. 33c“3Ii -0.5?E-Clf -0.l6oE+CCt - C . 33 :jL OG t -0*4:>4t*00f -0 • 3Ef OC ♦ 

-C.^9db + 0jf -0.?75E + CCf -0.7bE'-Cl» U.^j^t + OUt O.ob£-01f 0.77b-02f 

-Otot-Glf -0-191E+'^?f -'^.32dE'i-C Ct -0.^2E + O0f - 0* 4:>i Lt^CO# -0*4lE+0Or 
-0*<:^2t*^0t 0*66b-*0Lf 0.:>b7E^00f 0.4b»t-uif -u*2:)t-0i» -U*luabt-00f 

-'^.321E^^^i -O.3bi»E^'?0f -0*:597b + 00f -0*3HiE^0Oi 

^c.2llb+LJf 0.198Et-CCf 0*63t)b^-00f u.Iut-Q2f -0*o4t-01, -0*136bt-00f 

-0*227c + C}f .3lIE + '''"f ^3.3:>lE'i''^0i -C*347£*00f “0*Z9b£+0Cf 

-u* Lii9ct-Uuf O^?X2Et'O0f 0*6dE-*-0Cf -O.^ob-OIt -O.lOBb+'OOt -O-loiE+OOf 
-: .23 jbfO;>f -':.298Et-^0, -).321E + 00, -0.3Cot^^00» -0*^OEt00f -0-34t-01t 
G.403Lt'C:t r.719E + ^0, -0.d4E-Clf -C.i3.::c + C0f -0.i74t + 00i 
-0.23obt jJf -3.234E + ?r, -0*298E + OQt -Q*27^tt-0Qf -0*l9b£ + j0f a*Ut-OU 
G.^T'fctCOf '^* 757E + ':0f -0.i0!>E+^0f -O.lHot+OOf -0. t 

-'^*269E + CC» -0.27bE<’00f -O.Z^tbfrUUt -w*od£-^Oit 0.2a9t+00f 

0. :i3:>bA'J0f '^.798E<-'''^f -)*114E + 0Cf -0.i:>lcfC0i ’-0* i S*ibt^OO » 

-w^.233tf'uu, -0.256F+C0f -.)*256e + 00f -U.19t<-00t 0*^bc-Ulf 0*337fc + 00i 

J.i^SrJtf'QO, ’'*84?F + '^^f -0.113E + 00» -Q*l:^4tf00f -0.19bt+0G# 
-0.232£f00i -'^.24'^Et-CCf -0* 2 2 3E*i' 00 t -u.UE*-00f 0.io3t+00# 

O.^+^ubfLJf '^.6b9E + rCf 0.d94E+O0f -0 . L^. Tb + OO t - 0* lo b bt-0 Ct 

-J.20:^Et-o0» -^*?32EA-C0f -0-2il£+00t -0.i7bt*00f 0- Qt 0*29:^t^00f 

3.3j.c^uLf J. 727Et*0Cf 0.946Et'pCt -0.i<tL^-UUt --U.loLt-OOf -0 * 2 13 E+ J Ot 

-0.^3iic*-^Df -"'.219E+'"'^t -"'.IZoE't'OOf 0*iA.**t*G0f C.3bbt*'0Ct O.bE + OOt 
0*/9b£+00f O.lb+OLf It If If It It It it I» it it it it It It it It 
I, it If if Tf If If If If If If It If i# if It if It if if If If 

It i, I, I, 1, If If If If If It If If if If if it It it if If It It 

it A, It if if If If If If It If If If If if it if it it it if If If 

It if If If it If If If It It If If If if it If if if It if it it if 

1, i, i, it it If If If It It I* I* It I* It At it if It If If If It 

It If if If if If If If If It It I* It it It At if if it it it It if 

it it 1 * I . If I . If If If if It It If If I f 1 f 1 f If If it It It If 

1, it it if if If If It If It It If If it if If it it It it If It It 

i, I, i, it If I* If If If If It If If At if if if It it if if If It 

it 1, 1, it it If It If If If It It If if if it if if it it if It If 

i, It If If if If If I* If It If If if if If At if If if if if It It 

It It It it At If If If If If If If If it it If if if if if If If It 

it It If If if I* It I* It It I* It It I' I> I* At If If If if If If 

1, It If If If It If It If If If It it if it it if if if If It if If 

it It it if if It If If If If If It If if if if it If It if if if If 

It If If it if If If If It It It It If if It If if I* If If It If If 

it At If if If If If If It if I* I* It if it if If it it if If if If 

i, it it It if It It If If I# If It It it It If It If If If it If It 

it If If If If If It If it If If It if It if if If if If if It If 

it It if If if If If If It If If It if if It if If it it it it It It 

1, it it If At If I* It If If If If If if it if it if if If If If If 

it It if if if If If It if If If If if If If If If It if It If It 

I, it It it At If If If If It If If If it it If if If it I# if It If 

it itlf It If If If It If It I* It It it It If if it If It it Ii If 

1 , I, I, it If If If It If If If if It it it If if if if if if if If 

A, i. It If It If If If If If If If if it It if if If If if If if 

i. It it it if If If ^t I* it It It If if if if AI If If if it It If 

It It if if If If If If It If It If If if it it If if it I* It it 

It If If it if If If If If If If I* It If it If if If if If If If It 
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CWi>PST 
CwSP I 
CrtCJI 
AlAd 

CWPK2 

C rt P rv 3 

CWPP .4 

C^Pkb 

CrtPk o 

C k b 

CwPRbX 

aCBmK 

AiC 

MNP 

i.>iTaR 

NPNCHvj 

JIJ 


MACH 

K 

TANl E 

TANLTt 

S 

NMjUS 

NSYM 

Sbt^iU 


It It it it It 

if I* it it I* 

if it If it if 

1 . It It It It 


Ot 
0 1 


It If 1 
It If I 
!t It I 

It It 


It It It I . 

It If It If 
It I t It If 


[« If If If if 
it if if if if 
if it If it if 


it it if If I 
if if If if 1 
« if If 1 1 1 


f It 
f It 
t It 


0*ut 0#2F + 00t 0.3E + 00f O-^fc+OOf U-^t^-OOf C#oE + 00f 

J.intoo. o.iE^oif i, i, if i, i, i, i, i, i. 


Ot 

Ot 

Of 

Ot 

Ot 

Ot 


■» it if If if If If If 

it It If If if i, I, If 1 , I, 1 , 
ii » 

2 . 

it 

Of 


( 

i 

i 

i 

{ 

0 . Ot 

J. Oj ♦ 

( "‘-'"t 

?. 0 I » 

( 

0 * 0 t 

0 . 0 ), 

I 

C • U 1 

0 * 0 i f 

i 

O.Or 

u.O) 

0 * 0 t 

0 .OJ t 

( ^•"'t 

'". 0 ) t 

( 

0 -% 

C.O), 

( 

0 • 0 f 

O.Ci • 

i 

0 . 0 , 

O.OJ 

0 • u t 

j.Ci t 

( O.Of 

0 * 0 ) , 

( 

O.Ct 

0 . 0 ), 

( 

U*U t 

U.Oi f 

( 

0 • U f 

0 . 0 } 

0 • »J f 

• 0 i f 

( ''•^f 

t 

1 

">.Cf 

0 . 0 *. 

i 

0 • C f 

C.OJ f 

( 

O.Of 

0 . 0 ) 

0 « 0 f 

0 • 0 ) t 

( 0 , 0 t 

O-r*) , 

i 

o.c, 

0 . 0 ) , 

i 

0.0 « 

0 « 0 ) f 

{ 

0.0 f 

0 . 0 ) 

1 

0 • 0 t 

0 » O' i t 

< o.'^t 

C.O) f 

i 

O.Ct 

0 . 0 ), 

i 

u *U t 

U.Oi f 

i 

0 « u f 

0 . 0 ) 

( 

I 

O.Ot 

0 » 0 t 

0 .0 J , 
C.O) t 



( 

O.Ot 

''. 0 ), 

i 

w/ • 0 f 

0 « wi f 

( 

O.Of 

0 . 0 ) 


0. 12t + 0lt 
L • t -I- 0 c » 

0 • 1 1 b c + 0 1 1 

0 ^ 352 b:>L>G 0 , 

O.blE+oOt 
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HACK* i-2000000Ufc+00 


K= A.ooooooooe-ci 


TANLt-* 1.2I600000E^OO 


TANLTE^r 3,525i0000£-0l 


S» 8,10000000£-01 


MMODS- 3 NSYM" 1 


INPUT UtfLECTIONS H SUB 1, IN C8DEH L.t- TU T-t-t J 

-2. IbOOOOOOE-Ul -I, 99000000E-01 -L.69U000GUfc-Ul 

4-t)OOOOOOOt-03 3. lf'rCf'0?0K-02 6*70G0GO00E-O2 

-L,9200OOUOe-0i -1.63CCOOOOE-01 -1,2 lOOOOOOfc-Ol 

v*^000000at-02 T.70CCC009E-02 1-loOOOOOOE-Ol 

-1 .baOOOOOOfc-ai -l.l 9OCOOOOt-0l -7,90000000£-02 

a ,2900000Cit-02 1.23000000E-01 I •abOOOOUOE->gi 

-I . lOOOOOOOt-Oi '5.57C0'^0^OE-02 -7-300DCOOO£-03 

l.31JODOUJk-Oi l.?6nrcC00E-Cl 2*25000000E-01 

-9. 7COOOOOOfc-02 -2 ,6?OOOOOOE-0 3 4 .OOOOOOOOE-02 

1 .96000000£-0i 2.39CCnOCOE-01 2* ojOOOOOOfc-Ol 

3 ,6O0Ok)OOUfc-Oc 7.9OOCCOOOE-02 1 .^OOUUUOUt-Ol 

2. 7300:000fc-Cl 3.15Tr''00'3b-Cl 3*boOOOOOO£-Ol 

I .37UUOOOOC-01 I . 79^>CCO?3E-01 2. l&OOOOOOt-Ol 

3*g5000000E-01 4 . 95':)'>OOC9E-01 <*.4o000000L-0 I 

2.&2OOQOO0b-01 2.900CCC00t-0l 3*27000000fc-0l 

4.7lG0000Jfc-0Z 5.O8OCO00OE-O1 5.4bUU0UU0t-Ol 

3<, BiOOOODOl-01 4 . ! 7'^'J'*0'''>E-ni 4 .slOCfOOOOE-Oi. 

5 .92000000C-01 6.27O0CO0OE-OI 6. o2000000E-0l 

r».2bOOOOOOt-Ol 5*50'):OOOOE-Ol 5,9l000000t-0 1 

7.i630000Jt-DA 7.4BT'f'0»''5E-Ol 7. b30DOOOOt-Ol 

0. 7bOOOUUOE-Oi 7<.06OCCO0OE-Ol 7.37JJOOO Jt-01 

b.^tiOOCDOOc-Ol 8.B90''f'0'>'^E-Cl 9.2>00')000t-0 i 

1 . LCJOJUOOt+DO 8.2POOOOOOE-OI 6.^OO0UOOOt-Ol 

1 ,20000000fc-0j -I.^IC'^OOOOE-Ol -2.130DuiCC-Ob-Cl 

?4210000u0£-0i e».87'>COCODt-01 4.o4000000c-01 

-^.oOOOOOOOt-O^ -l.47^00000E-01 *2 • /boOOOOO£-Ul 

b.4*>OOODOOL-0 1 4,?9n-rr^0E-0l 3 -3bOOOCOOE-Ol 

-b.OOUOaUUOt-U3 -1.78000000e-0l -3.2!>UU\jU0Gt~0l 

4. SbOOOOOOt-Oi 3.550^00''-)E-?1 2.bbQOOOOOt-Oi 

-O400000000fc-03 “l.97ncOOOOe-OI -3<.4700000t>t-Dl 

j.abOOOOODt-Ol 3.090'^i:C!:aE-01 2.or>UUUuGijE-Ol 

o.OOnODOOac-C3 -3.370OCOOOb-Ol 

3.f)QOOUUOuc-Ol 248a0CO0C06-Ul 2. b 7D00000c-01 

3.4CDOJOOOL-02 -I .4 5'^0CCf'?6-0 1 -2.O3000000L-01 

3*28000000L-Ul 2.840000CCE-01 2 .3b000000£-0l 

U09s)00030t“0l 'l.?0^?0000£-02 -1 .2 70Co0J0b-Ui 

3*2800000 3 L-Ol 3 1 2* 7i»00D00Ob-01 

2* lti30O0O’Jt'’0l 1,920COOOO£-01 I •bbOoOGUOt-0 I 

3. dbOOOOOOL-01 3.64''''?0onfc-0l 3 .4b0000aDE-01 

3.Ob0QOOOCJ<.-Ol 3*O20rconCE-0l 2 . 9«OOOOODb-0 1 

b.04DODOOOt-31 4,880''CO00e-0l 4. /3U0uUUUt-C» I 

4*2700000C^c- 01 ,r,i 4* 1 200 uOO 0£-0 1 

6*39000000l*u 1 6.210C0000E-&1 6* ObOOOOOOt-Ol 

b43bJOOOOOL-Oi S.IVO'^OCOOE-Ol 5 .06000000£-01 

l.220000COc-01 a,8O0C'C'OCOe-O2 4. ■>» 00 LI 00 00c -02 

-1 .O^OOOOOOfc'Ol -T. 1400COOQt-Oi -1 .lbOOOOOO£-Oi 

3.baOOODOLc-Di 7.7"‘''''r'5-i-^e**0 3 -2.b000c000£-02 

-1 .4bOaOOJOb-Ol -1, blOOCOOO£-Ol - I .p-tOwUOLlUt-Ol 

-3.300000001-02 -a.''OOrrnoOE-02 -l-QttOOOOOOfc-C I 

-1 .d^OOOOOOE-Ol -i. 94COr>OOOE-Dl -1 . 4o000000t-01 

-l.obOCOOOOL-Oi -l,9ir0C0^''e-01 -2.l300v000t-0i 

-2.340OQO00b-0x -2. 3 3'^f'^'' ''^t-0 1 -2.3iOCOOOOt-01 

-3*3>OOOOOOc-Ol -3* ?80COOOOt-Ol -3 -2i0000lju£-ul 

-2.o9)0j000t-0 1 -2.5b''^<''''^OE-Ol - 2 .43 JCOOOOt-C 1 

-4.340000001-01 -4.20CCC0‘^OE-0 I - 3 . dbOOOOO Ofc-01 

-2.7b000030t;-01 -?.5b0 0n!:o'>L-Cl -2 .2iOC-OOOOt-Ol 

-3.0b 000 000c-0l -4.blO''0 0''''E-0 1 -3. 97DOOOOOt-Cl 

-2*400000ugt;-Lll -1.90000000E-J1 -1 .1 JJUOuOUt-01 

-4 .9b000000t-3i -4. 1 ^'0^'E -01 - 3 .41 OOOOCOt-C i 

-d.dOOOOOOut-J2 4. aOOOOOOOt-02 l.a3000000E-Ol 

-3. 750000GOt-01 -2 . 9?0o^'n0''E-0 1 -2 .llOOOOOOE-01 

2.09000000E-0 1 3, *» 7CCCCCQE-0 1 4.4u000000£-0i 

- bCOOOOOOb-Oi. 646''000<'00e-0 2 1.9o0u000o£-ui 

3.3bOOODOOL-0 1 !>495''0''f’??E-'' I b, bVOOOOOOE-Ol 

b.3400UOjOt-Ol 5. 07OCCOOOC-O I 6.3t>000000£-0l 

?.9d000000t-0i 8.4‘»'‘"CCOOt-01 8. 4-*0Cc wOOE -0 1 


:»PAj\ 8TaI iUNS 

0. d.l-^^OOOOOC-C2 1.620cu000£-01 

4.d6000000t-0i b.i.TCCOt'OOE-Ol 644bOOCOOOfc-01 


NbUARO E1»ST« MODE 1 FIRST 

-1.24000000E-01 -T.20000000E-02 -2. 700000CCt-02 

i.22000000fc-01 l.ddODOOOO£-01 

-tt.3ODDO000E-02 -3 440000000E-02 9.00000000E-03 

1. b9000000£-0l 2.30000000E-01 

-3.70000000b-02 4 .lOOOODOOE-03 4.40000000E-02 

2-iaOOOOOOE-Ol 2.71000000£-0l 

3.100000COE-02 b. 200000006-02 9.30000000E-02 

2. 730000 OOb-Ol 3*2000000 0^01 

7-90000000t-02 1. 170000006-01 1. 550000006-01 

3.310000006-01 3-7toOOOOOO€-Ol 

I. 69000000E-01 l.9bOOOOOOE-Ol 2. 33CODOCCfc-Ol 

4-OOOOOOQOb-Ol 4.420000006-01 

2.550000006-01 2.900000006-01 3. 2bOOCOOCE- 0 1 

4.8O000000E-01 5.25000000E-01 

3.O30000006-01 3.9BOQOOOOE-01 4 . 34000000E-D1 

5.840000006-01 6.2I000000E-C1 

4.800000006-01 5.21000000E-01 5.570000006-01 

6.9S0OOOOOE-O1 7.35000000E-01 

b.23000000E-Ol 6.550000006-01 6. 86000000E-01 

B. 210000006-01 b.oOOOOOOOE-Ol 

7.670000006-01 7.96000000E-01 8.250000006-01 

9.62000000E-01 1. OOOOOOOOE+OO 

4.590000006- 01 2.dbOOOOOOE-Ol l.3eOOOOCOfc-Ol 

-3.360000C06-01 -4.O60C0000E-01 

3.53000000fc-0l 2.49000000E-01 1. 310000COE-01 

-3.870000006-01 -4.820000006-01 

2. b20000CQ6-01 2.280000006-01 1.36000000E-01 

-4.23000000E-01 -4.870000006-01 

2.61000000E-01 2.33000000E-01 1.49000000£-01 

-4.310000006-01 -4.O80030006-01 

2.590000006- 01 2. 45000000E- 01 1.85000000E-01 

-4.040000006-01 -4.180000006-01 

2.O10000006-01 2.54000000E-C1 1.820000006-01 

-3.430000C0E-01 -3. 510000006-01 

2.54000000E-01 2. 460000006-01 2. 020000006-01 

-2.10000000E-01 -2.710000006-01 

2.57000000E-01 2 .45000000E-0 I 2.340000006-01 

l.ObOOOOOOE-01 5.80000000E-02 

3.280000C06-01 3. loOCOOOOE-Ol 3.0900000ut-0l 

2.9100'>000E-01 2.84000000E-01 

4.590000006- 01 4.470000006-01 4. i b COOOOOE-0 I 

4.0bOOQOOOE-OL 4.000000006-01 

5.91000000E-01 5.7b000000E-01 5.58C000C0E-01 

5.00000000E-01 4*970000006-01 

l.b0000000£-03 -4 4o00L-0300t-02 -8. 400000006-02 

-U2/000000E-01 -1.400000006-01 

-b.4000?DOOE-02 -l*03000000fc-0l -I .32000000fc-0l 

-1.850000006-01 -1.800000006-01 

-1.560000C0E-01 -i.blOOODOOE-Cl -1. 790C00006-01 

-2.050000006-01 -c . 1 3O0D00OE-01 

-2. 270000006-01 -2.350000006-01 -2.360000006-01 

-2.320000006-01 -2 .320003006-01 

-3.UOUOOQOE-01 -2.980000006-01 -2. 84 COOO CCE-01 

-2.310000CCt-Cl -2.19000000L-01 

-3.510000006-01 -3.210000006-01 -2.980000006-01 

-l« 7bUOOOOOE-01 -l*2bOOOOOOE-01 

-3.4700COC06-C1 -3.0o000D00fc-01 -2. 7400C0CDE-01 

J. 1.24000000E-01 

-2.980000C06-01 -2.600030006-01 -I .950000006-01 

2.950000006-01 3.ba000000£-01 

-1.290000006-01 -3.4000000CE-02 8.OOC00O00E-O2 

5.3OCOOOC06-C1 o. 000000006-01 

3.120000006-01 4.O3OOJO0OE-01 4. 74000000E-01 

7.270000006-01 7.9800QOOCE-C1 

o. 800000006-01 7.190000006-01 7.57CC000D6-01 

9.4b0O000C6-01 1. 000000006*00 


2.430000C06-01 5.24000300E-01 4.050000006-01 

7.290000006-01 O.10C000006-O1 
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ANK MATp\lAr i>V MHOAL COLUMNS OK 
HtAL iMAr, 

^7,6006ouJ3L-0i ?. ■*806666 ^£-02 

-V. <JC03o32vt-01 ! ♦! ?7'59670e'0l 

l.X947os31fc+00 *?.04O219C9fc--‘3 

-H. 67J41111L-\;1 -3. 923l2979fc-03 

l,3‘*7391d0L-0^ -«.'^'39 7!9«5£-'>4 

-%-9t>d^7474L-0X 5,3972069iL-'C2 

-1. 7 3341^3 :i6t*-G^ -?,''50581 91E-01 
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-4.61480535E-02 - 4* 58958959c-03 

r.ocoooocct +00 o. 
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-\FU J * IMtl J , lwi2l,erc. 

1 . 56«3i4olt4 :i -1. ' 06A66b2E*TJ 4 *l VoL‘ 1 3t*uU 

t tNE«.LyutNi,lES 

i .O3idoii27t + 0<. 3.53923l94t+'>2 6.-»4?07975t +02 

UAMPlfvb l^utPhlC lEWTS G 
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STIF* PAKAMETEPS 
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-1-726628786 <-00 
4,04677!6(5t'02 
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-9.00£>b033 7E-02 
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-4*1524b938E-01 
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tli^ENVECTLkl 2J 
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Results for most alpha values are not shown 


alphas 9. CO. 
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13.13.33. TlMfc LO attach 

ia.ij.33,tNu» ftTUh 
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Langley Research Center, 

National Aeronautics and Space Administration, 
Hampton, Va,, December 28, 1973. 


63 



APPENDIX 


EVALUATION OF THE SUPERSONIC KERNEL FUNCTION 

Based on equation (16c) of reference 3, the kernel function can be expressed in a 
separated form 


K(M,k,XQ,YQ) 


2K(M,k,XQ,yQ) 

(2b)2y^2v 


(Al) 


where Xq = x - y^ = y - tj, and the reduced kernel K is a dimensionless part of K. 
For the steady case, 


-Xr 




(A2) 


For Yq = 0 (that is, ?? = y) 

K^M,k,XQ,0^ = -e ^^^0 
The general expression for rj ^ y is 


(A3) 


K=-e-'^^0; 


XqC 


-i2kxo/^2 


'Xq 2 - /32y^2 


COS 


^(xo2 - 


1/2 




+ ik 


y 1 \ e 

“I f77 


-i2k|yo|T 


dr> 


(A4) 


where the upper limit + T 2 and lower limit = '^1 ~ "^2’ 


T, = 


^0 




M 


^2 = 


\Zxq 2 - g2y^2 

^>o| 


(A5) 
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APPENDIX - Continued 


In the integrand in equation (A4) the quantity 
for positive values of t as in reference 4 


+ t2 


is closely approximated 


■■ ■ ^ « 1 + + N26^^ + N3e^'^sin tit (A 6) 

]jl+T^ 

where Nj = -0.101; N2 = -0.899; N3 = -0.0932307; a = -0.329; b = -1.4067; c = -2.90; 
and only N3 varies slightly from the corresponding constant in reference 4. With this 
approximation substituted in equation (A4), the integral can be evaluated in closed form. 

The upper limit is always positive. The lower limit is sometimes positive and 
sometimes negative. Since the approximation of equation (A6) holds only for positive r, 
negative is accounted for in the following discussion. 

Let the integral be defined 



\/l + 


dr 


(A7) 


where Qj is complex, Qj = Re^Qj^+ ilm^Qj). Since there are four terms on the right- 
hand side of equation (A 6), 


Ql «Q2 +Q3 +Q4 +Q5 


(A8) 


and each is complex as follows: 


_ 1 


cos 2k|y^|T - i sin 2k|yQ|T| dr 


'0 


lLs2k|y„|T^Tcos|yo||r,|-Q| 


(A9) 


where upper of the two signs + and the upper of the two quantities 
positive and the lower of the two are used for negative t^. 



are used for 
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APPENDIX A — Continued 


^3 = ^1 e^^^cos 2k|yQ|T - i sin 2k|yQ|Tj 


dr 


where 


N4J cos 2k|y(j|T^ + 2k|yQ|sin 2k|yQ|T^j 


a T, 

- e /a cos 2k; 


+ i 


{* ^|yo||^i| + 2k|yo|sln 2k|yo|| T,|j1 

"(-a sin 2k|y(,|T„ + 2k|y^|co6 2k|y^|T„j 


a T, 

± e ' sin 2k| 


(a sin 2k|y„||T,|- 2k|yg|cos 2k|yj||T;|^ 



N4 = 


Ni 


+ 


etoi) 


Q4 = N2 J e^'’’^cos 2kjyQ|T - i sin 2kjyQ|Tj dr 

is obtained from the right-hand side of equation (AlO) by replacing by 
by b. 

^5 - ^3 irr^cos 2k|y^|T - i sin 2k|y^|Tj dr 

8<|e^^'^[c sin - Ng cos NgT^^J - sin Ng|T^| - Ng cos NejT^^ 

+ Ng|e sin N^t^, - N7 cos _ e ^ sin N^jT^I - N7 cos N7 |t^^ 

- Ng<^ «|c COS + N7 Sin N7TJ re I ^Ijc cos N7I | + N, sin N7 |t^^ - 


= N, 


(AlO) 

(All) 

N2 and a 


(A12) 
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APPENDIX - Concluded 


where 


N0 = 7T + 2k|yQ| N7 = 7T - 2k|yQ| 



and where the upper or lower of two signs and two quantities are used as described with 
equation (A9). 

In the subprogram KERNEL, the number of calculations of sines and cosines has 
been minimized to only four of each with simple arguments, and the sines and cosines of 
the sum and difference arguments in Qg are obtained as their ”multiply-adds.” 
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